Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,673,534 3 bp→CTC intergenic (+83/+13) fabI → / ← PA14_RS16690 enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,673,5340TC99.8% 1749.9 / ‑4.2 439intergenic (+83/+15)fabI/PA14_RS16690enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (226/212);  total (227/212)
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.94e-01
*NC_0084633,673,5350CT95.8% 1880.4 / 37.8 457intergenic (+84/+14)fabI/PA14_RS16690enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (11/8);  new base T (226/212);  total (237/220)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01
*NC_0084633,673,5360TC95.8% 1677.6 / 49.4 457intergenic (+85/+13)fabI/PA14_RS16690enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (11/8);  new base C (226/212);  total (237/220)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 457 total aligned reads displayed.

ATCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGNGTAAACGGCCTTTTTCTTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTC  >  NC_008463/3673392‑3673680
                                                                                                                                              |||                                                                                                                                                 
 tCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCAtct                                                                                                                                           <  2:1207883/151‑1 (MQ=255)
  ccTCTACGTCGAGGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCAtctc                                                                                                                                          <  2:9743490/151‑1 (MQ=255)
   ctctACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCg                                                                                                                                         <  2:324025/151‑1 (MQ=255)
      tACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGACTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCGCGGtc                                                                                                                                      >  2:11405975/1‑151 (MQ=255)
       aCGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTAATGAAAAATGGCCGCTTCGAGCGGCCATTTTCATGGGCGGAATAACGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGtct                                                                                                                                     <  2:10178320/151‑1 (MQ=255)
         gTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCa                                                                                                                                   <  2:896050/151‑1 (MQ=255)
            gATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGACTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTa                                                                                                                                >  1:348989/1‑151 (MQ=255)
              tGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTAtt                                                                                                                              >  1:6080598/1‑151 (MQ=255)
                gcggcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTtc                                                                                                                            >  2:1103138/1‑151 (MQ=255)
                   gcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCg                                                                                                                         >  1:5272180/1‑151 (MQ=255)
                    cTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGa                                                                                                                        <  1:6682249/151‑1 (MQ=255)
                     ttCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGAc                                                                                                                       <  1:9406064/151‑1 (MQ=255)
                       cAACACTACCGCCCTGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACtt                                                                                                                     >  2:6163506/1‑151 (MQ=255)
                         acacTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAATGGCCGCGTAAGCGGCCTTTTCTCTGGAGCTCGGTCTCAGTATTTCTCGACTTcc                                                                                                                   >  2:8572232/1‑151 (MQ=255)
                            cTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCGTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGcc                                                                                                                <  1:3548932/151‑1 (MQ=255)
                            cTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGATCATTAAAATTGGCCGCTTCGAGCGGCATTTTTCATGGGCGTAATAAACAAACGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGcc                                                                                                                <  1:3029597/151‑1 (MQ=255)
                             tACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCt                                                                                                               >  2:253949/1‑151 (MQ=255)
                             tACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCt                                                                                                               <  2:9665586/151‑1 (MQ=255)
                              aCCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAACAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCtt                                                                                                              <  1:9193307/151‑1 (MQ=255)
                                cGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATACAGAACAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGt                                                                                                            <  1:7574109/151‑1 (MQ=255)
                                 gcctggggCCCGCTGGACGACGCCTGAACACCCCGCACATGAAAAAGGGCCGCTGCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAACGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTc                                                                                                           <  1:10781119/146‑1 (MQ=255)
                                   cATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTGTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCt                                                                                                         >  2:1485201/1‑151 (MQ=255)
                                    aTGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTg                                                                                                        >  2:3606864/1‑151 (MQ=255)
                                       ggCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCaa                                                                                                     <  2:8093592/151‑1 (MQ=255)
                                         cccGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACt                                                                                                   >  2:7971955/1‑151 (MQ=255)
                                            gCTGTACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTAATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAc                                                                                                <  2:6400919/151‑1 (MQ=255)
                                             cTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcg                                                                                               >  1:2366492/1‑151 (MQ=255)
                                             cTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcg                                                                                               >  2:5572986/1‑151 (MQ=255)
                                             cTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcg                                                                                               >  1:6901901/1‑151 (MQ=255)
                                             cTGGACGACGACTGATCACCCCGATCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGTAATAAAGAAACGGCCGCGTACGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcg                                                                                               <  1:10847491/151‑1 (MQ=255)
                                             cTGGACGACGACTGAACACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcg                                                                                               <  1:6486937/151‑1 (MQ=255)
                                              tGGACGAGGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGAAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTATCGACTTCCGCCTTGTCCTGCAACTGAcgg                                                                                              >  1:2554488/1‑151 (MQ=255)
                                               ggACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTACGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcggc                                                                                             <  2:2439720/151‑1 (MQ=255)
                                                    acgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAg                                                                                        >  1:8702995/1‑151 (MQ=255)
                                                      gacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAggc                                                                                      <  2:6269743/151‑1 (MQ=255)
                                                       acTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAggcg                                                                                     >  1:11590178/1‑151 (MQ=255)
                                                         tGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCATTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTACGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAggcggc                                                                                   <  1:1495984/151‑1 (MQ=255)
                                                            tCACCCCGCCCCTGAAAAATGGCCGCTTCGAGCGGCCTTTTTCTTGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCATTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGaa                                                                                <  1:9744304/151‑1 (MQ=255)
                                                                     tCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTg                                                                       <  2:4681951/151‑1 (MQ=255)
                                                                      cATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGt                                                                      >  2:6568326/1‑151 (MQ=255)
                                                                           aaaaTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCt                                                                 <  1:4080040/151‑1 (MQ=255)
                                                                           aaaaTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCt                                                                 >  1:4194613/1‑151 (MQ=255)
                                                                                  ccGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGaa                                                          <  2:5996881/151‑1 (MQ=255)
                                                                                   cGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAg                                                         <  2:2936518/151‑1 (MQ=255)
                                                                                   cGCTTCGAGCGGCCATTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAg                                                         <  2:2815557/151‑1 (MQ=255)
                                                                                    gCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGc                                                        >  1:1612716/1‑151 (MQ=255)
                                                                                     cTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGcc                                                       <  2:3388921/151‑1 (MQ=255)
                                                                                      ttCGAGCGGCCTTTTTCATGGGCGGAAGAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCa                                                      <  2:3759966/151‑1 (MQ=255)
                                                                                        cGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAgg                                                    >  2:108654/1‑151 (MQ=255)
                                                                                          aGCGGCCTTTTTCATGGGCGTAATACAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGaa                                                  <  1:4649245/151‑1 (MQ=255)
                                                                                          aGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGTCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGaa                                                  <  2:3088794/151‑1 (MQ=255)
                                                                                          aGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGaa                                                  >  1:7702436/1‑151 (MQ=255)
                                                                                          aGCGGCCTATTACATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGaa                                                  <  2:1972259/151‑1 (MQ=255)
                                                                                          aGCGGCCATTTTCATGGGCGGAATAAAGACAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGaa                                                  <  1:616567/151‑1 (MQ=255)
                                                                                            cGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgc                                                >  1:1377600/1‑151 (MQ=255)
                                                                                                cTTTTTCAGGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcggcg                                            <  2:7676793/151‑1 (MQ=255)
                                                                                                    ttCATGGGCGTAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTAc                                        <  1:11252030/151‑1 (MQ=255)
                                                                                                     tCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCCTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACa                                       >  1:7955602/1‑151 (MQ=255)
                                                                                                       aTGGGCGGAATACAGCAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATg                                     <  2:9893245/151‑1 (MQ=255)
                                                                                                        tGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATgg                                    >  2:9149512/1‑151 (MQ=255)
                                                                                                           gCGGACTAACGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCt                                 <  2:11581884/151‑1 (MQ=255)
                                                                                                           gCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCt                                 <  2:1596229/151‑1 (MQ=255)
                                                                                                            cGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCtt                                >  1:87494/1‑151 (MQ=255)
                                                                                                            cGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCtt                                >  1:10649020/1‑151 (MQ=255)
                                                                                                             ggAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtc                               >  1:10606111/1‑151 (MQ=255)
                                                                                                              gAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtct                              <  1:8149516/151‑1 (MQ=255)
                                                                                                               aaTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtctc                             >  1:5504037/1‑151 (MQ=255)
                                                                                                               aaTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtctc                             >  1:5154162/1‑151 (MQ=255)
                                                                                                                aTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcg                            >  2:574621/1‑151 (MQ=255)
                                                                                                                aTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcg                            >  1:8267092/1‑151 (MQ=255)
                                                                                                                aTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcg                            <  1:440515/151‑1 (MQ=255)
                                                                                                                aTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCAGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcg                            <  1:7406357/151‑1 (MQ=255)
                                                                                                                 tAAACAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgt                           <  1:193042/151‑1 (MQ=255)
                                                                                                                  aaaGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtc                          <  1:9557232/151‑1 (MQ=255)
                                                                                                                  aaaGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtc                          >  1:8421076/1‑151 (MQ=255)
                                                                                                                   aaGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcg                         >  2:7324410/1‑151 (MQ=255)
                                                                                                                    aGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgt                        >  1:11083199/1‑151 (MQ=255)
                                                                                                                     gAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtc                       <  2:4354172/151‑1 (MQ=255)
                                                                                                                      aaaaGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtcg                      >  1:10370011/1‑151 (MQ=255)
                                                                                                                        aaGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGga                    >  2:7187701/1‑151 (MQ=255)
                                                                                                                         aGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGgag                   >  1:7236793/1‑151 (MQ=255)
                                                                                                                          ggCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGgaga                  >  1:11459404/1‑151 (MQ=255)
                                                                                                                           gCCGCGTACGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGAt                 <  1:8949623/151‑1 (MQ=255)
                                                                                                                            ccGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATa                >  2:5118982/1‑151 (MQ=255)
                                                                                                                            ccGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATa                <  1:7843717/151‑1 (MQ=255)
                                                                                                                             cgcgTAAGCGGCCTTTTCTCTGCCTCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAg               <  2:8425838/151‑1 (MQ=255)
                                                                                                                               cgTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcc             <  1:5757816/150‑1 (MQ=255)
                                                                                                                                 tAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctc           <  2:1039072/151‑1 (MQ=255)
                                                                                                                                 tAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctc           >  1:10037894/1‑151 (MQ=255)
                                                                                                                                 tAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctc           >  1:2850917/1‑151 (MQ=255)
                                                                                                                                     cGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCtcctcc       >  2:4500221/1‑151 (MQ=255)
                                                                                                                                     cGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCtcctcc       >  2:1715204/1‑151 (MQ=255)
                                                                                                                                      ggCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCg      >  1:10965120/1‑151 (MQ=255)
                                                                                                                                      ggCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCg      >  2:10552398/1‑151 (MQ=255)
                                                                                                                                      ggCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTCCTCGTCGTCGGAGATAGCCTCCTCCg      >  1:11690644/1‑151 (MQ=255)
                                                                                                                                       gCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCgg     <  1:3771817/151‑1 (MQ=255)
                                                                                                                                           ttttCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTc  <  2:3982080/151‑1 (MQ=255)
                                                                                                                                           ttttCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTc  <  2:2314747/151‑1 (MQ=255)
                                                                                                                                           ttttCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTc  >  2:6962044/1‑151 (MQ=255)
                                                                                                                                           ttttCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTc  >  2:4933895/1‑151 (MQ=255)
                                                                                                                                              |||                                                                                                                                                 
ATCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGNGTAAACGGCCTTTTTCTTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTC  >  NC_008463/3673392‑3673680

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