Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084634,347,5800TG28.3% 704.0 / 28.7 348intergenic (‑79/+150)PA14_RS19855/PA14_RS19860hypothetical protein/zonular occludens toxin family protein
Reads supporting (aligned to +/- strand):  ref base T (80/138);  new base G (85/1);  total (208/140)
Fisher's exact test for biased strand distribution p-value = 2.83e-27
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 356 total aligned reads displayed.

GGGGGAGTCCAGCAGTTGTCCGTCGGGTTGCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCGCCCAGCCCATCGCCACCGGAACGAGGAATTCCCAATCGAT  >  NC_008463/4347430‑4347730
                                                                                                                                                      |                                                                                                                                                       
gggggAGTCCAGCAGTTGTCCGTCGGGTTGCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCt                                                                                                                                                         <  1:1646702/151‑1 (MQ=255)
    gAGTCCAGCAGTTGTCCGTCGGGTTGCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCATCgg                                                                                                                                                     >  2:3673625/1‑151 (MQ=255)
      gTCCAGCAGTTGTCCGTCGGGTTGCAACAGGAGGCGGTGTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGGTTATCTGACGGGTTACAGGGGCGTCGGCCGCGCCTTCGGcc                                                                                                                                                   >  2:3403206/1‑151 (MQ=255)
       tCCAGCAGTTGTCCGTCGGGTTGCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCt                                                                                                                                                  <  2:1111741/151‑1 (MQ=255)
              gTTGTCCGTCGGGTTGCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGcgc                                                                                                                                           >  1:7002309/1‑149 (MQ=255)
              gTTGTCCGTCGGGTTGCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCGATCGCTc                                                                                                                                           >  2:8622068/1‑151 (MQ=255)
                 gTCCGTCGGGTTGCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCATCGGCCGATCGCGCATg                                                                                                                                       >  2:9377140/1‑151 (MQ=255)
                  tCCGTCGGGTTGCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCGTCGGCCTAACGCGCATGc                                                                                                                                      >  1:10318342/1‑151 (MQ=255)
                   ccGTCGGGTTGCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCGTCGGCCTATCGCTCATGcc                                                                                                                                     >  2:8411343/1‑151 (MQ=255)
                   ccGTCGGGTTGCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCATCGGGCTATCGCTCAAGcc                                                                                                                                     >  1:438749/1‑151 (MQ=255)
                     gTCGGGTTGCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGGCGCGCCGGCGGCCGATCGCGCAAGCCtg                                                                                                                                   >  1:3799732/1‑150 (MQ=38)
                     gTCGGGTTGCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCtt                                                                                                                                   <  2:4259378/151‑1 (MQ=255)
                       cGGGTTGCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTGGGCCGATCGCTCAGGCCGTgc                                                                                                                                 >  1:2543028/1‑151 (MQ=255)
                          gTTGCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGGCGCGCCGTCGGCCGATAGCACAATCCtttgggt                                                                                                                              >  1:746682/1‑146 (MQ=38)
                          gTTGCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCCCCATCGGCCAATCGCTCATTCCATTCGCt                                                                                                                              >  1:6755143/1‑151 (MQ=38)
                           ttGCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTc                                                                                                                             <  2:2984364/151‑1 (MQ=255)
                            tGCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTcc                                                                                                                            <  1:4165712/151‑1 (MQ=255)
                             gCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTccc                                                                                                                           <  1:3384950/151‑1 (MQ=255)
                              cAACTGGAGGCGGATTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAACCACAGGGAGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCccg                                                                                                                          <  2:2053941/151‑1 (MQ=255)
                              cAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCCCCAAACGCCGATCGCACAATCCATTCGCTcccc                                                                                                                          >  1:2921430/1‑150 (MQ=255)
                               aaCAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTAGGCCTATAGCTCATTCCTTGCGCTacccg                                                                                                                         >  1:8789205/1‑149 (MQ=255)
                                aCAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCccggc                                                                                                                        <  2:9131925/151‑1 (MQ=255)
                                   ggaggCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCccggccgg                                                                                                                     <  1:8731251/151‑1 (MQ=255)
                                   ggaggCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCATGCGCTCccggccgg                                                                                                                     >  2:8727068/1‑151 (MQ=255)
                                     aggCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCcggcg                                                                                                                   <  2:4560554/151‑1 (MQ=255)
                                     aggCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCcggcg                                                                                                                   <  1:7614559/151‑1 (MQ=255)
                                       ggggTTTTTTTGATTCGGTTTCACTTTCTTCTTGGGCTTACCCTGCAATGTCCATTATTGTCCAAAACACAGGGCTTTTATCTGACGTGTTACTGGGGCGTCGGCGGCGCCTTGGGCCTATCGATCTTGCATTGCGCTCCCGGCCggcggc                                                                                                                 <  2:1870273/149‑1 (MQ=255)
                                       gCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCggcggc                                                                                                                 <  2:6800364/151‑1 (MQ=255)
                                        cGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGgcggcg                                                                                                                <  1:4616185/151‑1 (MQ=255)
                                         ggTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGGCGCGCCGTCGGCCGATAGCTCATGGCTTGCGCTCCCGGCCGgcggcgc                                                                                                               >  1:7706028/1‑151 (MQ=255)
                                         ggTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGGCGCGCCGTCGGCCAAAAGCACATTCCGTGCGCGCCCGGCCCgcggcgc                                                                                                               >  1:3136394/1‑151 (MQ=37)
                                          gTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGgcggcgcg                                                                                                              <  1:141194/151‑1 (MQ=255)
                                              tgatgaATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGgccc                                                                                                          <  2:5182608/151‑1 (MQ=255)
                                                atgaATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGGCTTCGGCCGATAGCTCATGCCGTGCGCTCCCGGCCCGCGGCGCGgccccc                                                                                                        >  1:9544242/1‑149 (MQ=37)
                                                   aaTCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATTCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGcccc                                                                                                     >  2:7468290/1‑151 (MQ=255)
                                                     tCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGGCGCGCCATAGGGCGATCGCGCATTCCGTGCGCTCCCGGCCGGCGGGGCGGCCCGccccac                                                                                                   >  1:3892825/1‑149 (MQ=37)
                                                      cGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGGGCCGTAGGCCTATAGCTCAAGCCTGGCGCGACCCGCCGGCGGCGCCGCCCGCCCcaca                                                                                                  >  1:3576086/1‑151 (MQ=25)
                                                       ggTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGGCGCGCCATAGGCCGATACCACAATGCCTTCGCGCCCCGCCCGGGGCGCGGCCCCCCCCCCaa                                                                                                 >  1:49486/1‑150 (MQ=25)
                                                       ggTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCAt                                                                                                 <  1:6943539/151‑1 (MQ=255)
                                                       ggTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCGTCGGCCTATCGCACATGCCTTGCGCACCCGGCCGGCGGCGCGGCCCGCCCCACAt                                                                                                 >  2:4972352/1‑151 (MQ=255)
                                                          ttCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCGTAGGCCTATAGCTCAAGCCATTCGCTCCCCGCCGGCGGCGCGGCCCCCCCcacaatgc                                                                                              >  1:9910792/1‑147 (MQ=25)
                                                           tCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGGCGCGCCGGCGGCCGAACGCGCAAGCCGTGCGCGCCCCGCCCGCGGCGCGGCCCCCCCCCCATggcg                                                                                             >  1:10056851/1‑151 (MQ=25)
                                                           tCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTAAAGCACAAGCCTTGGGCTCCCCGCCGGCGGCGCGGCCCCCCCCGCAAggcg                                                                                             >  1:9271239/1‑151 (MQ=25)
                                                             aCTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCGTCGGCCTATCGCTCAAGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATggcggc                                                                                           >  2:10604837/1‑151 (MQ=255)
                                                              cTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGGCGCGCCGTCGGCCAATAGCGCAAGCCATGCGCGCCCCGCCGGCGGCGCGGCCCCCCCCACAAggcgggg                                                                                          >  1:1831273/1‑149 (MQ=25)
                                                              cTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGgcggcg                                                                                          <  2:3841375/151‑1 (MQ=255)
                                                              cTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCGTCGGCCAATCGCTCATGCCATGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCCCATGgcggcg                                                                                          >  1:7684082/1‑151 (MQ=255)
                                                              cTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGGGTTACAGGGGCGTCGGCCGCGCCGTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCGCATGgcggcg                                                                                          >  2:9543976/1‑151 (MQ=255)
                                                                 tcttcATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGGCCCGCCGTAGGCCAATAGCTCATTCCATTCGCTCCCCGCCCGCCGCGCGGCCCCCCCCACAAGGCGGCGccc                                                                                       >  1:3196441/1‑151 (MQ=25)
                                                                 tcttcATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCGTCGGCCTATCGCTCATGGCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCCTGGCGGCGccc                                                                                       >  2:710998/1‑151 (MQ=255)
                                                                 tcttcATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCGTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCACATGGCGGCGccc                                                                                       >  2:6638928/1‑151 (MQ=255)
                                                                  cttcATGGGCTGACCCAGCAATGTCCATTAATGTCCAAACCACCGGGCGTTTATCTGACGAGTTCCAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCAGGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGcccc                                                                                      <  2:6170635/151‑1 (MQ=25)
                                                                     cATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGGCGCGCCGTCGGCCAATAGCTCAATACATGCGCGCCCCGCCGGCGGCGCGGGCCCCCCCACAAGGCGGGGCCCCAAc                                                                                   >  1:5451624/1‑151 (MQ=21)
                                                                      aTGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTAcc                                                                                  <  1:6844625/151‑1 (MQ=255)
                                                                      aTGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTAcc                                                                                  <  1:6068365/151‑1 (MQ=255)
                                                                           cTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCAATCGCTCATGCCTTTCGCTCCCGGCCGGCGGCGCGGCCCGCCCCACAAGGCGGCGCCCCAACccccgc                                                                             >  2:596537/1‑151 (MQ=25)
                                                                            ttACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCt                                                                            <  1:6403325/151‑1 (MQ=255)
                                                                             tACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCGTCGGCCTATCGCTCCTGCCTTGCGCGCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCAACCGCCGCTa                                                                           >  2:3216488/1‑151 (MQ=255)
                                                                              aCCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCCTCGGCCTAACGCTCAAGCCTAGCGCACCCCGCCGGCGGCGCCGCCCCCCCCACAATGCGGCGCCCCAACCGCCGCTAg                                                                          >  1:5873457/1‑151 (MQ=21)
                                                                               cccAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAgc                                                                         <  1:2815284/151‑1 (MQ=255)
                                                                                 cAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGGCGCGCCGTCGGCCTAACGCGCAAGGCATGCGCGCCCCGCCCGCGGCGCGGCCCCCCCCACAAGGCGGCGCCCCAACCCCCCCGAgccc                                                                       >  1:7290061/1‑149 (MQ=21)
                                                                                    cAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGGCGCGCCGTAGGCCGAAAGCGCATGCCGAGCGCGCCCCGCCCGCGGGGCGGGCCCCCCCGCCAGGCGGCGGCCCAACCCCCGCCAGCGcccg                                                                    >  1:5705884/1‑149 (MQ=16)
                                                                                     aaTGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCGATCGCTCATGCCTTGGGCTCCCGGCCGGCGGCGCGGCCCGCCCCACCTGGCGGCGCCCCGAACGCCGCAAGCGCCGtc                                                                   >  2:4660979/1‑151 (MQ=25)
                                                                                     aaTGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCGTCGGGCGATAGCGCAGGCCGTGCGCGCCCCGCCGGCGGCGCGGCCCCCCCCACATGGGGGCGCCCCGACCCCCCCGCGCGcccgc                                                                   >  1:8348518/1‑148 (MQ=21)
                                                                                      aTGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGtca                                                                  <  2:9098897/151‑1 (MQ=255)
                                                                                        gTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGtcatc                                                                <  2:8900749/151‑1 (MQ=255)
                                                                                         tCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGGCGAACGCTCCTGCCGGGGGCTCCCGGGCGGCGGGGCGGCCCGCCCCTCAATGCGGCGCCCCGCCCGCCGCGAGCGCCGgcctca                                                               >  2:9015626/1‑151 (MQ=21)
                                                                                         tCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCGAGCGCCGtcatcc                                                               >  2:2553202/1‑150 (MQ=255)
                                                                                         tCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCCCCGTAGGCCAATCGCGCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCACATGGGGGCGCCCCAACCGCCCCTAGCGCCGgcaaaa                                                               >  1:5991274/1‑148 (MQ=21)
                                                                                               aaTGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTcc                                                         <  2:4310687/151‑1 (MQ=255)
                                                                                                aTGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCGTCGGCCTATCGCGCATGCCTTGCGCGCCCGGCCGGCGGCGCGGCCCGCCCCGCATGGCGGCGCCCCTACCGCCCCTAGCGCCGAACTAACCGTcca                                                        >  2:6601406/1‑151 (MQ=21)
                                                                                                  gTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGGCTATCGCTCATGCCTTGCGCACCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCAACCGCCGCTAGCGCCGGCATCACCGTccccc                                                      >  2:10084980/1‑151 (MQ=255)
                                                                                                  gTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTccacc                                                      <  1:7925982/151‑1 (MQ=255)
                                                                                                  gTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCATCGGCCAATAGCTCATGCCGTTCGCTCCCGGCCGGCGGCGCGGCCCGCCCCACAATGCGGCGCCCCCAACGCCGCTCTCGCCGTCATCACCCGccacc                                                      >  2:9888706/1‑151 (MQ=18)
                                                                                                   tCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTccacca                                                     <  2:2717591/151‑1 (MQ=255)
                                                                                                   tCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCCCCTTCGGCCAATCGCACATGCCATGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCCGCGCCCCAACCGCCGCTAGCGCCCTCATCACCGTccacca                                                     >  2:5549068/1‑151 (MQ=25)
                                                                                                     cAAAACACCGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGt                                                   <  2:2352545/151‑1 (MQ=255)
                                                                                                     cAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGt                                                   <  2:2834908/151‑1 (MQ=255)
                                                                                                        aaCACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCAATCGCTCATGCCTTGCGCTCCCGGCCCGCGGCGCGGCCCGCCCCTCAAGGCGGCGCCCCAACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGAt                                                >  2:7150371/1‑151 (MQ=255)
                                                                                                           acaGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGACTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCATACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCa                                             <  1:6598024/151‑1 (MQ=255)
                                                                                                           acaGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCGATCGCTCATGCCGTGCGCTCCCCGCCGGCGGCGCGGCCCCCCCCACAAGGCGGCGCCCCAACCCCCGCGAGCGCCGGCCTCACCCTCCACCCCTGATGCa                                             >  2:2967438/1‑151 (MQ=18)
                                                                                                                  gTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCgccca                                      <  1:1452023/151‑1 (MQ=255)
                                                                                                                   tttATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCgcccag                                     <  1:9334503/151‑1 (MQ=255)
                                                                                                                      aTCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCgcccagccc                                  <  1:2262362/151‑1 (MQ=255)
                                                                                                                       tCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCgcccagccca                                 <  2:5755192/151‑1 (MQ=255)
                                                                                                                         tGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCGCCCAGCCCATc                               <  2:7930521/151‑1 (MQ=255)
                                                                                                                         tGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTAACGCGCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCGCCTGGCGGCGCCCCAACCGCCGCTAGCGCCGGCCTCACCCTCCCCCAGGGATGCCGCgcccagcccaag                               >  1:8829842/1‑149 (MQ=21)
                                                                                                                                 tACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCGCCCAGCCCATCGCCACCgg                       <  2:2649065/151‑1 (MQ=255)
                                                                                                                                       ggCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCGCCCAGCCCATCGCCACCGGAACGAg                 <  2:2602964/151‑1 (MQ=255)
                                                                                                                                       ggCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCGCCCAGCCCATCGCCACCGGAACGAg                 <  1:2120574/151‑1 (MQ=255)
                                                                                                                                          gTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCGCCCAGCCCATCGCCACCGGAACGAGGaa              >  2:10353009/1‑151 (MQ=255)
                                                                                                                                            cGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCGCCCAGCCCATCGCCACCGGAACGAGGAAtt            <  1:8679706/151‑1 (MQ=255)
                                                                                                                                                cgcgCCTTCGGCCTCTCGCTCATGCCTAGCGCGCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCCCCGTCCACCAGTGATGCAGCGCCCAGCCCATCGCCACCGGAACGAGGAATTCCCa        <  2:7307678/151‑1 (MQ=255)
                                                                                                                                                cgcgCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCGCCCAGCCCATCGCCACCGGAACGAGGAATTCCCa        <  1:9419353/151‑1 (MQ=255)
                                                                                                                                                 gcgcCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCGCCCAGCCCATCGCCACCGGAACGAGGAATTCCCaa       >  2:1577217/1‑151 (MQ=255)
                                                                                                                                                  cgcCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCGCCCAGCCCATCGCCACCGGAACGAGGAATTCCCAAt      >  1:3300406/1‑151 (MQ=255)
                                                                                                                                                  cgcCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCGCCCAGCCCATCGCCACCGGAACGAGGAATTCCCAAt      <  2:3799732/151‑1 (MQ=255)
                                                                                                                                                     cTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCGCCCAGCCCATCGCCACCGGAACGAGGAATTCCCAATCGa   >  1:72255/1‑151 (MQ=255)
                                                                                                                                                      ttCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCGCCCAGCCCATCGCCACCGGAACGAGGACTTCCCAATCgct  >  2:6300404/1‑149 (MQ=255)
                                                                                                                                                      ttCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCGCCCAGCCCATCGCCACCGGAACGAGGAATTCCCAATCGAt  >  2:7745346/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
GGGGGAGTCCAGCAGTTGTCCGTCGGGTTGCAACAGGAGGCGGTTTTGATGAATCGGTTTCACTTTCTTCATGGGCTTACCCAGCAATGTCCATTAATGTCCAAAACACAGGGCGTTTATCTGACGTGTTACAGGGGCGTCGGCCGCGCCTTCGGCCTATCGCTCATGCCTTGCGCTCCCGGCCGGCGGCGCGGCCCGCCCCTCATGGCGGCGCCCCTACCGCCGCTAGCGCCGTCATCACCGTCCACCAGTGATGCAGCGCCCAGCCCATCGCCACCGGAACGAGGAATTCCCAATCGAT  >  NC_008463/4347430‑4347730

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