Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,147,021 T→C 100% M174V (ATG→GTG)  alkA ← 3‑methyl‑adenine DNA glycosylase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,147,0210TC96.9% 2884.5 / 41.8 792M174V (ATG→GTG) alkA3‑methyl‑adenine DNA glycosylase II
Reads supporting (aligned to +/- strand):  ref base T (10/15);  new base C (358/409);  total (368/424)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 801 total aligned reads displayed.

AAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCACCAGTT  >  NC_000913/2146872‑2147170
                                                                                                                                                     |                                                                                                                                                     
   ataatTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACAcc                                                                                                                                                    >  1:7699588/1‑150 (MQ=255)
    taatTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCt                                                                                                                                                   <  2:2245187/150‑1 (MQ=255)
     aatTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTa                                                                                                                                                  <  2:9113884/150‑1 (MQ=255)
     aatTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACAGCGCCCGGTATGGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTa                                                                                                                                                  >  2:287332/1‑150 (MQ=255)
       ttCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAAc                                                                                                                                                >  1:6441788/1‑150 (MQ=255)
       ttCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAAc                                                                                                                                                >  2:7958480/1‑150 (MQ=255)
              tCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTa                                                                                                                                         >  2:2773825/1‑150 (MQ=255)
              tCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTa                                                                                                                                         >  2:11715765/1‑150 (MQ=255)
              tCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTa                                                                                                                                         <  1:1469672/150‑1 (MQ=255)
                     cccGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTg                                                                                                                                  <  2:538454/150‑1 (MQ=255)
                     cccGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTg                                                                                                                                  <  2:10639344/150‑1 (MQ=255)
                        gATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGTgg                                                                                                                               <  1:4235672/150‑1 (MQ=255)
                        gATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCgg                                                                                                                               <  2:2498694/150‑1 (MQ=255)
                        gATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCgg                                                                                                                               <  2:8827319/150‑1 (MQ=255)
                            cccGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTct                                                                                                                           >  2:5047153/1‑150 (MQ=255)
                            cccGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTct                                                                                                                           <  2:6939230/150‑1 (MQ=255)
                              cGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTctgc                                                                                                                         <  2:3263572/150‑1 (MQ=255)
                               ggAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTctgct                                                                                                                        <  2:6106168/150‑1 (MQ=255)
                               ggAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTctgct                                                                                                                        >  1:10964007/1‑150 (MQ=255)
                               ggAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATGGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTctgct                                                                                                                        >  2:11345604/1‑150 (MQ=255)
                               ggAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTAATTGGTAAGGTGCCCTCCAGCGCCGCAATTGCCAGATGAATCAGCGCCACTGCCCGTTTCAACGGCACACCAAACGCAATTAATGCAAGAGGGTctgct                                                                                                                        >  1:3370361/1‑150 (MQ=255)
                                  aaaGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTctgctgct                                                                                                                     <  1:7895320/150‑1 (MQ=255)
                                   aaGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTctgctgctg                                                                                                                    >  1:3074409/1‑150 (MQ=255)
                                     gTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGcc                                                                                                                  >  1:4811637/1‑150 (MQ=255)
                                      tttGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCa                                                                                                                 >  1:10005397/1‑150 (MQ=255)
                                      tttGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCa                                                                                                                 >  1:8809611/1‑150 (MQ=255)
                                      tttGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCa                                                                                                                 <  1:11949605/150‑1 (MQ=255)
                                      tttGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCAAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCa                                                                                                                 >  2:7284010/1‑150 (MQ=255)
                                          cAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCg                                                                                                             <  2:2397065/150‑1 (MQ=255)
                                          cAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCg                                                                                                             >  2:9478599/1‑150 (MQ=255)
                                           aGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCGGc                                                                                                            >  1:1264967/1‑150 (MQ=255)
                                              gTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGa                                                                                                         <  2:4889600/150‑1 (MQ=255)
                                               ttttCATAGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGACAGATGAATCAGCGCCTCTGCCCGTTACAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAg                                                                                                        >  2:3689570/1‑150 (MQ=255)
                                                 ttCATCGCCTGCTCCACATCGCCCGGTATTGTCAGTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTACCGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGc                                                                                                      <  2:8369879/150‑1 (MQ=255)
                                                    aTCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTc                                                                                                   <  1:587308/150‑1 (MQ=255)
                                                    aTCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTc                                                                                                   <  1:7098237/150‑1 (MQ=255)
                                                     tCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCg                                                                                                  <  2:7495464/150‑1 (MQ=255)
                                                      cGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCgg                                                                                                 <  1:4725983/150‑1 (MQ=255)
                                                      cGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCgg                                                                                                 <  1:5436208/150‑1 (MQ=255)
                                                      cGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCgg                                                                                                 <  1:5004673/150‑1 (MQ=255)
                                                         cTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGaa                                                                                              <  1:6617345/150‑1 (MQ=255)
                                                             tCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAg                                                                                          <  1:9780290/150‑1 (MQ=255)
                                                              ccACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGa                                                                                         <  1:3703827/150‑1 (MQ=255)
                                                              ccACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGa                                                                                         <  2:866090/150‑1 (MQ=255)
                                                               cacaTCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGat                                                                                        >  1:179937/1‑150 (MQ=255)
                                                                acaTCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGata                                                                                       >  1:2278608/1‑150 (MQ=255)
                                                                acaTCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGata                                                                                       >  1:8309660/1‑150 (MQ=255)
                                                                 caTCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGatat                                                                                      <  1:1647292/150‑1 (MQ=255)
                                                                     gCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTc                                                                                  <  1:2036204/150‑1 (MQ=255)
                                                                       ccGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCg                                                                                >  1:3937033/1‑150 (MQ=255)
                                                                           tATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGaaa                                                                            <  2:11523931/150‑1 (MQ=255)
                                                                                tCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAatcat                                                                       >  2:11221328/1‑150 (MQ=255)
                                                                                 cATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAatcatc                                                                      <  1:3606184/150‑1 (MQ=255)
                                                                                  aTTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGACGGGACGCAGATATACTCCGGAAAATCATcc                                                                     >  1:380508/1‑150 (MQ=255)
                                                                                   ttGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCa                                                                    <  2:4426927/150‑1 (MQ=255)
                                                                                    tGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAg                                                                   <  2:10392210/150‑1 (MQ=255)
                                                                                     ggTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGc                                                                  >  1:7241272/1‑150 (MQ=255)
                                                                                        aaGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGt                                                               >  2:656792/1‑150 (MQ=255)
                                                                                           gTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCg                                                            <  1:11354615/150‑1 (MQ=255)
                                                                                            tGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGc                                                           >  1:6259396/1‑150 (MQ=255)
                                                                                                cTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATa                                                       <  2:1644209/150‑1 (MQ=255)
                                                                                                cTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATa                                                       >  1:5239609/1‑150 (MQ=255)
                                                                                                  ccAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATaaa                                                     <  2:11925699/150‑1 (MQ=255)
                                                                                                   cAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAg                                                    >  2:5774850/1‑150 (MQ=255)
                                                                                                     gcgcCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGACAATCATCCAGCCGTTCGCCATAAAGCt                                                  >  2:1818237/1‑150 (MQ=255)
                                                                                                       gccgcATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTACCGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCtgt                                                <  2:1464658/150‑1 (MQ=255)
                                                                                                       gccgcATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCtgt                                                <  1:5577863/150‑1 (MQ=255)
                                                                                                        ccgcATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCtgtg                                               <  2:4224156/150‑1 (MQ=255)
                                                                                                        ccgcAATTGGCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATACAGCtgtg                                               <  2:3968093/150‑1 (MQ=255)
                                                                                                         cgcATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGc                                              <  1:2469539/150‑1 (MQ=255)
                                                                                                             tttGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCAct                                          >  1:8444491/1‑150 (MQ=255)
                                                                                                                gCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTg                                       <  2:10590054/150‑1 (MQ=255)
                                                                                                                 ccAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTACTGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTgg                                      <  2:12275930/150‑1 (MQ=255)
                                                                                                                 ccAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTgg                                      <  1:2781714/150‑1 (MQ=255)
                                                                                                                  cAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGc                                     <  2:211938/150‑1 (MQ=255)
                                                                                                                   aGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCg                                    <  1:9226269/150‑1 (MQ=255)
                                                                                                                   aGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCg                                    >  1:9397777/1‑150 (MQ=255)
                                                                                                                   aGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCg                                    >  2:9007048/1‑150 (MQ=255)
                                                                                                                    gATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCgg                                   <  2:8228058/150‑1 (MQ=255)
                                                                                                                          tCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAAtt                             <  1:10307321/150‑1 (MQ=255)
                                                                                                                           cAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAAttt                            <  2:7626672/150‑1 (MQ=255)
                                                                                                                               gcCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTgcc                        <  2:5264058/150‑1 (MQ=255)
                                                                                                                                 ctctGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTgccgc                      >  1:6763403/1‑150 (MQ=255)
                                                                                                                                  tctGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTgccgcc                     >  1:11104207/1‑150 (MQ=255)
                                                                                                                                  tctGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTgccgcc                     >  1:898675/1‑150 (MQ=255)
                                                                                                                                   ctGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCa                    <  1:11727571/150‑1 (MQ=255)
                                                                                                                                   ctGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCa                    >  2:10262010/1‑150 (MQ=255)
                                                                                                                                    tGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCAt                   >  2:10927241/1‑150 (MQ=255)
                                                                                                                                    tGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCAt                   >  2:7829826/1‑150 (MQ=255)
                                                                                                                                      cccGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCg                 >  2:9828274/1‑150 (MQ=255)
                                                                                                                                       ccGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGc                >  1:4440583/1‑150 (MQ=255)
                                                                                                                                       ccGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGc                >  2:5601505/1‑150 (MQ=255)
                                                                                                                                        cGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGcc               <  2:852885/150‑1 (MQ=255)
                                                                                                                                         gTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCa              <  1:3178348/150‑1 (MQ=255)
                                                                                                                                          tttCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCAc             >  1:9569799/1‑150 (MQ=255)
                                                                                                                                            tCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGc           <  2:4627596/150‑1 (MQ=255)
                                                                                                                                            tCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGc           >  2:9284210/1‑150 (MQ=255)
                                                                                                                                              aaCGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTc         <  2:3074409/150‑1 (MQ=255)
                                                                                                                                              aaCGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTc         >  2:3390855/1‑150 (MQ=255)
                                                                                                                                                 ggCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCAcc      >  2:6329318/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCACCAGTT  >  NC_000913/2146872‑2147170

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