Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084636,206,0690TG42.8% 958.1 / inf 760intergenic (‑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 (225/210);  new base G (163/162);  total (388/372)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.89e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 762 total aligned reads displayed.

CGACCGATTCGATTGTCCTGCCTGGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGCCCGGGTCC  >  NC_008463/6205920‑6206218
                                                                                                                                                     |                                                                                                                                                     
  aCCGATTCGATTGTCCTGCCTGGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTgg                                                                                                                                                     >  2:8851788/1‑150 (MQ=34)
          gattGTCCTGCCTGGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTgc                                                                                                                                             <  1:716241/150‑1 (MQ=32)
          gattGTCCTGCCTGGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTgc                                                                                                                                             >  2:5085354/1‑150 (MQ=32)
          gattGTCCTGCCTGGCGCGGGCGAGACAAGGGCGCGCAGGATCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTgc                                                                                                                                             <  2:1747887/150‑1 (MQ=17)
            ttGTCCTGCCTGGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACACAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTgcgg                                                                                                                                           <  1:8173627/150‑1 (MQ=17)
              gTCCTGCCTGGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTgcggcg                                                                                                                                         <  1:246061/150‑1 (MQ=32)
                        gcgcGGGCGCGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACACCACACCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTc                                                                                                                               <  2:1956355/150‑1 (MQ=11)
                           cgGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGc                                                                                                                            <  1:7010744/150‑1 (MQ=17)
                             ggCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCg                                                                                                                          <  1:3545387/150‑1 (MQ=37)
                             ggCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTGTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCg                                                                                                                          <  1:7831222/150‑1 (MQ=17)
                                 agaCAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTCCAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGa                                                                                                                      <  2:5525152/150‑1 (MQ=17)
                                 agaCAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGa                                                                                                                      <  1:8433675/150‑1 (MQ=37)
                                 agaCAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGa                                                                                                                      >  2:10827573/1‑150 (MQ=37)
                                  gaCAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGAc                                                                                                                     <  1:3104976/150‑1 (MQ=17)
                                     aaGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCgg                                                                                                                  <  2:11027904/150‑1 (MQ=37)
                                      aGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCggg                                                                                                                 >  2:8509751/1‑150 (MQ=37)
                                      aGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGACGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCggg                                                                                                                 <  1:10788600/150‑1 (MQ=17)
                                          cgcgcAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAgc                                                                                                             <  1:5406130/150‑1 (MQ=37)
                                             gcAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGcc                                                                                                          >  1:6886217/1‑150 (MQ=37)
                                              cAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCa                                                                                                         >  2:821729/1‑150 (MQ=37)
                                                 gCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAgtc                                                                                                      >  1:4608651/1‑150 (MQ=17)
                                                   tCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAgtcgt                                                                                                    <  2:8050510/150‑1 (MQ=37)
                                                     gATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGtcgtcg                                                                                                  >  1:2716458/1‑150 (MQ=37)
                                                     gATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGtcgtcg                                                                                                  <  1:2272784/150‑1 (MQ=37)
                                                         ggAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGc                                                                                              >  2:1003226/1‑150 (MQ=17)
                                                             gCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGa                                                                                          <  1:11108534/150‑1 (MQ=17)
                                                             gCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGa                                                                                          >  1:7428037/1‑150 (MQ=17)
                                                                gCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGt                                                                                       <  1:47835/150‑1 (MQ=37)
                                                                gCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGt                                                                                       <  1:7553146/150‑1 (MQ=17)
                                                                gCCAACAACAAAGAGGCATCGACCCGGCGCGGGCTTTCCACGCCACGGCGCGGAGGGTCATCCATGGCACGACTGCCAGCGAGGCTGGCCCCGTGCCGCCGTCGAGGCTCAGCCGGCGACCACGGCGCGCAAGTCGACGGAGCCGAAGGt                                                                                       >  2:3975641/1‑150 (MQ=11)
                                                                  cTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGc                                                                                     <  1:6498236/150‑1 (MQ=37)
                                                                   tACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGcc                                                                                    <  2:6008073/150‑1 (MQ=17)
                                                                    acaacaAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCg                                                                                   >  2:9342982/1‑150 (MQ=17)
                                                                      aacaaAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGga                                                                                 >  2:9882006/1‑150 (MQ=37)
                                                                      aacaaAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGga                                                                                 >  1:9160219/1‑150 (MQ=17)
                                                                        caaAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGgaga                                                                               <  1:9759966/150‑1 (MQ=37)
                                                                        caaAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGgaga                                                                               >  2:4959790/1‑150 (MQ=17)
                                                                         aaaGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGTGTCATCGAGGGCACGACGGCCAGCGAGCCTGGCGGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGGCGGAGAc                                                                              <  1:2838406/150‑1 (MQ=255)
                                                                          aaGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACt                                                                             >  2:3110377/1‑150 (MQ=37)
                                                                           aGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACtt                                                                            <  2:9156026/150‑1 (MQ=37)
                                                                           aGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACtt                                                                            >  1:4088464/1‑150 (MQ=17)
                                                                                cATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTgg                                                                       <  2:7376583/150‑1 (MQ=37)
                                                                                cATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTgg                                                                       <  2:5771674/150‑1 (MQ=17)
                                                                                cATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTgg                                                                       >  2:2949976/1‑150 (MQ=17)
                                                                                 aTCGCCCCGGCGCGGGCATTCCTGGCCACGGCGCTGAGGGTCAACGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTggg                                                                      <  2:6334644/150‑1 (MQ=255)
                                                                                 aTCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTggg                                                                      >  2:8742371/1‑150 (MQ=17)
                                                                                  tCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCCGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGt                                                                     >  1:3257780/1‑150 (MQ=37)
                                                                                      cccGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAg                                                                 >  2:3192320/1‑150 (MQ=37)
                                                                                          gcgcGGGCTTTCCTGGCCCCGGCGCTGTGGGACATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGa                                                             <  2:747011/150‑1 (MQ=12)
                                                                                          gcgcGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGa                                                             <  1:1363034/150‑1 (MQ=37)
                                                                                          gcgcGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGa                                                             <  2:7665304/150‑1 (MQ=17)
                                                                                          gcgcGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGa                                                             >  1:115173/1‑150 (MQ=17)
                                                                                           cgcgGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGGGGGTCCAGGTGAt                                                            >  2:1646533/1‑150 (MQ=17)
                                                                                            gcgGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATc                                                           <  2:2852653/150‑1 (MQ=37)
                                                                                            gcgGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATc                                                           <  1:293502/150‑1 (MQ=37)
                                                                                              gggCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCtt                                                         <  1:43563/150‑1 (MQ=17)
                                                                                               ggCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTg                                                        <  1:8300581/150‑1 (MQ=37)
                                                                                                gCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGc                                                       <  1:9393430/150‑1 (MQ=17)
                                                                                                  cttCCTGCCCACGGCGCTGACGCCCAGCGCGGGCACGACGGCCAGCGAGGCTGGCCGCGAGCGGCGGTCGCCGCTCCGCCGGCGCCCGGGGAGCCCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTGCGTGGGTCCAGGTGATCTTGCgg                                                     <  2:6142084/149‑1 (MQ=11)
                                                                                                   ttCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGt                                                    >  1:1556426/1‑150 (MQ=37)
                                                                                                   ttCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGt                                                    >  1:2160200/1‑150 (MQ=17)
                                                                                                    tCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTa                                                   <  2:6888539/150‑1 (MQ=37)
                                                                                                    tCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTa                                                   <  1:101792/150‑1 (MQ=37)
                                                                                                     ccTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTGGCGGTAg                                                  >  2:6748910/1‑150 (MQ=37)
                                                                                                      cTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAgg                                                 <  2:7789488/150‑1 (MQ=37)
                                                                                                      cTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAgg                                                 <  2:4165771/150‑1 (MQ=17)
                                                                                                        ggCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTg                                               >  2:10938061/1‑150 (MQ=17)
                                                                                                           cACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAg                                            <  2:1306274/150‑1 (MQ=17)
                                                                                                            aCGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGt                                           <  2:4379521/150‑1 (MQ=37)
                                                                                                              ggCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGc                                         <  1:8630761/150‑1 (MQ=17)
                                                                                                                cgcTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCAc                                       <  2:2608237/150‑1 (MQ=17)
                                                                                                                cgcTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCAc                                       >  1:201215/1‑150 (MQ=17)
                                                                                                                    gtgGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCt                                   <  1:5970072/150‑1 (MQ=17)
                                                                                                                     tgGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCtt                                  <  2:7252864/150‑1 (MQ=17)
                                                                                                                      gggTCATCGAGGGCACGACGGCCAGCGAGGCTGGGCGGGTGCGGCGGTCGCGGCTCAGCGGGCGACCGGGGAGCGCCAGTCGTCGGAGCGGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGGTCGTCTTc                                 <  2:9208612/150‑1 (MQ=11)
                                                                                                                       ggTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTcc                                <  1:5781106/150‑1 (MQ=37)
                                                                                                                       ggTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTcc                                >  1:1777240/1‑150 (MQ=17)
                                                                                                                          cATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACGTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAgg                             <  2:2188349/150‑1 (MQ=37)
                                                                                                                            tCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCTGGTGATCTTGCGGTAGGTGATGTGCACGTCTGCCAgggg                           >  2:5362826/1‑148 (MQ=255)
                                                                                                                            tCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGAGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAggtg                           >  2:10732894/1‑150 (MQ=37)
                                                                                                                              gAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAggtggg                         >  1:5181027/1‑150 (MQ=17)
                                                                                                                               aGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAggtgggt                        <  2:1469165/150‑1 (MQ=37)
                                                                                                                                 ggCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGGGa                      >  1:9512189/1‑150 (MQ=37)
                                                                                                                                     cgacgGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTg                  >  2:3308571/1‑150 (MQ=37)
                                                                                                                                     cgacgGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTg                  >  1:10389066/1‑150 (MQ=17)
                                                                                                                                     cgacgGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTg                  >  1:629422/1‑150 (MQ=17)
                                                                                                                                       acgGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTggg                >  1:10004277/1‑150 (MQ=17)
                                                                                                                                        cgGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGGGAAGTGGGc               >  2:10533728/1‑150 (MQ=37)
                                                                                                                                          gCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGt             <  1:8994567/150‑1 (MQ=17)
                                                                                                                                           ccAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGtt            >  2:10274019/1‑150 (MQ=37)
                                                                                                                                              gCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGcc         >  2:353943/1‑150 (MQ=37)
                                                                                                                                              gCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGcc         >  2:83155/1‑150 (MQ=37)
                                                                                                                                               cGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGccc        <  1:3693035/150‑1 (MQ=17)
                                                                                                                                               cGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGccc        <  1:10520385/150‑1 (MQ=17)
                                                                                                                                                gAGGCGGGCCGAGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGACGTGCACGTCTGCCAGGTGGGTGACGTGGGCGTTGCCCg       >  2:47779/1‑150 (MQ=11)
                                                                                                                                                 aGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGGGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGCCCgg      >  1:4493654/1‑150 (MQ=37)
                                                                                                                                                  ggCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGCCCggg     >  1:7946578/1‑150 (MQ=37)
                                                                                                                                                   gCTGGCCGCGTGCGGCGGTCGCGGCTCCGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGGGCCGGAGACTTCGTGCGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGCCCGGGt    <  2:9463996/150‑1 (MQ=255)
                                                                                                                                                   gCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGCCCGGGt    <  1:5541476/150‑1 (MQ=37)
                                                                                                                                                     tGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGCCCGGGTcc  <  1:6141900/150‑1 (MQ=35)
                                                                                                                                                     |                                                                                                                                                     
CGACCGATTCGATTGTCCTGCCTGGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGCCCGGGTCC  >  NC_008463/6205920‑6206218

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