Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084636,206,0690TG48.6% 664.1 / inf 745intergenic (‑174/+23)PA14_RS28345/PA14_RS28350type VI secretion system tip protein VgrG/Hcp family type VI secretion system effector
Reads supporting (aligned to +/- strand):  ref base T (186/197);  new base G (181/181);  total (367/378)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 750 total aligned reads displayed.

GACCGATTCGATTGTCCTGCCTGGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGCCCGGGTCC  >  NC_008463/6205921‑6206218
                                                                                                                                                    |                                                                                                                                                     
gACCGATTCGATTGTCCTGCCTGGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTg                                                                                                                                                      >  1:8379122/1‑150 (MQ=33)
 aCCGATTCGATTGTCCTGCCTGGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTgg                                                                                                                                                     >  1:843231/1‑150 (MQ=34)
  ccGATTCGATTGTCCTGCCTGGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGc                                                                                                                                                    >  1:1833086/1‑150 (MQ=34)
    gattcgcttGTCCTGCCTGGCGCGGGCGAGTCAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACTACAAAGCGGCATCGACCCGGCGCGGGCGTTCCTGGGCTCGGGGCGGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCcg                                                                                                                                                  >  2:4501989/1‑150 (MQ=11)
       tcgattGTCCTGCCTGGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGt                                                                                                                                               >  1:5781896/1‑150 (MQ=32)
         gattGTCCTGCCTGGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTgc                                                                                                                                             <  1:6838885/150‑1 (MQ=32)
               cctgcctgGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAACGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGt                                                                                                                                       <  2:10626441/150‑1 (MQ=17)
               cctgcctgGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGt                                                                                                                                       <  2:5577288/150‑1 (MQ=32)
                     tgGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGc                                                                                                                                 >  1:10506494/1‑150 (MQ=35)
                       gcgcGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTc                                                                                                                               <  1:4642937/150‑1 (MQ=37)
                       gcgcGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTc                                                                                                                               >  1:2726795/1‑150 (MQ=17)
                          cgGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGc                                                                                                                            <  2:289581/150‑1 (MQ=37)
                          cgGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGc                                                                                                                            >  1:6142671/1‑150 (MQ=17)
                           gggCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCTGcc                                                                                                                           >  1:7340083/1‑150 (MQ=37)
                           gggCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGcc                                                                                                                           >  2:5783187/1‑150 (MQ=17)
                            ggCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCg                                                                                                                          >  1:3083863/1‑150 (MQ=17)
                             gCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCgg                                                                                                                         >  1:5410646/1‑150 (MQ=37)
                             gCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCgg                                                                                                                         >  1:8146497/1‑150 (MQ=17)
                              cGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGc                                                                                                                        <  2:6272218/150‑1 (MQ=37)
                                agaCAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGa                                                                                                                      >  1:2593056/1‑150 (MQ=17)
                                agaCAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGa                                                                                                                      >  1:9022948/1‑150 (MQ=17)
                                  aCAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGAcc                                                                                                                    >  1:405728/1‑150 (MQ=37)
                                     aGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCggg                                                                                                                 <  1:10255147/150‑1 (MQ=17)
                                          gcgcAGGCTCGATCGGAAGCTGCCTACAACACAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAgcg                                                                                                            >  2:2156979/1‑150 (MQ=17)
                                             cAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCa                                                                                                         >  1:709559/1‑150 (MQ=17)
                                             cAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCa                                                                                                         >  2:859053/1‑150 (MQ=17)
                                              aGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAg                                                                                                        <  1:2359880/150‑1 (MQ=17)
                                                    gATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGtcgtcg                                                                                                  <  2:6403532/150‑1 (MQ=17)
                                                      tCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGa                                                                                                >  1:10350508/1‑150 (MQ=37)
                                                              tGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAgg                                                                                        >  1:8102589/1‑150 (MQ=17)
                                                               gCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGt                                                                                       >  1:4241346/1‑150 (MQ=37)
                                                                 cTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGc                                                                                     <  1:6461932/150‑1 (MQ=17)
                                                                 cTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGc                                                                                     >  2:10337536/1‑150 (MQ=17)
                                                                     aacaaAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGga                                                                                 >  1:8243564/1‑150 (MQ=17)
                                                                      acaaAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGgag                                                                                <  2:9583509/150‑1 (MQ=17)
                                                                         aaGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACt                                                                             >  2:5470916/1‑150 (MQ=37)
                                                                          aGCGGCATCGACCCGGCGCGGGCTTGCCTGGCCACGGCGCTGTGGGTCAGCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACtt                                                                            <  2:8453247/150‑1 (MQ=255)
                                                                           gCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTc                                                                           >  1:5565359/1‑150 (MQ=37)
                                                                            cGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCg                                                                          <  1:3406100/150‑1 (MQ=37)
                                                                              gCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTg                                                                        >  1:5952980/1‑150 (MQ=17)
                                                                                aTCGCCCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTggg                                                                      <  2:6886559/150‑1 (MQ=37)
                                                                                aTCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTGATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTggg                                                                      >  2:930931/1‑150 (MQ=37)
                                                                                aTCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGTGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTggg                                                                      <  1:8876609/150‑1 (MQ=37)
                                                                                  cGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTc                                                                    <  1:10104870/150‑1 (MQ=17)
                                                                                   gACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTcc                                                                   <  1:3640725/150‑1 (MQ=37)
                                                                                   gACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTcc                                                                   >  1:3864612/1‑150 (MQ=37)
                                                                                    aCCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCa                                                                  >  1:10572445/1‑150 (MQ=37)
                                                                                    aCCCGGCGCGGGCTTTCCTGGCCACGGCGCGGTGGGGCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGGCGGCGCCGGAGGTGCCGGAGACTTCGTGGGTCCa                                                                  >  2:9713717/1‑150 (MQ=11)
                                                                                          cgcgGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGAt                                                            >  1:6649631/1‑150 (MQ=37)
                                                                                          cgcgGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGAt                                                            >  1:10014652/1‑150 (MQ=17)
                                                                                           gcgGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATc                                                           <  2:8086612/150‑1 (MQ=37)
                                                                                           gcgGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATc                                                           <  2:763267/150‑1 (MQ=17)
                                                                                           gcgGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATc                                                           >  1:9392819/1‑150 (MQ=17)
                                                                                            cgGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACGTCGTGGGTCCAGGTGATCt                                                          >  2:4278354/1‑150 (MQ=37)
                                                                                               gCTTTCCTGGCCTCGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGc                                                       <  2:2668546/150‑1 (MQ=17)
                                                                                               gCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGc                                                       <  1:8765779/150‑1 (MQ=17)
                                                                                                  ttCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGt                                                    >  2:4093739/1‑150 (MQ=37)
                                                                                                   tCCTGGCCACGGCGCTGTGGGTGTTCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCGCAGCCGGCGACCGGGGAGCGCCTGTCGTCGGCGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGAGCGTGCGGTa                                                   >  2:6621471/1‑150 (MQ=11)
                                                                                                    ccTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAg                                                  <  2:7084947/150‑1 (MQ=37)
                                                                                                     cTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAgg                                                 >  2:5238344/1‑150 (MQ=17)
                                                                                                     cTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAgg                                                 <  1:2807717/150‑1 (MQ=17)
                                                                                                      tGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGt                                                >  1:9517055/1‑150 (MQ=37)
                                                                                                       ggCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTg                                               <  2:1948966/150‑1 (MQ=37)
                                                                                                       ggCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTg                                               <  1:4686976/150‑1 (MQ=37)
                                                                                                        gCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGa                                              <  2:1564965/150‑1 (MQ=37)
                                                                                                          cACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAg                                            <  1:2305435/150‑1 (MQ=17)
                                                                                                            cGGCGCTGTGGTTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTg                                          <  2:7577487/150‑1 (MQ=37)
                                                                                                             ggCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGc                                         <  1:3275215/150‑1 (MQ=17)
                                                                                                              gcgcTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCa                                        <  1:8018382/150‑1 (MQ=17)
                                                                                                                 cTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGt                                     <  2:5079287/150‑1 (MQ=37)
                                                                                                                  tgtgGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTc                                    >  2:7922290/1‑150 (MQ=37)
                                                                                                                  tgtgGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTc                                    <  2:9115894/150‑1 (MQ=17)
                                                                                                                  tgtgGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTc                                    <  2:10422833/150‑1 (MQ=17)
                                                                                                                   gtgGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCt                                   <  1:5043998/150‑1 (MQ=17)
                                                                                                                     gggTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTc                                 <  2:869415/150‑1 (MQ=37)
                                                                                                                     gggTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTc                                 >  2:4848810/1‑150 (MQ=17)
                                                                                                                      ggTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTcc                                <  1:3748159/150‑1 (MQ=37)
                                                                                                                       gTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCa                               >  2:3592429/1‑150 (MQ=37)
                                                                                                                        tCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAg                              <  1:9917632/150‑1 (MQ=37)
                                                                                                                             gAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAggtggg                         >  1:8624308/1‑150 (MQ=17)
                                                                                                                             gAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAggtggg                         <  2:10129982/150‑1 (MQ=17)
                                                                                                                              aGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAggtgggt                        >  1:7630152/1‑150 (MQ=17)
                                                                                                                               gggCACGACGGCCAGCGAGGCTGGCCTCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGATGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCAAggtgggtg                       <  2:5996335/150‑1 (MQ=255)
                                                                                                                                ggCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGa                      <  1:7168445/150‑1 (MQ=37)
                                                                                                                                 gCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGaa                     <  1:1534633/150‑1 (MQ=37)
                                                                                                                                    cgacgGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCCGTAGGTGAAGTGCACGTCCTCCAGGTGGGTGAAGTg                  <  2:3344464/150‑1 (MQ=12)
                                                                                                                                      acgGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTggg                <  2:7276270/150‑1 (MQ=37)
                                                                                                                                      acgGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTggg                <  2:2244686/150‑1 (MQ=17)
                                                                                                                                          ccAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGtt            >  1:5204901/1‑150 (MQ=37)
                                                                                                                                          ccAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGtt            >  1:1198258/1‑150 (MQ=17)
                                                                                                                                           ccgCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGCCCGGGGAGCGCCAGTCGACGGCGCCGGAGGGGCCGGAGACTTCGAGGGTCCAGGTGATCTCGCCGTAGGAGACGTGCACGTCTTCCAGGTGGGAGAAGTGGGCGTTg           <  2:4757352/148‑1 (MQ=255)
                                                                                                                                           cAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTg           >  2:2363739/1‑150 (MQ=37)
                                                                                                                                            aGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGTTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGc          <  1:919865/150‑1 (MQ=37)
                                                                                                                                            aGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGc          <  2:1138603/150‑1 (MQ=37)
                                                                                                                                            aGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGc          <  1:9725313/150‑1 (MQ=37)
                                                                                                                                            aGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGGGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGc          >  1:1408875/1‑150 (MQ=17)
                                                                                                                                             gCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGcc         <  2:82614/150‑1 (MQ=37)
                                                                                                                                             gCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGcc         >  1:452400/1‑150 (MQ=17)
                                                                                                                                                aGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGGAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGCCCgg      >  1:7735114/1‑150 (MQ=17)
                                                                                                                                                    tGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGCCCGGGTcc  >  2:3157970/1‑150 (MQ=35)
                                                                                                                                                    |                                                                                                                                                     
GACCGATTCGATTGTCCTGCCTGGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGCCCGGGTCC  >  NC_008463/6205921‑6206218

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: