Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 943,126 A→G 5.7% E723G (GAA→GGA)  dmsA → dimethyl sulfoxide reductase, anaerobic, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913943,1260AG 5.7% 214.9 / 2.5 70E723G (GAA→GGA) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base A (36/30);  new base G (2/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

TGCCGATCTACACGCCGGGCTTTGAAAGTTATCAGGATCCGCTGAACAAACAGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCC  >  NC_000913/942982‑943274
                                                                                                                                                |                                                                                                                                                    
tGCCGATCTACACGCCGGGCTTTGAAAGTTATCAGGATCCGCTGAACAAACAGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAAtg                                                                                                                                                  >  7:22522/1‑149 (MQ=255)
        tACACGCCGGGCTTTGAAAGTTATCAGGATCCGCTGAACAAACAGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGACGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCaa                                                                                                                                          >  7:299980/1‑149 (MQ=255)
                              aTCAGGATCCGCTGAACAAACAGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAAcgc                                                                                                                    >  5:130446/1‑149 (MQ=255)
                               tCAGGATCCGCTGAACAAACAGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGAt                                                                                                                                <  2:20303/136‑1 (MQ=255)
                               tCAGGATCCGCTGAACAAACAGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGAt                                                                                                                                >  1:20303/1‑136 (MQ=255)
                                   gATCCGCTGAACAAACAGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAAcgcg                                                                                                                   <  5:240633/145‑1 (MQ=255)
                                   gATCCGCTGAACAAACAGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAAcgcg                                                                                                                   >  6:240633/1‑145 (MQ=255)
                                   gATCCGCTGAACAAACAGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTAt                                                                                                               >  8:174399/1‑149 (MQ=255)
                                   gATCCGCTGAACAAACAGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTAt                                                                                                               >  4:68979/1‑149 (MQ=255)
                                          tGAACAAACAGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAAc                                                                                                        >  1:211823/1‑149 (MQ=255)
                                           gAACAAACAGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACg                                                                                                       >  1:248564/1‑149 (MQ=255)
                                           gAACAAACAGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACg                                                                                                       >  4:278281/1‑149 (MQ=255)
                                            aacaaacaGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACgg                                                                                                      >  3:300026/1‑149 (MQ=255)
                                               aaacaGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGa                                                                                                   >  1:295208/1‑149 (MQ=255)
                                               aaacaGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGa                                                                                                   >  4:350372/1‑149 (MQ=255)
                                               aaacaGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGa                                                                                                   >  1:194575/1‑149 (MQ=255)
                                               aaacaGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGACGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGa                                                                                                   <  8:299980/149‑1 (MQ=255)
                                                 acaGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTAGATCCCCAACACCGCGGTATCCACAACGGCGAt                                                                                                  >  1:172750/1‑148 (MQ=255)
                                                   aGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATaaa                                                                                               <  2:194575/149‑1 (MQ=255)
                                                    gTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAg                                                                                              <  3:278281/149‑1 (MQ=255)
                                                       tCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCa                                                                                           <  1:300767/149‑1 (MQ=255)
                                                         cGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACCCGTTATCCACAAGGGGAATaaatgaag                                                                                          >  8:34062/1‑143 (MQ=255)
                                                              cAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGGAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCtt                                                                                    >  8:115725/1‑149 (MQ=255)
                                                              cAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGGAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCtt                                                                                    >  8:121175/1‑149 (MQ=255)
                                                                 cTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTaa                                                                                 <  2:211823/149‑1 (MQ=255)
                                                                       ggTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCtt                                                                                                                                   >  2:76786/1‑93 (MQ=255)
                                                                       ggTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCtt                                                                                                                                   <  1:76786/93‑1 (MQ=255)
                                                                         tttCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCgtg                                                                         <  2:295208/149‑1 (MQ=255)
                                                                         tttCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCgtg                                                                         <  4:300026/149‑1 (MQ=255)
                                                                                tataAATCTCGCGTTCCCTCAACTTACTGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGACATGTGTATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAg                                                                                          <  7:267481/125‑1 (MQ=255)
                                                                                tataAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATaaa                                                                                               >  5:294391/1‑120 (MQ=255)
                                                                                tataAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATaaa                                                                                               <  6:294391/120‑1 (MQ=255)
                                                                                tataAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGt                                                                  >  1:146647/1‑149 (MQ=255)
                                                                                  taAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGtt                                                                 <  3:68979/148‑1 (MQ=255)
                                                                                  taAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGt                                                                  <  2:248564/147‑1 (MQ=255)
                                                                                    aaTCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGatcacctt                                                              >  2:57907/1‑142 (MQ=255)
                                                                                      tctcGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATaaa                                                                                               <  4:234009/114‑1 (MQ=255)
                                                                                       ctcGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATAtt                                                           >  8:231935/1‑149 (MQ=255)
                                                                                                  tCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATAtt                                                           >  7:141083/1‑138 (MQ=255)
                                                                                                  tCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATAtt                                                           <  8:141083/138‑1 (MQ=255)
                                                                                                    aaCTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGGAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGa                                              <  7:121175/149‑1 (MQ=255)
                                                                                                    aaCTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGGAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGa                                              <  7:115725/149‑1 (MQ=255)
                                                                                                    aaCTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGa                                                    >  4:168601/1‑143 (MQ=255)
                                                                                                    aaCTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGa                                                    <  3:168601/143‑1 (MQ=255)
                                                                                                    aaCTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGt                                                <  6:130446/147‑1 (MQ=255)
                                                                                                           ggCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCAc                                        >  5:154744/1‑148 (MQ=255)
                                                                                                                 gTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAAtgt                                                                                                                                                 <  2:192276/37‑1 (MQ=255)
                                                                                                                 gTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAAtgt                                                                                                                                                 >  1:192276/1‑37 (MQ=255)
                                                                                                                 gTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTaa                                                                                 <  4:138834/101‑1 (MQ=255)
                                                                                                                 gTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAAc                                                                                >  8:188253/1‑102 (MQ=255)
                                                                                                                 gTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAAc                                                                                <  7:188253/102‑1 (MQ=255)
                                                                                                                 gTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCgtggtg                                                                      >  3:239315/1‑112 (MQ=255)
                                                                                                                 gTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCgtggtg                                                                      <  4:239315/112‑1 (MQ=255)
                                                                                                                  ttGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAatgatg                                >  6:257430/1‑149 (MQ=255)
                                                                                                                      tgtgCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCgg                            <  2:146647/149‑1 (MQ=255)
                                                                                                                        tgCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATTCCCAAAAACGCGGTATCCACAACGTCGATAAAGTCAGGATCTTTAACGATCGTGGTTAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGgt                          <  1:35283/149‑1 (MQ=255)
                                                                                                                          cTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCGACGAATGATGCCGGgtgt                        >  8:343570/1‑149 (MQ=255)
                                                                                                                           tGAACGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGGGAGGTTCATATTGAGGCGAAAGTGACGCCCCGAAGGATGCCGGGtgtg                       >  7:366245/1‑149 (MQ=255)
                                                                                                                           tGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCGTATTGAGGCGAAAGTGACGCCACGAatga                                  <  2:71667/138‑1 (MQ=255)
                                                                                                                           tGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCGTATTGAGGCGAAAGTGACGCCACGAatga                                  >  1:71667/1‑138 (MQ=255)
                                                                                                                           tGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGtgtg                       >  3:321085/1‑149 (MQ=255)
                                                                                                                           tGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGtgtg                       >  1:151651/1‑149 (MQ=255)
                                                                                                                            gAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAAcgcg                                                                                                                   <  6:214540/56‑1 (MQ=255)
                                                                                                                            gAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAAcgcg                                                                                                                   >  5:214540/1‑56 (MQ=255)
                                                                                                                             aaaGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGt                     <  8:22522/149‑1 (MQ=255)
                                                                                                                                     ttGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTgg             >  1:314678/1‑149 (MQ=255)
                                                                                                                                       gCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTgg             <  6:154744/147‑1 (MQ=255)
                                                                                                                                           tCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTgg             <  4:198255/143‑1 (MQ=255)
                                                                                                                                           tCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTgg             >  3:198255/1‑143 (MQ=255)
                                                                                                                                           tCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAg       <  3:350372/149‑1 (MQ=255)
                                                                                                                                                aaaTGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGcc  <  7:174399/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TGCCGATCTACACGCCGGGCTTTGAAAGTTATCAGGATCCGCTGAACAAACAGTATCCGCTGCAGCTTACAGGTTTCCACTATAAATCTCGCGTTCACTCAACTTACGGCAACGTTGATGTGCTGAAAGCGGCTTGCCGTCAGGAAATGTGGATCAACCCGCTTGATGCCCAAAAACGCGGTATCCACAACGGCGATAAAGTCAGGATCTTTAACGATCGTGGTGAGGTTCATATTGAGGCGAAAGTGACGCCACGAATGATGCCGGGTGTGGTCGCACTGGGTGAAGGTGCC  >  NC_000913/942982‑943274

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: