Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,371,932 A→C 5.4% 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.4% 371.2 / 7.5 124E303A (GAA→GCA) groLCpn60 chaperonin GroEL, large subunit of GroESL
Reads supporting (aligned to +/- strand):  ref base A (46/71);  new base C (3/4);  total (49/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 = 9.80e-01

Only 100 of 129 total aligned reads displayed.

TGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGA  >  NC_000913/4371783‑4372081
                                                                                                                                                     |                                                                                                                                                     
tGTAGACGGCGAAGCGCCGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCGTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGa                                                                                                                                                       <  2:1481348/150‑1 (MQ=255)
tGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGa                                                                                                                                                       <  1:558451/150‑1 (MQ=255)
tGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGa                                                                                                                                                       >  1:682557/1‑150 (MQ=255)
tGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGa                                                                                                                                                       >  2:1051143/1‑150 (MQ=255)
   aGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGCCTGGCGGTACCGTGATCTCTGAaga                                                                                                                                                    >  1:1489016/1‑150 (MQ=255)
   aGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGGGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAaga                                                                                                                                                    <  2:1223980/150‑1 (MQ=255)
      aGGCGAAGCGCTGGCAACTCTGGTTGTTAACGCCATGCGTGGCATCGTGAAAGTCGCTGCGGTAAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGTATATCGCCACCCGGACTGGCCGTACCGTTATCTCTGAAgagag                                                                                                                                                 >  2:283075/1‑149 (MQ=255)
       ggCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATc                                                                                                                                                >  1:1098432/1‑150 (MQ=255)
       ggCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATc                                                                                                                                                <  1:219176/150‑1 (MQ=255)
       ggCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGGTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATc                                                                                                                                                <  1:1793122/150‑1 (MQ=255)
          gAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTTCAGAGATcacg                                                                                                                                             <  2:890052/150‑4 (MQ=255)
          gAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGt                                                                                                                                             >  1:686511/1‑150 (MQ=255)
           aaGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTTCAGAGATcacgg                                                                                                                                            <  2:456529/150‑5 (MQ=255)
             gcgcTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATg                                                                                                                                          <  2:1886724/150‑1 (MQ=255)
               gcTGGCAACTCTGGTGGTTAACACCATGCGTGGCATCGGGAAAGTCGCGGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGa                                                                                                                                        <  1:1282110/150‑1 (MQ=255)
                cTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATgagc                                                                                                                                       <  1:1734224/150‑4 (MQ=255)
                  ggCAACTCTGGTTGGTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCt                                                                                                                                     <  2:344306/150‑1 (MQ=255)
                    cAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTgg                                                                                                                                   <  1:910776/150‑1 (MQ=255)
                     aaCTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGa                                                                                                                                  <  2:1411989/150‑1 (MQ=255)
                      aCTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTAGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGaa                                                                                                                                 <  1:1144092/150‑1 (MQ=255)
                       ctctGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGaaa                                                                                                                                <  2:1621567/150‑1 (MQ=255)
                       ctctGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGAGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGaaa                                                                                                                                >  1:1457310/1‑150 (MQ=255)
                        tctGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGaaaa                                                                                                                               >  1:1152974/1‑150 (MQ=255)
                         ctGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGaaaaa                                                                                                                              >  1:657874/1‑150 (MQ=255)
                              tgttAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAAc                                                                                                                         <  1:464690/150‑1 (MQ=255)
                                  aaCACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTg                                                                                                                     >  2:889809/1‑150 (MQ=255)
                                  aaCACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTg                                                                                                                     <  1:1461073/150‑1 (MQ=255)
                                   acacCATGCGTGGCCTCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTgg                                                                                                                    <  2:1724418/150‑1 (MQ=255)
                                    cacCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGa                                                                                                                   >  1:1079430/1‑150 (MQ=255)
                                    cacCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGa                                                                                                                   >  2:688312/1‑150 (MQ=255)
                                        aTGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGAc                                                                                                               >  1:946287/1‑150 (MQ=255)
                                           cGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTg                                                                                                            <  2:1656781/150‑1 (MQ=255)
                                             tGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTggg                                                                                                          <  1:1346614/150‑1 (MQ=255)
                                             gggCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTggg                                                                                                          <  1:1885728/149‑1 (MQ=255)
                                               gCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTc                                                                                                        >  1:1218994/1‑150 (MQ=255)
                                               gCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTc                                                                                                        >  1:1057274/1‑150 (MQ=255)
                                                cATCGTGAAAGTCGCTGCGGTTTAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATTGAGCTGGAAAAAGCAACCCTGGAAGACCAGGGTcc                                                                                                       >  2:1939202/1‑149 (MQ=255)
                                                cATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCa                                                                                                       <  2:61095/150‑1 (MQ=255)
                                                  tCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAgg                                                                                                     >  1:1607023/1‑150 (MQ=255)
                                                   cGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGc                                                                                                    <  2:1133258/150‑1 (MQ=255)
                                                     tGAAATTCGCTGCGGTTAAAGCACCGGGCGTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCACCCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTa                                                                                                  <  2:203986/150‑1 (MQ=255)
                                                      gAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCACCCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAGCAACCCTGGAAGACCTGGGTCAGGCTaaa                                                                                                >  1:76901/1‑150 (MQ=255)
                                                         aGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACg                                                                                              <  2:520864/150‑1 (MQ=255)
                                                           tCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGAGTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAGGCCTGGAAGACCTGGGTCAGGCTAAACgtg                                                                                            >  2:277418/1‑150 (MQ=255)
                                                             gCTTCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGtgtt                                                                                          <  1:218872/150‑1 (MQ=255)
                                                              cTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGtgttg                                                                                         <  1:72511/150‑1 (MQ=255)
                                                                gCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTtgtg                                                                                       <  1:1174038/150‑1 (MQ=255)
                                                                     tAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATcaa                                                                                  >  1:740881/1‑150 (MQ=255)
                                                                       aaGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATcaaca                                                                                >  2:1702273/1‑150 (MQ=255)
                                                                        cgCACCGGGCTTCGGCGCTCGGCGTAAAGCTATGCTGCAGCATATCGCAACCCTGACTGGCGGTACCGTGATCTCAGAAGTGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGTCCTGGGTCAGGCTAAACGTGTTGTGATcaacaa                                                                               <  2:1303415/149‑1 (MQ=255)
                                                                          cACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGTATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTCGAAGACCTGGGTCAGGCTACACGTGTTGTGATCAACAAAg                                                                             <  2:1374257/150‑1 (MQ=255)
                                                                             cGGGGTTCGGCGAGCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGaca                                                                          <  1:606662/150‑1 (MQ=255)
                                                                              gggCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGacac                                                                         >  1:1138194/1‑150 (MQ=255)
                                                                               ggCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGGGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGAcacc                                                                        <  2:1152973/150‑1 (MQ=255)
                                                                                  ttCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGAcaccacc                                                                     <  2:1590143/150‑1 (MQ=255)
                                                                                     ggCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACATGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACt                                                                  >  2:1560102/1‑150 (MQ=255)
                                                                                      gCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTa                                                                 <  2:626975/150‑1 (MQ=255)
                                                                                       cGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTat                                                                >  1:1432151/1‑150 (MQ=255)
                                                                                         aTCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTatca                                                              <  2:843659/150‑1 (MQ=255)
                                                                                               gtAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATg                                                        <  1:889810/150‑1 (MQ=255)
                                                                                               gtAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTAATCAACAAAGACACCACCACTATCATCGATg                                                        >  2:63838/1‑150 (MQ=255)
                                                                                                    gCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTg                                                   <  1:2020438/150‑1 (MQ=255)
                                                                                                    gCTATGCTGAAGGTTATCGCTACGCTGATTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCTTTCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTg                                                   <  2:1071359/150‑1 (MQ=255)
                                                                                                      tATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTggg                                                 <  2:682557/150‑1 (MQ=255)
                                                                                                      tATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTggg                                                 <  1:2041893/150‑1 (MQ=255)
                                                                                                       aTGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGt                                                <  1:921953/150‑1 (MQ=255)
                                                                                                       aTGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGt                                                <  1:970203/150‑1 (MQ=255)
                                                                                                        tgctgcAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTg                                               <  1:618294/150‑1 (MQ=255)
                                                                                                         gctgcAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTga                                              <  1:257101/150‑1 (MQ=255)
                                                                                                           tgcAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGATAAAGCAACCCTGGTAGACCTGGGTCAGGCTAAACGTGTTGTGATCTACAAAGACACCACCACTATCATCGATGGCGTGGGTgaag                                            >  1:1235214/1‑150 (MQ=255)
                                                                                                               ggATATCGCCACCCTGACAGGCGGTACCGTGATCACTGAAGAGAGCGGTAAGGAGCTGGAAAAAGCAACCCTGGCAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGCCACCACCACTATCATCGATGGCGTGGGTGAAGAAGc                                        <  2:1342225/150‑1 (MQ=255)
                                                                                                                 atatCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTg                                      >  2:1102960/1‑150 (MQ=255)
                                                                                                                 atatCGCAACCCTGACTGGCGGAACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTg                                      <  2:278714/150‑1 (MQ=255)
                                                                                                                  tatCGCAACCCTGCCTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGc                                     <  1:293432/150‑1 (MQ=255)
                                                                                                                  tatCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGc                                     <  1:1548073/150‑1 (MQ=255)
                                                                                                                  tatCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGc                                     >  1:780278/1‑150 (MQ=255)
                                                                                                                    tCGCAACCCTTACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTTGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCaa                                   <  2:148690/150‑1 (MQ=255)
                                                                                                                     cGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAAt                                  >  2:230980/1‑150 (MQ=255)
                                                                                                                      gCAACCCTGACTGGCGGTACCGTGATCTCAAAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTCTCATCGATGGCGTGGGTGAACAAGCTGCAATc                                 >  1:819313/1‑150 (MQ=255)
                                                                                                                      ccAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGCAAAAAGCAACCCTGGAAGACCTGGTTCAGGCTAAACGTGTTGTGATCAACAAAAACACCACCACTATCAGCGATGGCGTGGCTGAAGAAGAAGCAATc                                 >  2:394043/2‑150 (MQ=255)
                                                                                                                           ccTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGAGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAggg                            >  2:434281/1‑150 (MQ=255)
                                                                                                                            cTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGc                           <  2:1489015/150‑1 (MQ=255)
                                                                                                                            cTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGc                           <  2:1198200/150‑1 (MQ=255)
                                                                                                                            cTGACTGGCGGTAACGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCAAGGGc                           >  1:1294085/1‑150 (MQ=255)
                                                                                                                             tGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGcc                          <  1:272149/150‑1 (MQ=255)
                                                                                                                                cTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCgtg                       >  2:221040/1‑150 (MQ=255)
                                                                                                                                 tGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCgtgt                      <  2:686511/150‑1 (MQ=255)
                                                                                                                                    cGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAGGCTGCAATCCAGGGCCGTGTTGc                   <  1:1883430/150‑1 (MQ=255)
                                                                                                                                       tACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCa                <  1:1702274/150‑1 (MQ=255)
                                                                                                                                       tACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACACCTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCa                >  1:1805160/1‑150 (MQ=255)
                                                                                                                                         ccGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACAACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGa              >  2:637545/1‑150 (MQ=255)
                                                                                                                                         ccGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAAACCTGGAAGACCTGGGTCAGGCGAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGGGATGATGCTGCTATCACGGGCCGTGTGGCTCAgc              >  2:361154/1‑149 (MQ=255)
                                                                                                                                          cGTGATCTCTGAAGAGATCGGTATGGTTCTGGAAAAAGCAACCCTGGAAGACCTGTGTCAGGCTAAACTTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCTTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGAt             <  1:1176429/150‑1 (MQ=255)
                                                                                                                                          cGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGAt             >  2:383790/1‑150 (MQ=255)
                                                                                                                                          accgATCTCTTCAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGAt             >  1:673591/4‑150 (MQ=255)
                                                                                                                                           gTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATc            >  1:1524857/1‑150 (MQ=255)
                                                                                                                                               tctcTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGGGATCAACAAAGACACCACCACTATCATCGACGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTc        <  1:1107948/150‑1 (MQ=255)
                                                                                                                                                ctctGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTca       >  1:1882662/1‑150 (MQ=255)
                                                                                                                                                 tctGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTcag      <  2:1138193/150‑1 (MQ=255)
                                                                                                                                                     aaGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGa  <  1:443532/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATCGGTATGGAGCTGGAAAAAGCAACCCTGGAAGACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGA  >  NC_000913/4371783‑4372081

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: