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% 1668.8 / ‑6.2 437intergenic (+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 (261/175);  minor base G (1/0);  total (262/175)
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.99e-01
*NC_0084633,673,5350CT95.8% 1845.7 / 37.7 456intergenic (+84/+14)fabI/PA14_RS16690enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (18/1);  new base T (261/175);  total (280/176)
Fisher's exact test for biased strand distribution p-value = 1.32e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01
*NC_0084633,673,5360TC95.8% 1659.2 / 49.4 456intergenic (+85/+13)fabI/PA14_RS16690enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (18/1);  new base C (262/175);  total (280/176)
Fisher's exact test for biased strand distribution p-value = 1.33e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 456 total aligned reads displayed.

ATCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGNGTAAACGGCCTTTTTCTTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTC  >  NC_008463/3673392‑3673680
                                                                                                                                              |||                                                                                                                                                 
    tctACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGTAATAACGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCgg                                                                                                                                        <  1:7229420/151‑1 (MQ=255)
    tctACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCgg                                                                                                                                        <  2:9365125/151‑1 (MQ=255)
       aCGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGtct                                                                                                                                     <  2:438075/151‑1 (MQ=255)
         gTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCa                                                                                                                                   >  2:2725703/1‑151 (MQ=255)
          tCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAg                                                                                                                                  <  2:11162442/151‑1 (MQ=255)
             aTGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTAt                                                                                                                               >  2:7359942/1‑151 (MQ=255)
                  ggcTTCAACACGACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTtctc                                                                                                                          <  2:8895282/151‑1 (MQ=255)
                   gcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCg                                                                                                                         <  1:3579367/151‑1 (MQ=255)
                   gatacAACACTACCGCCTTGGGCCCGCTGGACGACTACTTATCACCCCGAGCATGAAAATGGGCCGCTTCGAGGGGCCTTTTTCATGGGGGGAATAAACAAAAGGCCGCGTACGGGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCg                                                                                                                         <  1:1165263/147‑1 (MQ=255)
                    cTTCAACACTCCCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCGTCGAGCGGCCTTTTTCATGGGCGGAATAACGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGa                                                                                                                        <  2:3791327/151‑1 (MQ=255)
                     tacAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGTAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGAc                                                                                                                       <  1:4150044/149‑1 (MQ=255)
                       cAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTCTTTCTCGACtt                                                                                                                     >  1:9677495/1‑151 (MQ=255)
                        aaCACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAACAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTc                                                                                                                    <  1:10849581/151‑1 (MQ=255)
                         acacTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATTAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTcc                                                                                                                   <  1:2639881/151‑1 (MQ=255)
                             tACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCt                                                                                                               >  1:1485263/1‑151 (MQ=255)
                              aCCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCtt                                                                                                              <  1:10608799/151‑1 (MQ=255)
                                cGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGt                                                                                                            >  1:3344883/1‑151 (MQ=255)
                                   cATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATTAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGTAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCt                                                                                                         <  1:878792/151‑1 (MQ=255)
                                   cATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCt                                                                                                         >  2:2561269/1‑151 (MQ=255)
                                   cATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCt                                                                                                         >  2:6573127/1‑151 (MQ=255)
                                    aTGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTg                                                                                                        <  2:8635482/151‑1 (MQ=255)
                                     tGGGCCCGCTGGACGACGACGGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGACAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGc                                                                                                       <  1:6540376/151‑1 (MQ=255)
                                            gCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAc                                                                                                <  1:3515758/151‑1 (MQ=255)
                                             cTGGACGCCGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATACAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcg                                                                                               <  1:7719118/151‑1 (MQ=255)
                                               ggACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcggc                                                                                             >  1:1699239/1‑151 (MQ=255)
                                                gacgacgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAACGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcggcg                                                                                            <  1:5063302/151‑1 (MQ=255)
                                                    acgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAg                                                                                        >  2:11002063/1‑151 (MQ=255)
                                                    acgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAg                                                                                        <  2:10781946/151‑1 (MQ=255)
                                                    acgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGAGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCTGTATTTGTCGAACTCCGCCTTGTCCTGCACCTGACGGCGGAAg                                                                                        >  1:7482863/1‑151 (MQ=255)
                                                     ctacTGATCACCCCGCTCATTAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAgg                                                                                       <  2:9649262/149‑1 (MQ=255)
                                                         tGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAggcggc                                                                                   <  2:3248223/151‑1 (MQ=255)
                                                         tGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAggcggc                                                                                   >  1:2322758/1‑151 (MQ=255)
                                                         tGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAggcggc                                                                                   >  1:7936312/1‑151 (MQ=255)
                                                         tGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAggcggc                                                                                   >  1:4253842/1‑151 (MQ=255)
                                                             cACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAg                                                                               >  1:9998167/1‑151 (MQ=255)
                                                               ccccGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGGCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGCAGTc                                                                             >  2:10516309/1‑151 (MQ=255)
                                                                cccGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGTAGGCGGCGAAGTCt                                                                            >  1:9790521/1‑151 (MQ=255)
                                                                cccGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCt                                                                            <  1:10010364/151‑1 (MQ=255)
                                                                    cTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGTAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCt                                                                        <  1:3941034/151‑1 (MQ=255)
                                                                     tCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTg                                                                       >  2:9265912/1‑151 (MQ=255)
                                                                     tCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTg                                                                       >  2:3612974/1‑151 (MQ=255)
                                                                      cATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGt                                                                      >  1:10697334/1‑151 (MQ=255)
                                                                        tGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTcc                                                                    >  1:6359311/1‑151 (MQ=255)
                                                                         gAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCg                                                                   <  2:10338252/151‑1 (MQ=255)
                                                                              aTGGCCGCTTCGAGCGGCCTTTGTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCTGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTgcg                                                              >  2:4022646/1‑151 (MQ=255)
                                                                               tGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGACACGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTACCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTgcgc                                                             <  1:5281360/151‑1 (MQ=255)
                                                                               tGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTgcgc                                                             >  1:5977017/1‑151 (MQ=255)
                                                                               tGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTgcgc                                                             >  1:5568088/1‑151 (MQ=255)
                                                                                   cGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAg                                                         >  1:7304112/1‑151 (MQ=255)
                                                                                   cGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAg                                                         >  2:2164397/1‑151 (MQ=255)
                                                                                   cGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAg                                                         >  2:1323094/1‑151 (MQ=255)
                                                                                   cGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCGGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAg                                                         >  1:3150291/1‑151 (MQ=255)
                                                                                       tCGAGCGGCCTTTTTCATGGGCGGAATAAAGAACAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAg                                                     <  2:4634069/151‑1 (MQ=255)
                                                                                       tCGAGCGGCCTATTTCATGGGCGGAATAAAGACAAGGCCGCGTACGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAg                                                     <  1:10899541/151‑1 (MQ=255)
                                                                                           gCGGCCTTTTTCATGGGCGGAATACAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAg                                                 <  1:6865074/151‑1 (MQ=255)
                                                                                                cTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcggcg                                            >  1:3412669/1‑151 (MQ=255)
                                                                                                cTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCAAGGAAgcggcg                                            >  1:8196166/1‑151 (MQ=255)
                                                                                                   tttCATGGGCGGACTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTa                                         <  2:6788246/151‑1 (MQ=255)
                                                                                                   tttCATGGGCGGAATACAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTa                                         <  1:11167929/151‑1 (MQ=255)
                                                                                                   tttCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGACGCCAGGAAGCGGCGGTa                                         >  2:2276896/1‑151 (MQ=255)
                                                                                                   attCATGGGCGCAATAACGAAAAGGCCGCGTAAGCGGCATATTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGACCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTa                                         <  2:5547984/150‑1 (MQ=255)
                                                                                                    ttCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTAc                                        >  1:2718400/1‑151 (MQ=255)
                                                                                                    ttCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTAc                                        >  1:10104417/1‑151 (MQ=255)
                                                                                                      cATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACAt                                      >  1:440442/1‑151 (MQ=255)
                                                                                                        tGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATgg                                    >  2:7753792/1‑151 (MQ=255)
                                                                                                        tGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATgg                                    >  1:3718136/1‑151 (MQ=255)
                                                                                                         gggCGGAATAAACAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGc                                   <  1:1978152/151‑1 (MQ=255)
                                                                                                          ggCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGcc                                  >  1:7000491/1‑151 (MQ=255)
                                                                                                           gCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCt                                 <  1:4030105/151‑1 (MQ=255)
                                                                                                             ggAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtc                               >  1:669823/1‑151 (MQ=255)
                                                                                                             cgaagaaAGAAAAGCCCGCGTAAGCCGCCTTTTCTCAGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACGGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtc                               <  2:6047283/150‑1 (MQ=255)
                                                                                                              gAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtct                              <  2:2987626/151‑1 (MQ=255)
                                                                                                              gAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtct                              <  2:639388/151‑1 (MQ=255)
                                                                                                               aaTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtctc                             >  2:7915500/1‑151 (MQ=255)
                                                                                                               aaTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtctc                             >  1:9072025/1‑151 (MQ=255)
                                                                                                                aTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcg                            <  2:7595490/151‑1 (MQ=255)
                                                                                                                 tAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgt                           <  1:11027210/151‑1 (MQ=255)
                                                                                                                 tAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgt                           <  2:8296218/151‑1 (MQ=255)
                                                                                                                   aaGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcg                         >  2:6853931/1‑151 (MQ=255)
                                                                                                                    aGAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgt                        <  2:3452335/151‑1 (MQ=255)
                                                                                                                     gAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtc                       >  2:5324969/1‑151 (MQ=255)
                                                                                                                     gAAAAGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCACAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtc                       <  2:5765763/151‑1 (MQ=255)
                                                                                                                      aaaaGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtcg                      <  2:8077733/151‑1 (MQ=255)
                                                                                                                       aaaGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCgg                     >  1:8357957/1‑151 (MQ=255)
                                                                                                                        aaGGCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCGCAGTAATTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGACGTCTGCCTGTCAGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGga                    >  2:847720/1‑151 (MQ=255)
                                                                                                                           gCCGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGAt                 >  2:1584978/1‑151 (MQ=255)
                                                                                                                           gCCGCGTAAGCGGCCATTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCATGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGAt                 >  2:5638113/1‑151 (MQ=255)
                                                                                                                            ccGCGTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGCAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATa                >  2:7179725/1‑151 (MQ=255)
                                                                                                                             cgcgTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAg               >  1:1408456/1‑151 (MQ=255)
                                                                                                                               cgTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcc             <  2:10284731/150‑1 (MQ=255)
                                                                                                                                gTAAGCGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcct            >  2:4544140/1‑151 (MQ=255)
                                                                                                                                    gcGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCctcctc        >  2:6985685/2‑151 (MQ=255)
                                                                                                                                     cGGCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCtcctcc       >  1:312484/1‑151 (MQ=255)
                                                                                                                                      ggCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCg      >  1:5320155/1‑151 (MQ=255)
                                                                                                                                       gCCTTTTCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCgg     >  2:1958132/1‑151 (MQ=255)
                                                                                                                                         ccttttCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGc    >  1:6174933/3‑151 (MQ=255)
                                                                                                                                          cttttCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACTGCGGAAGTCGGAGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGTCGGTACATGGACTTCTCGTCGTCGGAGCTAGCCTCCTCCGGCt   >  2:9049222/2‑151 (MQ=255)
                                                                                                                                          cttttCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCt   >  2:5146592/2‑151 (MQ=255)
                                                                                                                                          cttttCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCt   >  1:472792/2‑151 (MQ=255)
                                                                                                                                           ttttCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTc  >  1:2563744/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: