Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 888,038 2 bp→TG coding (1562‑1563/1878 nt) PA14_RS04165 → sigma‑54‑dependent Fis family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_008463888,0380CT99.7% 2419.9 / NA 572P521M (CCT→CTG) PA14_RS04165sigma‑54‑dependent Fis family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (322/248);  minor base G (0/2);  total (322/250)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_008463888,0390TG99.5% 2233.6 / NA 572P521M (CCT→CTG) PA14_RS04165sigma‑54‑dependent Fis family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (321/248);  total (322/250)
Fisher's exact test for biased strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 572 total aligned reads displayed.

GAGGATCTCTACTTCCGCCTCGGCGGCGCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCCTGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCAGAGCTGCGTGGCGAAGGCCGGACACCGGCC  >  NC_008463/887893‑888183
                                                                                                                                                 ||                                                                                                                                                
gAGGATCTCTACTTCCGCCTCGGCGGCGCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAAt                                                                                                                                              <  1:2322615/151‑1 (MQ=255)
 aGGATCTCTACTTCCGCCTCGGCGGCGCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATg                                                                                                                                             <  1:6653718/151‑1 (MQ=255)
 aGGATCTCTACTTCCGCCTCGGCGGCGCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATg                                                                                                                                             >  2:917747/1‑151 (MQ=255)
 aGGATCTCTACTTCCGCCTCGGCGGCGCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATg                                                                                                                                             >  2:7462345/1‑151 (MQ=255)
   gATCTCTACTTCCGCCTCGGCGGCGCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGcc                                                                                                                                           >  1:7013666/1‑151 (MQ=255)
    aTCTCTACTTCCGCCTCGGCGGCGCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCt                                                                                                                                          >  2:1923646/1‑151 (MQ=255)
     tctcTACTTCCGCCTCGGCGGCGCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTg                                                                                                                                         <  2:11293104/151‑1 (MQ=255)
      ctctACTTCCGCCTCGGCGGCGCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGc                                                                                                                                        <  2:5672532/151‑1 (MQ=255)
           cTTCCGCCTCGGCGGCGCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGc                                                                                                                                   >  2:500952/1‑151 (MQ=255)
            ttCCGCCTCGGCGGCGCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCt                                                                                                                                  >  2:6951369/1‑151 (MQ=255)
             tCCGCCTCGGCGGCGCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTa                                                                                                                                 <  1:10264210/151‑1 (MQ=255)
             tCCGCCTCGGCGGCGCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTa                                                                                                                                 <  1:7384106/151‑1 (MQ=255)
              ccGCCTCGGCGGCGCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTAt                                                                                                                                >  2:10249072/1‑151 (MQ=255)
                gCCTCGGCGGCGCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCg                                                                                                                              >  1:6799224/1‑151 (MQ=255)
                  cTCGGCGGCGCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCgat                                                                                                                            <  1:7923126/151‑3 (MQ=255)
                    cggcggCGCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTgg                                                                                                                          <  2:9523364/151‑1 (MQ=255)
                      gcggcgCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGcc                                                                                                                        <  1:82167/151‑1 (MQ=255)
                       cggcgCCCGCTTCGAACTGCCGCCGCTACGCGACCGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCg                                                                                                                       <  2:4383846/151‑1 (MQ=255)
                       cggcgCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCg                                                                                                                       <  1:10543381/151‑1 (MQ=255)
                        ggcgCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCgg                                                                                                                      >  2:5120031/1‑151 (MQ=255)
                        ggcgCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCgg                                                                                                                      >  2:5237988/1‑151 (MQ=255)
                        ggcgCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCgg                                                                                                                      >  1:10418844/1‑151 (MQ=255)
                        ggcgCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCgg                                                                                                                      >  1:403639/1‑151 (MQ=255)
                         gcgcCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCggg                                                                                                                     >  1:9874840/1‑151 (MQ=255)
                             ccGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAAt                                                                                                                 <  1:10639987/151‑1 (MQ=255)
                              cGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATg                                                                                                                >  1:4999385/1‑151 (MQ=255)
                                 ttCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTcc                                                                                                             >  1:9746286/1‑151 (MQ=255)
                                   cGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGc                                                                                                           >  1:5276597/1‑151 (MQ=255)
                                     aaCTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCa                                                                                                         >  2:6293872/1‑151 (MQ=255)
                                     aaCTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCa                                                                                                         >  1:9010604/1‑151 (MQ=255)
                                       cTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGt                                                                                                       >  2:3702866/1‑151 (MQ=255)
                                       cTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGt                                                                                                       >  1:10769198/1‑151 (MQ=255)
                                         gccgccGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTg                                                                                                     >  2:8691276/1‑151 (MQ=255)
                                          ccgccgCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGc                                                                                                    <  2:3468836/151‑1 (MQ=255)
                                          ccgccgCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGc                                                                                                    >  2:8915737/1‑151 (MQ=255)
                                          ccgccgCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGc                                                                                                    >  2:10704337/1‑151 (MQ=255)
                                           cgccgcTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCg                                                                                                   >  2:10419815/1‑151 (MQ=255)
                                           cgccgcTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCg                                                                                                   <  2:3576963/151‑1 (MQ=255)
                                            gccgcTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCgg                                                                                                  <  1:2454167/151‑1 (MQ=255)
                                             ccgcTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGc                                                                                                 >  2:6654382/1‑151 (MQ=255)
                                             ccgcTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGc                                                                                                 <  1:9329892/151‑1 (MQ=255)
                                               gcTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCAt                                                                                               >  1:9063659/1‑151 (MQ=255)
                                                cTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCAtg                                                                                              <  1:3331059/151‑1 (MQ=255)
                                                 tACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGACTGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCAtgt                                                                                             >  1:11291004/1‑151 (MQ=255)
                                                 tACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCAtgt                                                                                             >  1:4369883/1‑151 (MQ=255)
                                                  aCGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCAtgtg                                                                                            <  2:7657467/151‑1 (MQ=255)
                                                  aCGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCAtgtg                                                                                            <  1:2689061/151‑1 (MQ=255)
                                                  aCGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCAtgtg                                                                                            <  2:10877030/151‑1 (MQ=255)
                                                      gAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTgc                                                                                        >  1:5111621/1‑151 (MQ=255)
                                                           cAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTAc                                                                                   >  2:5485224/1‑151 (MQ=255)
                                                            aGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACg                                                                                  >  2:9650660/1‑151 (MQ=255)
                                                               gACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCt                                                                               >  2:6802475/1‑151 (MQ=255)
                                                               gACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAAACGCGCACTGCTGCGTACGCATCGAACTCGGAGAGTCAGCCATTGAATGCCTTCTCTGATATCGCTGGCCTGGCAATGTCCGCCAGTTGCGGCATGGGTTGCGCTACGCCt                                                                               >  2:140085/1‑151 (MQ=255)
                                                                aCCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTg                                                                              >  1:2927896/1‑151 (MQ=255)
                                                                  cGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGGCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCAGGTGTTGCGCTACGCCTgcg                                                                            >  1:3872339/1‑151 (MQ=255)
                                                                    gCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTgcgcg                                                                          >  1:8664647/1‑151 (MQ=255)
                                                                      tGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCAGGTGTTGCGCTACGCCTGCGCGCt                                                                        >  1:1954116/1‑151 (MQ=255)
                                                                         cgcTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCtgtg                                                                     <  1:565880/151‑1 (MQ=255)
                                                                          gcTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCCTCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTgc                                                                    >  2:1262247/1‑151 (MQ=255)
                                                                             gATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTgcggc                                                                 >  2:428972/1‑151 (MQ=255)
                                                                                ccGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTgcggcggc                                                              <  2:7078489/151‑1 (MQ=255)
                                                                                 cGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTAAGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTgcggcggcg                                                             <  1:4679920/151‑1 (MQ=255)
                                                                                       aTCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCAccc                                                       >  1:1459125/1‑151 (MQ=255)
                                                                                         ccTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCTGGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTg                                                     >  1:4614164/1‑151 (MQ=255)
                                                                                            cgacgaGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCaa                                                  >  2:8027423/1‑151 (MQ=255)
                                                                                             gacgaGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAAc                                                 >  1:9478155/1‑151 (MQ=255)
                                                                                              acgaGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACt                                                >  2:1598505/1‑151 (MQ=255)
                                                                                                gaGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCg                                              >  1:719710/1‑151 (MQ=255)
                                                                                                gaGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCg                                              >  1:6043123/1‑151 (MQ=255)
                                                                                                  ggAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGcc                                            <  1:10908532/151‑1 (MQ=255)
                                                                                                    aaaCCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGa                                          >  2:9829704/1‑151 (MQ=255)
                                                                                                    aaaCCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGa                                          >  2:3845693/1‑151 (MQ=255)
                                                                                                       ccGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCt                                       >  1:8887252/1‑151 (MQ=255)
                                                                                                        cgcgCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCtt                                      <  1:4880104/151‑1 (MQ=255)
                                                                                                        cgcgCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCtt                                      >  1:4103327/1‑151 (MQ=255)
                                                                                                         gcgcACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTc                                     <  2:2490055/151‑1 (MQ=255)
                                                                                                           gcACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCg                                   >  2:3990120/1‑151 (MQ=255)
                                                                                                            cACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCgg                                  >  2:3135063/1‑151 (MQ=255)
                                                                                                              cTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCa                                <  1:11418617/151‑1 (MQ=255)
                                                                                                               tGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCag                               >  1:2104210/1‑151 (MQ=255)
                                                                                                                 cggcgTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCagag                             <  1:9242768/151‑1 (MQ=255)
                                                                                                                  ggcgTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCAGAGc                            >  2:5751296/1‑151 (MQ=255)
                                                                                                                    cgTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCAGAGCTg                          >  1:6391146/1‑151 (MQ=255)
                                                                                                                     gTACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCAGAGCTGc                         >  2:9109837/1‑151 (MQ=255)
                                                                                                                      tACGCATCGAACTCGGCGAGGCAGCCCTGGAATGCCGGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCAGAGCTGCg                        <  2:1807036/151‑1 (MQ=255)
                                                                                                                         gCATCGAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCAGAGCTGCGTgg                     >  2:2310166/1‑151 (MQ=255)
                                                                                                                              gAACTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCACTGTCCGCCAGTTGCGGCCGGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCCACTCGCCGACCTGCCGGCAGAGCTGCGTGGCGAAg                >  2:10176099/1‑151 (MQ=255)
                                                                                                                               aaCTCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCAGAGCTGCGTGGCGAAgg               >  2:9037941/1‑151 (MQ=255)
                                                                                                                                  tCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCAGAGCTGCGTGGCGAAGGCCg            >  2:8472910/1‑151 (MQ=255)
                                                                                                                                  tCGGCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCGACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCAGAGCTGCGTGGCGAAGGCCg            <  2:4720692/151‑1 (MQ=255)
                                                                                                                                    ggCGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCAGAGCTGCGTGGCGAAGGCCGGa          >  2:7402101/1‑151 (MQ=255)
                                                                                                                                      cGAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCAGAGCTGCGTGGCGAAGGCCGGaca        >  1:1787703/1‑151 (MQ=255)
                                                                                                                                       gAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCAGAGCTGCGTGGCGAAGGCCGGacac       >  1:7671502/1‑151 (MQ=255)
                                                                                                                                       gAGGCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCAGAGCTGCGTGGCGAAGGCCGGacac       >  1:10049094/1‑151 (MQ=255)
                                                                                                                                        agtcAGCACTTGAATGACTTCTCGGCTATCGCTGGCCGGGCATTGTCCGCCTGTTGCTTCTTGTGTTGCTCTACTCCTGCGCTCTGTGCGGCGGCGCCACCCTGCAACTATCAGACCTTACGGCAGAGCTTCGTGGAGAATGACGGacaaa      <  2:2922581/148‑3 (MQ=255)
                                                                                                                                         ggCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCAGAGCTGCGTGGCGAAGGCCGGACACCg     <  1:5432791/151‑1 (MQ=255)
                                                                                                                                         ggCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCGGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCAGAGCTGCGTGGCGAAGGCCGGACACCg     >  1:18794/1‑151 (MQ=255)
                                                                                                                                          gCAGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCAGAGCTGCGTGGCGAAGGCCGGACACCgg    <  2:5950460/151‑1 (MQ=255)
                                                                                                                                           cAGCCCTTGAAGGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGACAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGAGCGGCGGCGCCACCCAGCAACTCGCCGACCTTCCGGCAGAGCTGCGTGGCGAAGGCCGGACACCGGc   <  2:754366/151‑1 (MQ=255)
                                                                                                                                            aGCCCTGGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCAGAGCTGCGTGGCGAAGGCCGGACACCGGcc  >  1:3717270/1‑151 (MQ=255)
                                                                                                                                                 ||                                                                                                                                                
GAGGATCTCTACTTCCGCCTCGGCGGCGCCCGCTTCGAACTGCCGCCGCTACGCGAACGCAGCGACCGGCTGGCGCTGATCCGGCGCATCCTCGACGAGGAAACCGCGCACTGCGGCGTACGCATCGAACTCGGCGAGGCAGCCCCTGAATGCCTGCTCGGCTATCGCTGGCCGGGCAATGTCCGCCAGTTGCGGCATGTGTTGCGCTACGCCTGCGCGCTGTGCGGCGGCGCCACCCTGCAACTCGCCGACCTTCCGGCAGAGCTGCGTGGCGAAGGCCGGACACCGGCC  >  NC_008463/887893‑888183

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