Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,183,258 G→A 6.5% E672K (GAA→AAA)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,183,2580GA 6.5% 479.4 / 12.6 155E672K (GAA→AAA) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (74/69);  new base A (8/2);  total (82/73)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 159 total aligned reads displayed.

ACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGG  >  NC_000913/4183109‑4183406
                                                                                                                                                     |                                                                                                                                                    
aCTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGCCTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTgg                                                                                                                                                      >  2:1838699/1‑150 (MQ=255)
aCTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTgg                                                                                                                                                      >  1:824778/1‑150 (MQ=255)
   tGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAac                                                                                                                                                   <  2:1849764/150‑1 (MQ=255)
    ggATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGAACAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAaca                                                                                                                                                  >  2:32677/1‑150 (MQ=255)
     gATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAacac                                                                                                                                                 >  1:217131/1‑150 (MQ=255)
         aagaagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTTTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTTAAACACGATg                                                                                                                                             >  2:365256/1‑150 (MQ=255)
         aagaagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATg                                                                                                                                             <  2:509071/150‑1 (MQ=255)
         aagaagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATg                                                                                                                                             >  2:1057140/1‑150 (MQ=255)
             agGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCTAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGc                                                                                                                                         <  2:125804/150‑1 (MQ=255)
             agGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGc                                                                                                                                         <  1:1194561/150‑1 (MQ=255)
             agGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGc                                                                                                                                         <  1:1740771/150‑1 (MQ=255)
              ggCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGcc                                                                                                                                        <  2:2030840/150‑1 (MQ=255)
                 cACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGTATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAAc                                                                                                                                     <  2:1501475/150‑1 (MQ=255)
                     tCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTCCATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCggg                                                                                                                                 >  2:386208/1‑148 (MQ=255)
                      cGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGc                                                                                                                                <  2:1284358/150‑1 (MQ=255)
                                 tGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGACCACGATGACGCCAACCGTGCATTGATGGGTg                                                                                                                     >  1:1283626/1‑150 (MQ=255)
                                  ggTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGCTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAGCACGATGACGCCACCCCTGCATTGATGGGAGc                                                                                                                    >  1:816558/1‑150 (MQ=255)
                                   gTAACTTGCCGTAGCATAGGCGAATCCAGGTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCg                                                                                                                   <  2:1572193/150‑1 (MQ=255)
                                   gTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCg                                                                                                                   <  2:57092/150‑1 (MQ=255)
                                    tAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCCCCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCACGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGa                                                                                                                  >  2:1126941/1‑150 (MQ=255)
                                         tGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATg                                                                                                             <  1:1579660/150‑1 (MQ=255)
                                         tGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATg                                                                                                             <  2:2033915/150‑1 (MQ=255)
                                         tGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATg                                                                                                             >  1:230554/1‑150 (MQ=255)
                                         tGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCACGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGGTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATg                                                                                                             >  2:1391142/1‑150 (MQ=255)
                                          gCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGCACAAGc                                                                                                            >  1:1117985/1‑150 (MQ=255)
                                            cGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCaa                                                                                                          <  1:103529/150‑1 (MQ=255)
                                            cGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCaa                                                                                                          >  1:1754249/1‑150 (MQ=255)
                                               aGCAAAGGCGAATCCAGCTTGTTCAGCCGCGCCCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGt                                                                                                       >  2:117766/1‑150 (MQ=255)
                                               aGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGCCGCCAACCGTGCATTGATGGGTGCGAACATGCCACGt                                                                                                       >  1:1799929/1‑150 (MQ=255)
                                               aGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGCCGCCAACCGTGCATTGATGGGTGCGAACATGCAACGt                                                                                                       >  2:42270/1‑150 (MQ=255)
                                                 cAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCa                                                                                                     <  1:35238/150‑1 (MQ=255)
                                                 cAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCCTCGGTGCGTCCCTGATCCCGTTCCTGATACACGATGACGCCAACCGTGCATTGATGGGTGCGCACATGCAACGTCa                                                                                                     >  1:361937/1‑150 (MQ=255)
                                                  aaaGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAg                                                                                                    >  1:1837031/1‑150 (MQ=255)
                                                    aGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGc                                                                                                  <  1:1832495/150‑1 (MQ=255)
                                                    aGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGc                                                                                                  >  1:31286/1‑150 (MQ=255)
                                                    aGGCGAATCCAGCTTGTGCAGCCGCGTCCAGGCTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTGCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGc                                                                                                  <  1:18575/150‑1 (MQ=255)
                                                     ggCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGcc                                                                                                 >  1:851299/1‑150 (MQ=255)
                                                       cGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCgg                                                                                               >  2:1900549/1‑149 (MQ=255)
                                                       cGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGt                                                                                               >  2:792181/1‑150 (MQ=255)
                                                          aTCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTcc                                                                                            >  1:990353/1‑150 (MQ=255)
                                                            ccAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGa                                                                                          <  1:1804280/150‑1 (MQ=255)
                                                             cAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGAc                                                                                         >  2:748112/1‑150 (MQ=255)
                                                             cAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGAc                                                                                         >  1:1973997/1‑150 (MQ=255)
                                                               gCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGActc                                                                                       <  1:754900/150‑1 (MQ=255)
                                                                 ttGTTCAGCCGCGACCAGGTGGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCCGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTg                                                                                     <  2:1731946/150‑1 (MQ=255)
                                                                   gTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcg                                                                                   <  2:551622/150‑1 (MQ=255)
                                                                      cAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcgcgc                                                                                <  1:163724/150‑1 (MQ=255)
                                                                      cAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGGCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATCGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcgcgc                                                                                >  2:1275320/1‑150 (MQ=255)
                                                                       aGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGACCACGATGCCGCCAACCGAGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCt                                                                               >  1:1275919/1‑150 (MQ=255)
                                                                       aGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCt                                                                               <  1:1904486/150‑1 (MQ=255)
                                                                       aGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCt                                                                               >  1:1405158/1‑150 (MQ=255)
                                                                        gacgcgACCAGGTTGACTACATGAACTTATCCACCCAGCACGTGGTATCCCTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTCCATGGCTCGGAGAGAACATGCAACCTCCGGCCGATCCGCCTCTGCGCGCTg                                                                              >  2:823087/3‑150 (MQ=255)
                                                                           gcgTCCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCTACGTCAGGCCGTTCCGACTCTGCGCGCTGATa                                                                           <  2:373809/150‑1 (MQ=255)
                                                                               ccAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTCCGCGCTGATAAGcc                                                                       >  2:899015/1‑150 (MQ=255)
                                                                                 aGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACTGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGc                                                                     <  2:1564756/150‑1 (MQ=255)
                                                                                 aGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGc                                                                     >  1:883809/1‑150 (MQ=255)
                                                                                  ggTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCt                                                                    >  2:750289/1‑150 (MQ=255)
                                                                                  ggTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGCCGGTGCGTCCCTGATCCCGTTCCTGTAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCATCGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCt                                                                    <  2:1085069/150‑1 (MQ=255)
                                                                                   gTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCACTGATAAGCCGCtg                                                                   <  2:955539/150‑1 (MQ=255)
                                                                                     tGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggt                                                                 >  1:607792/1‑150 (MQ=255)
                                                                                     tGACTACATGGACGTATCCACCCAGCAGGGGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggt                                                                 <  2:1245741/150‑1 (MQ=255)
                                                                                      gACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGACGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggat                                                                >  2:1782127/1‑148 (MQ=255)
                                                                                           cATGCACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTAc                                                           >  1:959326/1‑150 (MQ=255)
                                                                                             tGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGGGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTg                                                         <  2:68598/150‑1 (MQ=255)
                                                                                              ggACGTATCCACCCAGCCGGTGGTATCCGTCGGTGCGACCCTGATCCCGTTCCTGGAACACGATCACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTCCGCGCGGACAAGCCGCTGGTTGGTACTgg                                                        <  2:385276/150‑1 (MQ=255)
                                                                                              ggACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTgg                                                        <  2:1370226/150‑1 (MQ=255)
                                                                                               gACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGt                                                       <  2:1101191/150‑1 (MQ=255)
                                                                                                 cGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTAGGTACTGGTAt                                                     >  2:842044/1‑150 (MQ=255)
                                                                                                  gtctCCACCCAGCTGGTGGTATCCGTCGGGGCGTCCCGGATCGCGGTGCTGGAACACGCTCCCGCCAACCGTGCCGTGAGGGGGGAGAACATGCCACGTCAGGCCGTTCCGTCCCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTCTg                                                    <  2:1013463/147‑1 (MQ=255)
                                                                                                  gTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATg                                                    <  1:32677/150‑1 (MQ=255)
                                                                                                   tATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCACCGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATgg                                                   >  1:129657/1‑150 (MQ=255)
                                                                                                   tATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATgg                                                   >  1:1356464/1‑150 (MQ=255)
                                                                                                       cACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACg                                               >  2:1774187/1‑150 (MQ=255)
                                                                                                       cACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACg                                               >  1:747700/1‑150 (MQ=255)
                                                                                                         cccAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCGGGAACACGATGACGCCAACCGTGCATTGAGGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTAGGGAACGTg                                             >  2:433141/1‑150 (MQ=255)
                                                                                                          ccAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGc                                            <  2:1767663/150‑1 (MQ=255)
                                                                                                          ccAGCAGGTGGTATCCGTCGGAGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGc                                            <  1:1101096/150‑1 (MQ=255)
                                                                                                                ggtggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGc                                      <  2:642921/150‑1 (MQ=255)
                                                                                                                     tATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGGTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTg                                 <  2:675158/150‑1 (MQ=255)
                                                                                                                     tATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTg                                 <  1:929571/150‑1 (MQ=255)
                                                                                                                      aTCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGa                                >  1:915251/1‑150 (MQ=255)
                                                                                                                       tCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGtagcc                               >  1:1486131/1‑146 (MQ=255)
                                                                                                                        ccGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGCCGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACt                              >  2:520245/1‑150 (MQ=255)
                                                                                                                         cGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTc                             >  2:879567/1‑150 (MQ=255)
                                                                                                                            cGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCgg                          <  1:792182/150‑1 (MQ=255)
                                                                                                                            cGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCgg                          <  1:619227/150‑1 (MQ=255)
                                                                                                                                gCGGCCCTGATCCCTTTCCTGCAACACGCTGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGCTCCGTCTCTGCGCGCTGATAAGCCGCGGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTa                      <  1:365256/150‑1 (MQ=255)
                                                                                                                                  gTCCCTGATCCCGTTCNTGGAACACGATGACGCCAACCGTGCATAGATGGGTGCGCACATGCAACGTCAGGCCGGTCCGACTCTGCGCGCTGATAAGCCGCTGTTTGGTACTGGTATGGAACGTGCTGTTGCCCTGGACTCCGGTGTAAc                    >  1:849667/1‑150 (MQ=255)
                                                                                                                                  gTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAAc                    >  1:921108/1‑150 (MQ=255)
                                                                                                                                      cTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGCTGGGTGCGAACATGCACCGTCTGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCGGTTGCCGTTGACTCCGGTGTAACTGCg                <  2:472653/150‑1 (MQ=255)
                                                                                                                                        gATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCGGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGt              <  1:1450807/150‑1 (MQ=255)
                                                                                                                                         aTCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTa             >  1:465779/1‑150 (MQ=255)
                                                                                                                                          tCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGTTGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAg            <  2:1615860/150‑1 (MQ=255)
                                                                                                                                          tCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAg            >  1:1577207/1‑150 (MQ=255)
                                                                                                                                          tCCCGTTCCTGGAACACGATGACGCCAACCGTGCAATGATGCGTGCGAACATGCAACGTCAGGCCGTTCCGACACTCAGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAg            >  1:1095351/1‑150 (MQ=255)
                                                                                                                                          ccccGTTCCTGGAACACGATGACGCCAACCGTGCAGTGAGGGGTGCGAACATGCAACGGCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAg            <  1:1360625/149‑1 (MQ=255)
                                                                                                                                           cccGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGc           >  1:100187/1‑150 (MQ=255)
                                                                                                                                           cccGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGGTGCCGTTGACTCCGGTGTAACTGCGGTAGc           <  2:224084/150‑1 (MQ=255)
                                                                                                                                            ccGGTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCt          <  2:1986001/150‑1 (MQ=255)
                                                                                                                                                    ggAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGGCAGGCCGTTCCGACGCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACgtgg  <  2:1283625/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGG  >  NC_000913/4183109‑4183406

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: