Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,154,604 C→G 5.4% R197G (CGG→GGG)  mdtA → multidrug efflux system, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,154,6040CG 5.4% 321.3 / 5.1 102R197G (CGG→GGG) mdtAmultidrug efflux system, subunit A
Reads supporting (aligned to +/- strand):  ref base C (44/52);  new base G (3/3);  total (47/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 102 total aligned reads displayed.

CCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTACCCTGCCGGAAAGCGATATC  >  NC_000913/2154456‑2154750
                                                                                                                                                    |                                                                                                                                                  
ccTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCg                                                                                                                                                   >  2:1447120/1‑150 (MQ=255)
ccTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCg                                                                                                                                                   >  2:392181/1‑150 (MQ=255)
ccTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCg                                                                                                                                                   <  1:1349009/150‑1 (MQ=255)
 cTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCgg                                                                                                                                                  <  2:557425/150‑1 (MQ=255)
     cgcgTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTa                                                                                                                                              >  2:510823/1‑150 (MQ=255)
           aTCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCAc                                                                                                                                        <  2:1711548/150‑1 (MQ=255)
            tCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTTGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCAcc                                                                                                                                       <  1:1483739/150‑1 (MQ=255)
            tCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCAcc                                                                                                                                       >  1:382039/1‑150 (MQ=255)
             caacaaCTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCAccc                                                                                                                                      >  1:1694002/1‑149 (MQ=255)
              aacaacTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGTATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCAGATGAAGCAAGCGTTGCCAGCGAGCAGCTGCAAATCGACTGGAGCCGGATTACCGCACCAg                                                                                                                                     >  2:347661/1‑150 (MQ=255)
              aacaacTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCTGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCGGCATTACCTCACCAg                                                                                                                                     >  2:1672646/1‑150 (MQ=255)
                aaacTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTc                                                                                                                                   <  2:395741/149‑1 (MQ=255)
                  acTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGa                                                                                                                                 <  1:951203/150‑1 (MQ=255)
                  acTGGCAAAAACCAAACTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGa                                                                                                                                 >  1:586494/1‑150 (MQ=255)
                   cTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGAt                                                                                                                                >  1:1795783/1‑150 (MQ=255)
                     ggCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGTGCGGTAATCCGGCTCCAGTCGagtt                                                                                                                              <  2:635116/150‑4 (MQ=255)
                     ggCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATgg                                                                                                                              >  2:1130607/1‑150 (MQ=255)
                      gCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGt                                                                                                                             >  1:44316/1‑150 (MQ=255)
                       cAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGAGGGTc                                                                                                                            >  1:1522415/1‑150 (MQ=255)
                           aaCCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGt                                                                                                                        >  2:13900/1‑150 (MQ=255)
                             ccAATCTCGTTTCCCGCCAGGAGCTGGATGCCCCACCGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCCGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTg                                                                                                                      >  2:757384/1‑150 (MQ=255)
                             ccAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTg                                                                                                                      >  1:528658/1‑150 (MQ=255)
                                aTCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGtc                                                                                                                   <  2:418519/150‑1 (MQ=255)
                                aTCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATCGTCGCGTTGGtc                                                                                                                   >  2:791826/1‑150 (MQ=255)
                                 tctcGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGtct                                                                                                                  >  1:1207166/1‑150 (MQ=255)
                                   tcGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCa                                                                                                                >  1:565548/1‑150 (MQ=255)
                                   tcGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGGTGGTCTCa                                                                                                                >  2:1362821/1‑150 (MQ=255)
                                   ccGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCa                                                                                                                <  1:651794/149‑1 (MQ=255)
                                    cGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCaa                                                                                                               <  1:1816525/150‑1 (MQ=255)
                                    cGTTGCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCTAACCAGTCGATGGTCGCGTTGGTCTCaa                                                                                                               <  2:285382/150‑1 (MQ=255)
                                             ccAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGa                                                                                                      <  2:673009/150‑1 (MQ=255)
                                             ccAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGa                                                                                                      >  2:1608414/1‑150 (MQ=255)
                                               aGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATg                                                                                                    <  2:629272/150‑1 (MQ=255)
                                                 gAGCTGGATGCCCAACAGGCGCTGGGCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGtt                                                                                                  <  1:1715006/150‑1 (MQ=255)
                                                   gCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTgg                                                                                                <  1:1033186/150‑1 (MQ=255)
                                                   gCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTgg                                                                                                <  1:1821654/150‑1 (MQ=255)
                                                    cTGGATGCCCAACAGGCGCTGGTCAGTGAAACCTAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGt                                                                                               <  2:1453974/150‑1 (MQ=255)
                                                       gATGCCCAACAGGCGCTGGGCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAAc                                                                                            <  1:826177/150‑1 (MQ=255)
                                                         tGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCa                                                                                          >  1:800081/1‑150 (MQ=255)
                                                            ccAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAAt                                                                                       >  2:1774468/1‑150 (MQ=255)
                                                            ccAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAAt                                                                                       >  2:240839/1‑150 (MQ=255)
                                                                ccggCGCTGTCCAGTGAAACCGAAGGCCACATTCAGGCTGATGAAGCACGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGCGCCGGATTACCGCACCCGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGCTGTTGGTAACCAAATCTcc                                                                                   <  2:1806567/148‑1 (MQ=255)
                                                                cAGGCGCTGGTCACTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTcc                                                                                   <  2:138759/150‑1 (MQ=255)
                                                                  ggCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATTGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAg                                                                                 <  1:37726/150‑1 (MQ=255)
                                                                      cTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGGAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTTGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAgtggt                                                                             <  1:347661/150‑1 (MQ=255)
                                                                       ttgTCTGTGAAACCTATGGCACCATTAAGGCTGATGAAGCAAGCGTTTCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAgtggtg                                                                            <  1:1375671/148‑1 (MQ=255)
                                                                       tGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAgtggtg                                                                            <  2:737259/150‑1 (MQ=255)
                                                                         gTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGAt                                                                          <  1:1742574/150‑1 (MQ=255)
                                                                           cAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATac                                                                        >  2:675547/1‑150 (MQ=255)
                                                                             gTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATacca                                                                      >  1:635030/1‑150 (MQ=255)
                                                                             gTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGGCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATaccc                                                                      >  2:123286/1‑149 (MQ=255)
                                                                               gAAACCGAATGCACCATTAAGGCTGATGAAGCAAGAGTTGCCAGCGCGCAGCTGCAACTCGACTCGAGACTGATAACCGCACCAGTCGAAGGTCCCGTAGGTCTCAAGCAGGTTGATGTAGGTAACCAAATCTCCAGTGGTGATaccaac                                                                    >  2:65019/1‑148 (MQ=255)
                                                                               gAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATaccacc                                                                    >  2:784106/1‑150 (MQ=255)
                                                                               gAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATaccacc                                                                    >  2:88426/1‑150 (MQ=255)
                                                                                 aaCCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCgg                                                                  <  2:1133678/150‑1 (MQ=255)
                                                                                  cccGAAGGCACCATTACGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCggg                                                                 <  2:1429023/149‑1 (MQ=255)
                                                                                  aCCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCggg                                                                 <  1:1628721/150‑1 (MQ=255)
                                                                                     gAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATc                                                              <  2:803518/150‑1 (MQ=255)
                                                                                     gAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATc                                                              >  2:1655969/1‑150 (MQ=255)
                                                                                      aaGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCg                                                             <  1:1074408/150‑1 (MQ=255)
                                                                                      aaGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCg                                                             <  1:644000/150‑1 (MQ=255)
                                                                                       aGGCACCATTAAGGCGGATGAAGCAAGCGGTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGGTGGTCTCAAGCAGGTGGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCgt                                                            <  2:899252/150‑1 (MQ=255)
                                                                                         gCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCgtgg                                                          <  1:461116/150‑1 (MQ=255)
                                                                                         gCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCgtgg                                                          <  1:286452/150‑1 (MQ=255)
                                                                                         gCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCgtgg                                                          <  1:1658042/150‑1 (MQ=255)
                                                                                         gCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCgtgg                                                          <  2:252703/150‑1 (MQ=255)
                                                                                             cATTAAGGCTGATGACGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGAt                                                      <  2:291638/150‑1 (MQ=255)
                                                                                             cATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCCGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGAt                                                      <  2:1417366/150‑1 (MQ=255)
                                                                                             cATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGAt                                                      >  2:915504/1‑150 (MQ=255)
                                                                                              aTTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATc                                                     <  2:753221/150‑1 (MQ=255)
                                                                                              aTTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATc                                                     <  2:1072624/150‑1 (MQ=255)
                                                                                                  aGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGCTCAccc                                                 >  2:872905/1‑150 (MQ=255)
                                                                                                  aGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCAccc                                                 <  2:1795783/150‑1 (MQ=255)
                                                                                                   ggCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCa                                                >  1:1186625/1‑150 (MQ=255)
                                                                                                    gCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAg                                               >  1:521126/1‑150 (MQ=255)
                                                                                                      tgatgaAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGAc                                             >  1:1310521/1‑150 (MQ=255)
                                                                                                         tgaAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCa                                          >  1:1115217/1‑150 (MQ=255)
                                                                                                         tgaAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACAACCGGGATCGTGGTGATCAACCAGACGcc                                          >  2:911996/1‑149 (MQ=255)
                                                                                                              caaccGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCCGGTTGAGGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTa                                     <  1:627845/150‑1 (MQ=255)
                                                                                                               aagcGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTAt                                    >  2:459150/1‑150 (MQ=255)
                                                                                                                agcGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATc                                   <  2:1323993/150‑1 (MQ=255)
                                                                                                                agcGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGGCGCGTTGGTCTCAAGCAGGGTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATc                                   <  2:1064117/150‑1 (MQ=255)
                                                                                                                 gcGTTGCCAGCGCGCAGCTGCAACTCGACTGGATCCGGATTACCGCCCCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCg                                  >  1:93177/1‑150 (MQ=255)
                                                                                                                  cGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGa                                 <  1:13900/150‑1 (MQ=255)
                                                                                                                      gCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCTCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTCGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTa                             <  1:783424/150‑1 (MQ=255)
                                                                                                                       ccAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAg                            <  2:570853/150‑1 (MQ=255)
                                                                                                                       ccAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAg                            >  1:476645/1‑150 (MQ=255)
                                                                                                                         aGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTc                          >  2:1059668/1‑150 (MQ=255)
                                                                                                                            gcgcAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCttt                       <  2:1788297/150‑1 (MQ=255)
                                                                                                                             cgcAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTa                      >  2:1325648/1‑150 (MQ=255)
                                                                                                                                 gCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCCCCGGGATCGTGGTGATCACCCAGCCGCATCCTATCGATTTAGTCTTTACCCt                  >  1:256859/1‑150 (MQ=255)
                                                                                                                                  cTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTACCCTg                 <  1:822282/150‑1 (MQ=255)
                                                                                                                                      ccatCGACTGGTGCGGTAATCCGGCTCCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTACCCTGCCgg             >  1:913442/4‑150 (MQ=255)
                                                                                                                                      aaCTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGAATGTTGATGTTTGTAACCAAATCTCCAGTGGTTATACCACCGGGATCTTGGTGATCACCCATACGCATGGTATGTATTTAGTCTTTACCCTGCCgg             <  2:586494/150‑1 (MQ=255)
                                                                                                                                       catCGACTGGTGCGGTAATCCGGCTCCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGATGCATCCTATCGATTTAGTCTTTACCCTGCCGGa            >  2:499710/3‑150 (MQ=255)
                                                                                                                                             cTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTACCCTGCCGGAAAGCGa      <  2:1272580/150‑1 (MQ=255)
                                                                                                                                                gAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTACCCTGCCGGAAAGCGatat   >  2:235338/1‑150 (MQ=255)
                                                                                                                                                 aGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTACCCTGCCGGAAAGCGATATc  <  1:1612015/150‑1 (MQ=255)
                                                                                                                                                 aGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTACCCTGCCGGAAAGCGATATc  >  1:786620/1‑150 (MQ=255)
                                                                                                                                                 aGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTACCCTGCCGGAAAGCGATATc  <  1:510823/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGGTTGATGTTGGTAACCAAATCTCCAGTGGTGATACCACCGGGATCGTGGTGATCACCCAGACGCATCCTATCGATTTAGTCTTTACCCTGCCGGAAAGCGATATC  >  NC_000913/2154456‑2154750

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