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% 1563.5 / NA 404intergenic (+83/+15)fabI/PA14_RS16690enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (231/172);  minor base A (0/1);  total (231/173)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0084633,673,5350CT97.1% 1662.4 / 18.2 415intergenic (+84/+14)fabI/PA14_RS16690enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (6/6);  new base T (231/169);  total (237/178)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01
*NC_0084633,673,5360TC97.4% 1545.5 / 22.8 415intergenic (+85/+13)fabI/PA14_RS16690enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (6/5);  new base C (231/173);  total (237/178)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 415 total aligned reads displayed.

ATCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGNGTAAACGGCCTTTTTCTTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTC  >  NC_008463/3673392‑3673680
                                                                                                                                              |||                                                                                                                                                 
aTCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCAtc                                                                                                                                            >  1:10754187/1‑151 (MQ=255)
 tCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCAtct                                                                                                                                           >  1:5643305/1‑151 (MQ=255)
    tctACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCgg                                                                                                                                        >  1:4259612/1‑151 (MQ=255)
      tACGTCGATGGCGGCTACAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCTTGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGtc                                                                                                                                      <  1:500296/151‑1 (MQ=255)
           cGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGt                                                                                                                                 >  1:6852337/1‑151 (MQ=255)
               ggcggcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTAttt                                                                                                                             >  2:1405079/1‑151 (MQ=255)
                 cggcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTtct                                                                                                                           <  2:6645649/151‑1 (MQ=255)
                  ggcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTtctc                                                                                                                          <  2:8388275/151‑1 (MQ=255)
                   gcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCg                                                                                                                         <  1:10296910/151‑1 (MQ=255)
                   gcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCg                                                                                                                         <  2:249577/151‑1 (MQ=255)
                   gcTCCAACACTACCGCCTTGGGCCCGCTGGACGACGACTTATCACCCCGCTCATCAAAATTGGCCGCTACGAGCGGCATTTTTCTTGGGCGTAATAAAGAAACGGCCGCGTACGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCg                                                                                                                         <  1:7257418/151‑1 (MQ=255)
                          cacTACCGCCCTGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCGTCGAGCGGCCTTTTTCAGGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTAGTTCTCGACTTCCg                                                                                                                  >  2:2412456/1‑151 (MQ=255)
                          cacTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCg                                                                                                                  >  2:10095837/1‑151 (MQ=255)
                             aaCCGCCATGGGCCCGCTGTACGACGAATGAACACCCCGCTATTTAAAAATGGCCGATTCGAGCGGCCTTTTTCATGGGCGGAATAACCAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCt                                                                                                               <  1:874552/150‑1 (MQ=255)
                                  ccttGGGCCCGCTGTAAGACTACTTATAACCCCGCTCATTAAAATTGGCGGCTTCGAGCGGCCTTTTTCTTGGGCGTAATAAAGAAAAGGCCGCTTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTcc                                                                                                          <  1:8558921/148‑1 (MQ=255)
                                    aTGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTg                                                                                                        >  1:88692/1‑151 (MQ=255)
                                    aTGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTg                                                                                                        >  2:2226637/1‑151 (MQ=255)
                                    aTGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTg                                                                                                        <  2:9154492/151‑1 (MQ=255)
                                       ggCCCGCTGGACGACGACTTATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAAAAAAGACAAGGACGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCaa                                                                                                     <  2:6120244/151‑1 (MQ=255)
                                       ggCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCaa                                                                                                     <  1:7673916/151‑1 (MQ=255)
                                         cccGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATACAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACt                                                                                                   <  2:6231751/151‑1 (MQ=255)
                                          ccGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTg                                                                                                  >  2:1474788/1‑151 (MQ=255)
                                           cGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCCACTGa                                                                                                 >  1:4729488/1‑151 (MQ=255)
                                            gCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAc                                                                                                >  1:6059649/1‑151 (MQ=255)
                                             cTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcg                                                                                               <  2:5736623/151‑1 (MQ=255)
                                               ggACGACTACTTATCACCCCGCACCTGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAACGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcggc                                                                                             <  2:6757268/151‑1 (MQ=255)
                                                gacgacgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcggcg                                                                                            >  1:8296678/1‑151 (MQ=255)
                                                  cgacgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATCAAGAACAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGa                                                                                          <  2:7380738/151‑1 (MQ=255)
                                                  cgacgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGa                                                                                          >  2:8210207/1‑151 (MQ=255)
                                                   gacgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGaa                                                                                         >  1:7468229/1‑151 (MQ=255)
                                                     cgacTGAACACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAACGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAgg                                                                                       <  2:7262789/151‑1 (MQ=255)
                                                       acTGAACACCCCGCACATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAAAAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAggcg                                                                                     <  1:2004709/151‑1 (MQ=255)
                                                           aTCACCCCGCTCATTAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGa                                                                                 <  1:2260821/151‑1 (MQ=255)
                                                             cACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAg                                                                               >  2:10758142/1‑151 (MQ=255)
                                                              aCCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGACAACGCCGCGTACGCGGCCATTTCTCTGCCTCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGt                                                                              <  2:8869336/151‑1 (MQ=255)
                                                                      cATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGt                                                                      >  2:4797325/1‑151 (MQ=255)
                                                                        tGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTcc                                                                    >  1:6579031/1‑151 (MQ=255)
                                                                            aaaTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTg                                                                >  2:2443894/1‑151 (MQ=255)
                                                                                ggCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGcgcg                                                            <  1:1121897/151‑1 (MQ=255)
                                                                                ggCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGcgcg                                                            >  2:296508/1‑151 (MQ=255)
                                                                                   cGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAg                                                         >  2:3473204/1‑151 (MQ=255)
                                                                                     cTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGcc                                                       >  1:4821127/1‑151 (MQ=255)
                                                                                      ttCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCa                                                      >  2:4250921/1‑151 (MQ=255)
                                                                                       tCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAg                                                     >  1:2937984/1‑151 (MQ=255)
                                                                                          aGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGaa                                                  >  2:2085760/1‑151 (MQ=255)
                                                                                          aGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGaa                                                  >  1:3322253/1‑151 (MQ=255)
                                                                                          aGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGaa                                                  >  1:2486451/1‑151 (MQ=255)
                                                                                           gCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAg                                                 <  1:6218966/151‑1 (MQ=255)
                                                                                            cGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgc                                                >  1:9579549/1‑151 (MQ=255)
                                                                                            cGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgc                                                >  1:4904260/1‑151 (MQ=255)
                                                                                            cGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgc                                                <  1:2123505/151‑1 (MQ=255)
                                                                                            cGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTACCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGACTGCCTGTCCGCTGCGCGAAGCCAGGAAgc                                                >  1:7936089/1‑151 (MQ=255)
                                                                                             ggCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcg                                               >  1:917453/1‑151 (MQ=255)
                                                                                               ccTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcggc                                             >  2:3741398/1‑151 (MQ=255)
                                                                                                cTTTTTCTTGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcggcg                                            <  1:3123419/151‑1 (MQ=255)
                                                                                                cTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcggcg                                            >  1:105743/1‑151 (MQ=255)
                                                                                                 tttttCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGcggcgg                                           >  2:7939370/1‑151 (MQ=255)
                                                                                                  ttttCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGt                                          >  1:1323178/1‑151 (MQ=255)
                                                                                                  ttttCATGGGCGGAAACAAGAAAAGGCCGCGTACGCGGCCTTTTCTCTGCCTCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGt                                          <  1:9829597/151‑1 (MQ=255)
                                                                                                    ttCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTAc                                        <  2:8642594/151‑1 (MQ=255)
                                                                                                     tCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACa                                       <  2:3996135/151‑1 (MQ=255)
                                                                                                      cATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACAt                                      >  2:10274044/1‑151 (MQ=255)
                                                                                                      cATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACAt                                      >  2:5005143/1‑151 (MQ=255)
                                                                                                      cATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCGGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACAt                                      >  1:456833/1‑151 (MQ=255)
                                                                                                       aTGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGACGCGGCGGTACATg                                     >  2:6444279/1‑151 (MQ=255)
                                                                                                         gggCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGc                                   <  1:1652131/151‑1 (MQ=255)
                                                                                                         gggCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGc                                   >  1:5083529/1‑151 (MQ=255)
                                                                                                         gggCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGc                                   >  2:1337600/1‑151 (MQ=255)
                                                                                                          ggCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGcc                                  >  1:4174527/1‑151 (MQ=255)
                                                                                                           gCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCt                                 <  2:7750681/151‑1 (MQ=255)
                                                                                                           gCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCt                                 <  2:7637117/151‑1 (MQ=255)
                                                                                                               aaTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtctc                             >  2:9324914/1‑151 (MQ=255)
                                                                                                                aTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcg                            >  2:5174179/1‑151 (MQ=255)
                                                                                                                aTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCGTCtcg                            >  1:8115719/1‑151 (MQ=255)
                                                                                                                 tAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgt                           >  2:8066670/1‑151 (MQ=255)
                                                                                                                   aaGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTACtcgtcg                         >  2:1781516/1‑151 (MQ=255)
                                                                                                                      aaaaGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACTTGGCCTTCtcgtcgtcg                      >  2:5827680/1‑151 (MQ=255)
                                                                                                                      aaaaGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtcg                      >  2:9637200/1‑151 (MQ=255)
                                                                                                                        aaGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGga                    >  1:3853726/1‑151 (MQ=255)
                                                                                                                         aGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGgag                   <  2:7025323/151‑1 (MQ=255)
                                                                                                                         aGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGgag                   >  1:9197122/1‑151 (MQ=255)
                                                                                                                           gCCGCTTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGAt                 <  1:3705213/151‑1 (MQ=255)
                                                                                                                            ccGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATa                >  2:9348320/1‑151 (MQ=255)
                                                                                                                             cgcgTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGACGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAg               >  2:78991/1‑151 (MQ=255)
                                                                                                                             cgcgTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAg               >  1:8201783/1‑151 (MQ=255)
                                                                                                                              gCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGc              >  2:10517109/1‑151 (MQ=255)
                                                                                                                              gCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGc              >  1:1757298/1‑151 (MQ=255)
                                                                                                                              gCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGc              <  2:548811/151‑1 (MQ=255)
                                                                                                                               cgTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcc             >  1:5756762/2‑151 (MQ=255)
                                                                                                                               cgTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcc             >  1:567395/2‑151 (MQ=255)
                                                                                                                                gTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcct            >  2:3249979/1‑151 (MQ=255)
                                                                                                                                gTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcct            <  1:5138543/151‑1 (MQ=255)
                                                                                                                                gTAAGCGGCCCTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcct            >  2:903632/1‑151 (MQ=255)
                                                                                                                                  aagcGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctcc          <  1:9430569/148‑1 (MQ=255)
                                                                                                                                  aagcGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctcc          <  1:1007390/148‑1 (MQ=255)
                                                                                                                                     cGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCtcctcc       >  1:7771356/1‑151 (MQ=255)
                                                                                                                                      ggCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCTTCGTCGGAGATAGCCTCCTCCg      >  2:7048530/1‑151 (MQ=255)
                                                                                                                                      ggCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCg      <  1:5410687/151‑1 (MQ=255)
                                                                                                                                      ggCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGAACATGGCCTACTCGTCGTCGGAGATAGCCTCCTCCg      >  2:4625398/1‑151 (MQ=255)
                                                                                                                                           ttttCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTc  >  2:10016905/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: