Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,147,021 T→C 79.2% 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,0210TC79.2% 382.9 / 89.9 168M174V (ATG→GTG) alkA3‑methyl‑adenine DNA glycosylase II
Reads supporting (aligned to +/- strand):  ref base T (23/12);  new base C (49/84);  total (72/96)
Fisher's exact test for biased strand distribution p-value = 3.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 169 total aligned reads displayed.

AAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCACC  >  NC_000913/2146873‑2147166
                                                                                                                                                    |                                                                                                                                                 
aaaTAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATa                                                                                                                                                  >  1:1411500/1‑150 (MQ=255)
   taatTCGGCGTCCAGCGCCCGATTCCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCTCCGCTTTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCTACGGCAttcat                                                                                                                                               <  2:2338568/150‑5 (MQ=255)
   taatTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCt                                                                                                                                               >  2:1880716/1‑150 (MQ=255)
     atTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTaa                                                                                                                                             <  1:2501417/150‑1 (MQ=255)
      ttCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAAc                                                                                                                                            <  2:1334422/150‑1 (MQ=255)
      ttCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAAc                                                                                                                                            <  1:1032892/150‑1 (MQ=255)
      ttCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAAc                                                                                                                                            <  1:3085224/150‑1 (MQ=255)
       tCGCCGTCCAGCGCCCGATACCCGGAAAATTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACg                                                                                                                                           <  1:2116795/150‑1 (MQ=255)
        cGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGc                                                                                                                                          >  1:680431/1‑150 (MQ=255)
            gTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCtttt                                                                                                                                      <  2:2840641/150‑1 (MQ=255)
            gTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCtttt                                                                                                                                      >  2:1538032/1‑150 (MQ=255)
                aGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATg                                                                                                                                  <  1:76992/150‑1 (MQ=255)
                 gcgcCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGc                                                                                                                                 >  1:136104/1‑150 (MQ=255)
                   gcCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCt                                                                                                                               <  1:1819497/150‑1 (MQ=255)
                      cGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCg                                                                                                                            <  2:1761147/150‑1 (MQ=255)
                        aTACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCggg                                                                                                                          >  1:2992139/1‑150 (MQ=255)
                          aCCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTc                                                                                                                        >  1:3060497/1‑150 (MQ=255)
                          aCCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTc                                                                                                                        <  1:1368475/150‑1 (MQ=255)
                            ccGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTctg                                                                                                                      >  2:527046/1‑150 (MQ=255)
                            ccGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTctg                                                                                                                      <  1:790401/150‑1 (MQ=255)
                            ccGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTctg                                                                                                                      <  1:426810/150‑1 (MQ=255)
                              ggAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTctgct                                                                                                                    <  2:306241/150‑1 (MQ=255)
                               gAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTctgctg                                                                                                                   >  2:852569/1‑150 (MQ=255)
                                aaaaGTTTGCAGCGTTTTCATCGCCTGCTCCCCATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTctgctgc                                                                                                                  <  1:448297/150‑1 (MQ=255)
                                   aGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCtgctgctgc                                                                                                               <  1:359435/150‑1 (MQ=255)
                                    gTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTTATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGcc                                                                                                              >  1:1027744/1‑150 (MQ=255)
                                    gTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGcc                                                                                                              <  1:498546/150‑1 (MQ=255)
                                           gCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCt                                                                                                       <  2:3105723/150‑1 (MQ=255)
                                             gTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGa                                                                                                     <  2:22900/150‑1 (MQ=255)
                                             gTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGa                                                                                                     <  1:1602831/150‑1 (MQ=255)
                                             gTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGa                                                                                                     <  2:1074042/150‑1 (MQ=255)
                                             gTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGa                                                                                                     <  1:1797045/150‑1 (MQ=255)
                                              ttttCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAg                                                                                                    >  1:2909852/1‑150 (MQ=255)
                                                 tCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCg                                                                                                 >  2:521377/1‑150 (MQ=255)
                                                   aTCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTc                                                                                               <  1:904719/150‑1 (MQ=255)
                                                    tCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCg                                                                                              <  1:298870/150‑1 (MQ=255)
                                                    tCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCg                                                                                              >  2:1036327/1‑150 (MQ=255)
                                                    tCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCGGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCg                                                                                              >  2:2411144/1‑150 (MQ=255)
                                                       ccTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGa                                                                                           >  2:1494186/1‑150 (MQ=255)
                                                         tGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAg                                                                                         <  1:1050517/150‑1 (MQ=255)
                                                          gCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGc                                                                                        <  1:313924/150‑1 (MQ=255)
                                                          gCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGc                                                                                        <  1:2018593/150‑1 (MQ=255)
                                                           cTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCa                                                                                       >  1:1342441/1‑150 (MQ=255)
                                                           cTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCa                                                                                       >  2:718680/1‑150 (MQ=255)
                                                             ccACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCCGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTACCGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGa                                                                                     <  2:1149446/150‑1 (MQ=255)
                                                             ccACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGa                                                                                     <  2:1237194/150‑1 (MQ=255)
                                                               acaTCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGata                                                                                   <  1:795877/150‑1 (MQ=255)
                                                                 aTCGCCCGGTATTGGGATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTCACGCTTTTACTGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGAtata                                                                                 >  2:896115/1‑150 (MQ=255)
                                                                    gCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGTCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTc                                                                              <  1:657627/150‑1 (MQ=255)
                                                                    gCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTc                                                                              <  2:252001/150‑1 (MQ=255)
                                                                    gCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTc                                                                              <  2:252007/150‑1 (MQ=255)
                                                                     cccGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTcc                                                                             >  1:1610936/1‑150 (MQ=255)
                                                                            ttGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAat                                                                      >  2:2647980/1‑150 (MQ=255)
                                                                               tCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAatcat                                                                   >  2:2118545/1‑150 (MQ=255)
                                                                               tCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAatcat                                                                   <  1:3105093/150‑1 (MQ=255)
                                                                               tCAGTGGTAAGGGGCCCTCCCGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTACCGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAatcat                                                                   <  2:176810/150‑1 (MQ=255)
                                                                                 aTTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTACCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATcc                                                                 <  2:577141/150‑1 (MQ=255)
                                                                                    ggTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGc                                                              >  2:985031/1‑150 (MQ=255)
                                                                                     gTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGcc                                                             >  1:2239112/1‑150 (MQ=255)
                                                                                           tGCCCTCCAGCGCCGCATATGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACACCGGAAAATCATCCAGCCGTTCGc                                                       >  1:3069621/1‑150 (MQ=255)
                                                                                               cTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATa                                                   >  1:678358/1‑150 (MQ=255)
                                                                                                tCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGTCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATaa                                                  >  1:3007143/1‑150 (MQ=255)
                                                                                                tCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATaa                                                  <  1:1132725/150‑1 (MQ=255)
                                                                                                tCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATaa                                                  <  1:2681618/150‑1 (MQ=255)
                                                                                                  cAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAg                                                <  1:783282/150‑1 (MQ=255)
                                                                                                  cAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAg                                                <  1:1823802/150‑1 (MQ=255)
                                                                                                  cAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAg                                                <  2:2677707/150‑1 (MQ=255)
                                                                                                  cAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAg                                                <  2:2677705/150‑1 (MQ=255)
                                                                                                  cAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAg                                                <  2:939302/150‑1 (MQ=255)
                                                                                                  cAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAg                                                <  1:2268308/150‑1 (MQ=255)
                                                                                                    gcgcCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCt                                              <  2:2909683/150‑1 (MQ=255)
                                                                                                       ccgcATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCtgtg                                           >  1:2820133/1‑150 (MQ=255)
                                                                                                               gCCAGATGAATCAGCGCCTCTGCCCGTTTCAGCGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTg                                   >  2:1449537/1‑150 (MQ=255)
                                                                                                                 cAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGTAGCAGATATAGTCCGGTATTTCATCCAGCGGTTCGCGTTAAAGCTGTGCCACTCTggg                                 >  1:2138939/1‑149 (MQ=255)
                                                                                                                 cAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGc                                 >  1:1765697/1‑150 (MQ=255)
                                                                                                                 cAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGc                                 <  2:2984644/150‑1 (MQ=255)
                                                                                                                  aGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCg                                <  2:332687/150‑1 (MQ=255)
                                                                                                                  aGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCg                                >  1:1829561/1‑150 (MQ=255)
                                                                                                                   gATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCgg                               >  2:121052/1‑150 (MQ=255)
                                                                                                                     tGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTc                             >  1:1737905/1‑150 (MQ=255)
                                                                                                                      gAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCa                            <  2:671941/150‑1 (MQ=255)
                                                                                                                       aaTCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCaa                           <  2:1339052/150‑1 (MQ=255)
                                                                                                                        aTCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAAt                          <  2:1885929/150‑1 (MQ=255)
                                                                                                                        aTCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAAt                          <  2:2722517/150‑1 (MQ=255)
                                                                                                                         tCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAAtt                         <  1:2213307/150‑1 (MQ=255)
                                                                                                                         tCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAAtt                         >  2:2858927/1‑150 (MQ=255)
                                                                                                                           aGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAAtttt                       >  1:2536753/1‑150 (MQ=255)
                                                                                                                             cgcCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTgc                     <  1:2227370/150‑1 (MQ=255)
                                                                                                                                  ctGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCa                <  1:1193138/150‑1 (MQ=255)
                                                                                                                                   tGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCAt               <  1:2014243/150‑1 (MQ=255)
                                                                                                                                     cccGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCACATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCAGGCGGTCAAATTAGCCCCCATCg             >  2:2014984/1‑150 (MQ=255)
                                                                                                                                     cccGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCg             >  1:1796246/1‑150 (MQ=255)
                                                                                                                                     cccGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCg             >  2:249433/1‑150 (MQ=255)
                                                                                                                                     cccGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCg             >  1:1864946/1‑150 (MQ=255)
                                                                                                                                     cccGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCg             >  1:889465/1‑150 (MQ=255)
                                                                                                                                        gTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCa          >  2:1928317/1‑150 (MQ=255)
                                                                                                                                        gTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCa          <  2:1952454/150‑1 (MQ=255)
                                                                                                                                            cAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCt      >  2:1642215/1‑150 (MQ=255)
                                                                                                                                               cGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCAc   <  2:769878/150‑1 (MQ=255)
                                                                                                                                                ggCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCACCATCGCCACGCTCAcc  >  2:2606651/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCACC  >  NC_000913/2146873‑2147166

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: