Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,371,932 A→C 5.8% E303A (GAA→GCA) ‡ groL → Cpn60 chaperonin GroEL, large subunit of GroESL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,371,9320AC 5.8% 471.4 / 18.1 140E303A (GAA→GCA) ‡groLCpn60 chaperonin GroEL, large subunit of GroESL
Reads supporting (aligned to +/- strand):  ref base A (60/71);  new base C (4/4);  total (65/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.87e-01

Only 100 of 144 total aligned reads displayed.

TAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTCAGCA  >  NC_000913/4371785‑4372079
                                                                                                                                                   |                                                                                                                                                   
tAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGaa                                                                                                                                                    <  3:350396/149‑1 (MQ=255)
tAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGaa                                                                                                                                                    >  3:680540/1‑149 (MQ=255)
tAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGaa                                                                                                                                                    <  4:126380/149‑1 (MQ=255)
tAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGaa                                                                                                                                                    <  2:199086/149‑1 (MQ=255)
 aGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAag                                                                                                                                                   >  1:545450/1‑149 (MQ=255)
     ggCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGAt                                                                                                                                               <  4:743168/149‑1 (MQ=255)
       cGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCg                                                                                                                                             <  3:24333/149‑1 (MQ=255)
        gAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCgg                                                                                                                                            >  3:131596/1‑149 (MQ=255)
        gAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCgg                                                                                                                                            >  3:834223/1‑149 (MQ=255)
        gAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCgg                                                                                                                                            >  2:169798/1‑149 (MQ=255)
         aaGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTTCAGAGATcacg                                                                                                                                           <  2:90159/149‑4 (MQ=255)
         aaGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTTCAGAGATcacg                                                                                                                                           <  3:25568/149‑4 (MQ=255)
         aaGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGt                                                                                                                                           >  1:107353/1‑149 (MQ=255)
          aGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTa                                                                                                                                          <  1:51137/149‑1 (MQ=255)
           gcgcTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACC‑‑GATCTCTTCAGAGATcacggtac                                                                                                                                       <  1:498119/149‑8 (MQ=255)
               tGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAg                                                                                                                                     >  2:779304/1‑149 (MQ=255)
                 gCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCt                                                                                                                                   <  1:547853/149‑1 (MQ=255)
                 gCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCt                                                                                                                                   <  4:459694/149‑1 (MQ=255)
                   aaCTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTgg                                                                                                                                 >  3:781404/1‑149 (MQ=255)
                   aaCTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTgg                                                                                                                                 >  4:193003/1‑149 (MQ=255)
                      tctGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGaaa                                                                                                                              >  3:903605/1‑149 (MQ=255)
                        tGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGaaaaa                                                                                                                            <  2:713478/149‑1 (MQ=255)
                        tGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGaaaaa                                                                                                                            <  4:417360/149‑1 (MQ=255)
                         ggTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGTGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAg                                                                                                                           <  1:779436/149‑1 (MQ=255)
                          gttgttAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGc                                                                                                                          >  1:394148/1‑149 (MQ=255)
                             gttAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGTGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAAc                                                                                                                       >  4:343942/1‑149 (MQ=255)
                             gttAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAAc                                                                                                                       >  4:649009/1‑149 (MQ=255)
                                aaCACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCt                                                                                                                    >  4:177336/1‑149 (MQ=255)
                                aaCACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCt                                                                                                                    <  2:394148/149‑1 (MQ=255)
                                  cacCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTgg                                                                                                                  >  3:319762/1‑149 (MQ=255)
                                     cATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAaga                                                                                                                                                  >  4:757246/1‑114 (MQ=255)
                                     cATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAaga                                                                                                                                                  <  3:757246/114‑1 (MQ=255)
                                     cATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGaaaaa                                                                                                                            >  2:517573/1‑136 (MQ=255)
                                       tGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTCATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGAc                                                                                                             >  4:889435/1‑149 (MQ=255)
                                           tGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTgg                                                                                                         <  4:570101/149‑1 (MQ=255)
                                             gCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTTGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGt                                                                                                       >  4:330699/1‑149 (MQ=255)
                                               aTCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCa                                                                                                     >  2:406242/1‑149 (MQ=255)
                                               aTCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCa                                                                                                     >  2:106714/1‑149 (MQ=255)
                                                 cGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAgg                                                                                                   <  1:226382/149‑1 (MQ=255)
                                                    gAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTa                                                                                                >  2:875493/1‑149 (MQ=255)
                                                     aaaGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTaa                                                                                               <  1:755554/149‑1 (MQ=255)
                                                       aGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAAc                                                                                             <  4:319762/149‑1 (MQ=255)
                                                         tCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACgt                                                                                           <  3:762408/149‑1 (MQ=255)
                                                           gCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACgtgt                                                                                         <  4:521363/149‑1 (MQ=255)
                                                            cTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTGCCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGtgtt                                                                                        <  3:649009/149‑1 (MQ=255)
                                                               cGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTtgtg                                                                                     <  1:108178/149‑1 (MQ=255)
                                                                 gTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGtgtt                                                                                        >  4:858805/1‑144 (MQ=255)
                                                                 gTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGtgtt                                                                                        <  3:858805/144‑1 (MQ=255)
                                                                 gTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGAt                                                                                   >  1:699802/1‑149 (MQ=255)
                                                                  ttAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATc                                                                                  <  2:733880/149‑1 (MQ=255)
                                                                    ataGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGATTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATcaa                                                                                <  4:794678/147‑1 (MQ=255)
                                                                    aaaGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATcaa                                                                                <  1:222587/149‑1 (MQ=255)
                                                                     aaGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATcaac                                                                               >  2:632315/1‑149 (MQ=255)
                                                                         aCCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAg                                                                           <  1:374568/149‑1 (MQ=255)
                                                                          ccGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGa                                                                          <  2:268857/149‑1 (MQ=255)
                                                                             ggCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGacac                                                                       <  1:632315/149‑1 (MQ=255)
                                                                             ggCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGacac                                                                       <  2:322743/149‑1 (MQ=255)
                                                                              gCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGAcacc                                                                      >  4:703146/1‑149 (MQ=255)
                                                                                ttCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGAcaccac                                                                    >  2:208675/1‑149 (MQ=255)
                                                                                       aTCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTatc                                                             >  1:383713/1‑149 (MQ=255)
                                                                                            cgtAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGa                                                        <  3:182436/149‑1 (MQ=255)
                                                                                                 aGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGAcaccaccac                                                                 <  2:354727/135‑1 (MQ=255)
                                                                                                 aGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGAcaccaccac                                                                 >  1:354727/1‑135 (MQ=255)
                                                                                                 aGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCg                                                   >  2:76648/1‑149 (MQ=255)
                                                                                                    tATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTgg                                                >  2:271323/1‑149 (MQ=255)
                                                                                                    tATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTgg                                                <  3:177336/149‑1 (MQ=255)
                                                                                                    tATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTgg                                                <  1:875493/149‑1 (MQ=255)
                                                                                                     aTGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTggg                                               <  4:413496/149‑1 (MQ=255)
                                                                                                     aTGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTggg                                               <  3:178908/149‑1 (MQ=255)
                                                                                                        ctgcAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTga                                            <  3:434398/149‑1 (MQ=255)
                                                                                                         tgcAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTgaa                                           >  1:608466/1‑149 (MQ=255)
                                                                                                          gcAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGt                                                  <  4:257998/141‑1 (MQ=255)
                                                                                                          gcAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGt                                                  >  3:257998/1‑141 (MQ=255)
                                                                                                              gATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGc                                      <  3:129960/149‑1 (MQ=255)
                                                                                                               atatCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCt                                     <  1:664488/149‑1 (MQ=255)
                                                                                                                tatCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTg                                    <  1:749391/149‑1 (MQ=255)
                                                                                                                    ccgaCCCTGACTGTCTGTACCGTGATCTCTGAAGAGATCGGTATGGTGCTGGAAAAAGCAACCCTCGAAGACCTGGGGAAGGCTAAACGTGTTGTGATCAACGAGGTCACCACCACTATGATCGATGGCGGGGGTGAAGAAGCTGCAAt                                >  2:222619/4‑149 (MQ=255)
                                                                                                                         ccTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAgg                           >  3:311815/1‑149 (MQ=255)
                                                                                                                         ccTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAgg                           >  4:162124/1‑149 (MQ=255)
                                                                                                                          cTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAggg                          <  3:703146/149‑1 (MQ=255)
                                                                                                                          cTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAggg                          >  3:194638/1‑149 (MQ=255)
                                                                                                                           tGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGc                         <  1:292361/149‑1 (MQ=255)
                                                                                                                            gACTGGCGGTACCGTGATCTCTGAAGAGATCGTTATGGAGCTGGAAAAAGTAACCCTGGAAGACCTGGGTCAGGCTAAACGTATTTTGATTAACAAAGACACCACCACTATCGTCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGcc                        <  3:778017/149‑1 (MQ=255)
                                                                                                                             aCTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTgg                                                >  2:547955/1‑124 (MQ=255)
                                                                                                                               tGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCgtg                     >  4:766933/1‑149 (MQ=255)
                                                                                                                                 gCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGtt                   >  2:505571/1‑149 (MQ=255)
                                                                                                                                   ggTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAAc                                                                                             >  2:459622/1‑73 (MQ=255)
                                                                                                                                    gTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGc                                      <  4:145961/127‑1 (MQ=255)
                                                                                                                                     tACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTc               <  2:545450/149‑1 (MQ=255)
                                                                                                                                     tACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTc               <  3:448861/149‑1 (MQ=255)
                                                                                                                                       taccgATCTCTTCAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAg             >  2:632768/5‑149 (MQ=255)
                                                                                                                                         gTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGAt           >  2:755406/1‑149 (MQ=255)
                                                                                                                                          tGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATc          >  3:585417/1‑149 (MQ=255)
                                                                                                                                          cgATCTCTTCAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATc          >  4:384057/2‑149 (MQ=255)
                                                                                                                                          cgATCTCTTCAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATc          >  1:645373/2‑149 (MQ=255)
                                                                                                                                                ctGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTcag    >  4:17724/1‑149 (MQ=255)
                                                                                                                                                 tGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGc                 >  3:91553/1‑135 (MQ=255)
                                                                                                                                                  gAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACaaa                                                                            <  4:670829/75‑1 (MQ=255)
                                                                                                                                                  gAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACaaa                                                                            >  3:670829/1‑75 (MQ=255)
                                                                                                                                                  gAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTcagca  <  1:138422/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTCAGCA  >  NC_000913/4371785‑4372079

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: