Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 5,395,476 (G)5→4 intergenic (‑109/+44) PA14_RS24730 ← / ← PA14_RS24735 hypothetical protein/ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084635,395,4720G.97.3% 1716.1 / 29.1 366intergenic (‑105/+48)PA14_RS24730/PA14_RS24735hypothetical protein/ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/10);  new base . (177/179);  total (177/189)
Fisher's exact test for biased strand distribution p-value = 1.81e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 372 total aligned reads displayed.

GTCAGCCCCAGACAAGGGCGGATGAGGCAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGGCGGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTCCTTGTCGTCCTCTTCCATGCCGAGTTCGGA  >  NC_008463/5395324‑5395618
                                                                                                                                                    |                                                                                                                                                  
           aCAAGGGCGGATGAGGCAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGccc                                                                                                                                     >  2:11054736/1‑151 (MQ=255)
             aaGGGCGGATGAGGCAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCg                                                                                                                                   >  1:6768240/1‑151 (MQ=255)
             aaGGGCGGATGAGGCAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCg                                                                                                                                   <  1:10705676/151‑1 (MQ=255)
              aGGGCGGATGAGGCAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGc                                                                                                                                  >  2:3560163/1‑151 (MQ=255)
               gggCGGATGAGGCAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGcc                                                                                                                                 >  1:7834535/1‑151 (MQ=255)
               gggCGGATGAGGCAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGcc                                                                                                                                 >  1:1831921/1‑151 (MQ=255)
               gggCGGATGAGGCAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGGTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGcc                                                                                                                                 >  1:1729965/1‑151 (MQ=255)
                  cGGATGAGGCAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTc                                                                                                                              >  1:7264279/1‑151 (MQ=255)
                   ggATGAGGCAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCg                                                                                                                             <  1:6645722/151‑1 (MQ=255)
                     aTGAGGCAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGc                                                                                                                           <  1:8523041/151‑1 (MQ=255)
                      tGAGGCAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCa                                                                                                                          <  2:581775/151‑1 (MQ=255)
                      tGAGGCAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCa                                                                                                                          <  2:1999062/151‑1 (MQ=255)
                       gAGGCAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAg                                                                                                                         >  2:9752235/1‑151 (MQ=255)
                       gAGGCAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAg                                                                                                                         >  2:5164241/1‑151 (MQ=255)
                       gAGGCAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAg                                                                                                                         >  1:2491619/1‑151 (MQ=255)
                          gCAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGtt                                                                                                                      >  2:8437031/1‑151 (MQ=255)
                           cAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTc                                                                                                                     >  1:3039006/1‑151 (MQ=255)
                            agtggTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTcc                                                                                                                    <  2:6522099/148‑1 (MQ=255)
                            agagGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTcc                                                                                                                    <  2:5384473/151‑1 (MQ=255)
                            agagGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTcc                                                                                                                    >  1:3157821/1‑151 (MQ=255)
                             gagGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCg                                                                                                                   <  2:6909784/151‑1 (MQ=255)
                              agGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCgg                                                                                                                  >  2:10133630/1‑151 (MQ=255)
                              agGTTACGCGGCATGGCGGCAACCTCACGGGGGTTCCCAGTTTTGAGCCTCCACCTGCCGCCAGGTTCCGCATCGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCgg                                                                                                                  <  2:9309058/151‑1 (MQ=255)
                                gTTACGCGGCATGGCGGCAACCTCCCGGGGGTTCCCTGTTTGGAGCCTCCACCTGCCGCCTGGGTCCGCATCGCGCCGGGGGGAACTTGGCCACCGGGCTGAAACGCACGAGGG‑‑GGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCa                                                                                                                <  1:3259780/151‑1 (MQ=255)
                                  tACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAAc                                                                                                              >  1:6175065/1‑151 (MQ=255)
                                   aCGCGGCAGGGCGGCAGCCTCCCGGGGGTTCCCAGTTTTGCGCCTACCCCAGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACCGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAAcg                                                                                                             <  2:11371395/151‑1 (MQ=255)
                                       ggCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGcc                                                                                                         >  1:4753724/1‑151 (MQ=255)
                                           tGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGAc                                                                                                     >  2:1283874/1‑151 (MQ=255)
                                            ggcggcATCCTCCCGGGGGCGCCCTGTTTTGAGCCACCCCCGGCCGCCGGGTGCCGAATCGCGCCGGGGGGAACTGCACCACCGGGCTGAACCGCCCGTGGGCGGGGTTTCCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACt                                                                                                    <  1:3847672/151‑1 (MQ=255)
                                            ggcggcATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACt                                                                                                    >  1:10181993/1‑151 (MQ=255)
                                             gcggcATCCTCCCGGGGGTTCCCAGTTTTGAGCCTCCCCCTGCCGCCCGGTACCGCATACCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGTGGGCGGGGCTTTCCCGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTc                                                                                                   <  1:9995411/151‑1 (MQ=255)
                                               ggcATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTACAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAg                                                                                                 >  2:11145/1‑151 (MQ=255)
                                                gcATCCTCCCGGGGGTTCCCAGTTTTGCGCCTCCACCTGCCGCCAGGTTCCGCATCGCGCCGGCGGGAACTTGACCACAGGGCTGAAACGCACGCGGGCGGGGCTTACCCGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGt                                                                                                <  1:4467994/151‑1 (MQ=255)
                                                   tCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCct                                                                                             >  2:9611243/1‑151 (MQ=255)
                                                   tCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACGCAGTCct                                                                                             >  2:4503700/1‑151 (MQ=255)
                                                     cTCAGGGGGGTTCCCTGTTTTGAGCATACACCTGCCGCCTGGTTCCGCATAGCGCCGGCGGTAACTTGTCCACAGGGCTGAAACGCACGTGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCcttc                                                                                           <  1:1481189/151‑1 (MQ=255)
                                                     cTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCcttc                                                                                           >  2:10309404/1‑151 (MQ=255)
                                                     cTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTGGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCcttc                                                                                           <  2:8692177/151‑1 (MQ=255)
                                                       cACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCcttctt                                                                                         >  2:923797/1‑151 (MQ=255)
                                                       cACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCcttctt                                                                                         >  1:8902333/1‑151 (MQ=255)
                                                        aCGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCTGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTCCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCttcttc                                                                                        <  2:2891999/151‑1 (MQ=255)
                                                          gggggTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTcgc                                                                                      <  2:2097938/151‑1 (MQ=255)
                                                            gggTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGAATAGCGCCGGCGGGAACATGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTACCGTCAACGCGCCGGACTCAGTCCTTCTTcgccg                                                                                    >  2:6396889/1‑151 (MQ=255)
                                                              ggtCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCcgcg                                                                                  >  1:229932/3‑151 (MQ=255)
                                                               ttCCCTTTTTTGAGCCTACACCTGCCGAAAGGTTCCTCATAGCGCCGGCGGGAATTTGTCAACAGGGCTGAAACGCACTTGGGCGGGGCTTTCCCGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCgg                                                                                 <  1:6952374/151‑1 (MQ=255)
                                                                tCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGc                                                                                >  2:9487366/1‑151 (MQ=255)
                                                                 cccAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCTTAGGGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGcc                                                                               <  2:3032317/151‑1 (MQ=255)
                                                                  ccAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCt                                                                              <  2:4368492/151‑1 (MQ=255)
                                                                         tGAGCCAACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCtt                                                                       <  1:921607/151‑1 (MQ=255)
                                                                            gCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTcgc                                                                    >  2:3772510/1‑151 (MQ=255)
                                                                             ccTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTcgcc                                                                   >  1:10811969/1‑151 (MQ=255)
                                                                                acacCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCcgcg                                                                >  1:9488493/1‑151 (MQ=255)
                                                                                  cccTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTTAAACGCACTAGGGCGGGGCTTACCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGAc                                                              <  1:8191107/150‑1 (MQ=255)
                                                                                  acCTGCCGCCTGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGAc                                                              <  2:6752827/151‑1 (MQ=255)
                                                                                  acCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGAc                                                              >  2:7916443/1‑151 (MQ=255)
                                                                                   ccTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACg                                                             <  1:7522680/151‑1 (MQ=255)
                                                                                     tGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGt                                                           >  2:7255201/1‑151 (MQ=255)
                                                                                       ccgccAGGTTCCGCATAGCGCCGGCGGGAACGTGTCCACAGGGCTGAAACGCCCGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGtct                                                         <  2:6155800/151‑1 (MQ=255)
                                                                                          ccAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACCGGGCTGAAACGAACGTGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCtt                                                      <  1:5641023/151‑1 (MQ=255)
                                                                                           cAGGTTCCGCATCGCGCCGGCGGGAACTTGTCCACAGGGCTGAACCGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCttt                                                     <  1:3214033/151‑1 (MQ=255)
                                                                                               ttCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTcc                                                 <  2:1234387/151‑1 (MQ=255)
                                                                                                tCCTCATCGCGCCGCCGGGAACTTGTCCACCGGGCTGAACCGCAGGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCTGGTTCCGGCACCGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTccg                                                <  1:6089510/151‑1 (MQ=255)
                                                                                                tCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTccg                                                <  1:2736534/151‑1 (MQ=255)
                                                                                                tCCGCATAGCGCCGGCGGGAACTTGACCACAGGGCTGAAACGCACGAGGGCGGGGATTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTccg                                                <  1:4146266/151‑1 (MQ=255)
                                                                                                 ccGCCTACCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTccgc                                               <  2:7379885/151‑1 (MQ=255)
                                                                                                  cGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTccgcc                                              >  1:11204974/1‑151 (MQ=255)
                                                                                                  cGCATAGCGCCGGCGGGAACTGGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTccgcc                                              <  2:10641372/151‑1 (MQ=255)
                                                                                                    cATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCgg                                            >  1:5606761/1‑151 (MQ=255)
                                                                                                     aTAGCGCCGGCGGGACCTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGc                                           >  2:4812047/1‑151 (MQ=255)
                                                                                                     aTAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGCCGGTCTCTTTCTCCGCCGGc                                           <  2:7007332/151‑1 (MQ=255)
                                                                                                     aTAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGc                                           >  1:636163/1‑151 (MQ=255)
                                                                                                            cggcggGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCt                                    >  2:3803257/1‑151 (MQ=255)
                                                                                                             ggcggGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCtt                                   >  1:10016948/1‑151 (MQ=255)
                                                                                                             ggcggGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCtt                                   >  1:10279987/1‑151 (MQ=255)
                                                                                                              gcggGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCttt                                  >  1:7785696/1‑151 (MQ=255)
                                                                                                              gcggGAACTTGGCCACAGGGCTGAAACGCACGAGGGCGGGGCTGACCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCttt                                  <  2:5982307/151‑1 (MQ=255)
                                                                                                               cggGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCtttt                                 <  1:7401184/151‑1 (MQ=255)
                                                                                                               cggGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCtttt                                 <  2:1165275/151‑1 (MQ=255)
                                                                                                               cggGAACTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTGTCTCCGCCGGCTTGTCCtttt                                 >  1:2365473/1‑151 (MQ=255)
                                                                                                                 ggAACTTGTCCACCGGGCTGAAACGCACGTGGG‑‑GGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTcc                               <  1:5246531/151‑1 (MQ=255)
                                                                                                                  gAACTTGGCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTCCt                              <  2:552325/151‑1 (MQ=255)
                                                                                                                     cTTGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTCCTTgt                           <  2:1371864/151‑1 (MQ=255)
                                                                                                                      ttGTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTCCTTgtc                          <  2:2516735/151‑1 (MQ=255)
                                                                                                                        gTCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTCCTTgtcgt                        <  2:7532767/151‑1 (MQ=255)
                                                                                                                         tCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTCCTTgtcgtc                       <  2:5441178/151‑1 (MQ=255)
                                                                                                                         tCCACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTGTTCCTTgtcgtc                       >  1:10592708/1‑151 (MQ=255)
                                                                                                                          ccACAGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTCCTTGTCGTcc                      <  2:3541701/151‑1 (MQ=255)
                                                                                                                           cacaGGGCTGAAACGCCCGAGGGCGGGGCTGTCCAGCCCCCCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTCCTTGTCGTCct                     <  1:4228389/151‑1 (MQ=255)
                                                                                                                           cacaGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTCCTTGTCGTCct                     <  2:7800958/151‑1 (MQ=255)
                                                                                                                             caGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTCCTTGTCGTCctct                   <  2:5350875/151‑1 (MQ=255)
                                                                                                                              aGGGCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTCCTTGTCGTCCTCtt                  <  2:5081447/151‑1 (MQ=255)
                                                                                                                                ggCTGAAACGCACGAGGGCGCGGCTTACCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTAATCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTACGCCGGCTTGTCCTTTTACTTGTCGTCCTCTTcc                >  2:4057668/1‑151 (MQ=255)
                                                                                                                                 gCTGAAACGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTCCTTGTCGTCCTCTTCCa               <  2:4602449/151‑1 (MQ=255)
                                                                                                                                  cTGAAACGCACGAGGGCGGGGCGTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTCCTTGTCGTCCTCTTCCAt              <  2:8938136/151‑1 (MQ=255)
                                                                                                                                     aaaCGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTCCTTGTCGTCCTCTTCCATGcc           >  2:4255619/1‑151 (MQ=255)
                                                                                                                                      aaCGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTCCTTGTCGTCCTCTTCCATGCCg          >  1:2999767/1‑151 (MQ=255)
                                                                                                                                       aCGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTCCTTGTCGTCCTCTTCCATGCCGa         <  2:4684402/151‑1 (MQ=255)
                                                                                                                                       aCGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTATCCGCCGGCTTGTCCTTTTCCTTGTCGTCCTCTTCCATGCCGa         >  1:8400020/1‑151 (MQ=255)
                                                                                                                                        cGCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTACCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACTGTCTCTTTCTCAGCCGTCTTGTCCATTTCCTTGTCGTCCTCTTCCATTCCGAg        <  2:1753297/151‑1 (MQ=255)
                                                                                                                                         gCACGAGGGCGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTCCTTGTCGTCCTCTTCCATGCCGAGt       >  2:10488730/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTCAGCCCCAGACAAGGGCGGATGAGGCAGAGGTTACGCGGCATGGCGGCATCCTCACGGGGGTTCCCAGTTTTGAGCCTACACCTGCCGCCAGGTTCCGCATAGCGCCGGCGGGAACTTGTCCACAGGGCTGAAACGCACGAGGGGCGGGGGCTTTCCAGCCCCGCCCTCGGCAGGTTCCGGCAACGCGCCGGACTCAGTCCTTCTTCGCCGCGGCCTTGTCCTTCGCCGCGACGGTCTCTTTCTCCGCCGGCTTGTCCTTTTCCTTGTCGTCCTCTTCCATGCCGAGTTCGGA  >  NC_008463/5395324‑5395618

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: