Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084636,206,0690TG48.7% 618.0 / inf 692intergenic (‑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 (184/171);  new base G (161/176);  total (345/347)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.01e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 696 total aligned reads displayed.

GCGACCGATTCGATTGTCCTGCCTGGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGCCCGGGTCCT  >  NC_008463/6205919‑6206219
                                                                                                                                                      |                                                                                                                                                      
         tcgattGTCCTGCCTGGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCGGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTg                                                                                                                                               >  2:2738039/1‑151 (MQ=14)
                 cctgcctgGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTc                                                                                                                                       <  1:8203332/151‑1 (MQ=31)
                 cctgcctgGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTc                                                                                                                                       <  1:1704386/151‑1 (MQ=31)
                   tgcctgGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTcgc                                                                                                                                     >  1:1102304/1‑151 (MQ=32)
                         gcgcGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCCTCGACCCGGCGCGGGCTTTCCGGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCa                                                                                                                               <  2:3761278/151‑1 (MQ=9)
                               gCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGc                                                                                                                         <  2:8223645/151‑1 (MQ=14)
                                cGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCg                                                                                                                        <  2:2397431/151‑1 (MQ=35)
                                 gagaCAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGa                                                                                                                       >  2:5047921/1‑151 (MQ=14)
                                  agaCAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGAc                                                                                                                      >  1:9231555/1‑151 (MQ=35)
                                  agaCAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGAc                                                                                                                      >  1:2291489/1‑151 (MQ=35)
                                   gaCAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGAcc                                                                                                                     >  2:5793671/1‑151 (MQ=14)
                                    aCAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCg                                                                                                                    >  1:3349185/1‑151 (MQ=14)
                                      aaGGGCGCGCAGGCTCGATCGGAAGCTGCCTACACCAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCggg                                                                                                                  <  2:3179512/151‑1 (MQ=25)
                                        gggCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGa                                                                                                                >  2:2319890/1‑151 (MQ=35)
                                        gggCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGa                                                                                                                <  2:1409075/151‑1 (MQ=14)
                                         ggCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAg                                                                                                               >  2:5982744/1‑151 (MQ=14)
                                           cgcgcAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAgcg                                                                                                             >  1:9101918/1‑151 (MQ=35)
                                           cgcgcAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAgcg                                                                                                             <  1:9901625/151‑1 (MQ=14)
                                           cgcgcAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAgcg                                                                                                             <  2:9743877/151‑1 (MQ=14)
                                            gcgcAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAgcgc                                                                                                            <  2:5327435/151‑1 (MQ=14)
                                               cAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAg                                                                                                         >  1:8748917/1‑151 (MQ=14)
                                                 ggCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCGGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAgtc                                                                                                       >  2:5473164/1‑151 (MQ=25)
                                                   cTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAgtcgt                                                                                                     <  1:7550585/151‑1 (MQ=14)
                                                     cGATCGGAAGCTGCCTACAACAACGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGtcgtcg                                                                                                   <  1:6462161/151‑1 (MQ=25)
                                                      gATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCgg                                                                                                  >  1:10012456/1‑151 (MQ=35)
                                                       aTCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGa                                                                                                 >  2:923754/1‑151 (MQ=35)
                                                       aTCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGa                                                                                                 >  1:10508929/1‑151 (MQ=14)
                                                        tCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAg                                                                                                >  1:9421266/1‑151 (MQ=35)
                                                         cGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGc                                                                                               >  1:9601524/1‑151 (MQ=14)
                                                          ggAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGcc                                                                                              >  2:8777530/1‑151 (MQ=35)
                                                          ggAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGcc                                                                                              <  2:10315064/151‑1 (MQ=35)
                                                          ggAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGcc                                                                                              <  1:6706611/151‑1 (MQ=35)
                                                          ggAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGcc                                                                                              <  1:4146306/151‑1 (MQ=14)
                                                          ggAAGCTGCCAACAACTAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGcc                                                                                              <  2:6783129/151‑1 (MQ=25)
                                                           gAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCg                                                                                             >  1:10252321/1‑151 (MQ=35)
                                                           gAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCg                                                                                             <  2:1883259/151‑1 (MQ=14)
                                                           gAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCg                                                                                             <  2:749667/151‑1 (MQ=14)
                                                             aGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGa                                                                                           <  2:4342213/151‑1 (MQ=35)
                                                             aGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGa                                                                                           >  1:5708639/1‑151 (MQ=35)
                                                               cTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAgg                                                                                         >  2:1137095/1‑151 (MQ=14)
                                                                 gCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTg                                                                                       >  1:1477513/1‑151 (MQ=35)
                                                                 gCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTg                                                                                       <  1:2960717/151‑1 (MQ=35)
                                                                  ccTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGc                                                                                      >  1:7157191/1‑151 (MQ=14)
                                                                   cTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGTCCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGcc                                                                                     >  1:899225/1‑151 (MQ=25)
                                                                    tACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCg                                                                                    >  1:9772459/1‑151 (MQ=14)
                                                                    tACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGCCCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCg                                                                                    >  2:9379312/1‑151 (MQ=9)
                                                                      caacaaAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGga                                                                                  >  1:5531822/1‑151 (MQ=35)
                                                                       aacaaAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGgag                                                                                 <  1:1852752/151‑1 (MQ=14)
                                                                       aacaaAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGgag                                                                                 <  1:2714954/151‑1 (MQ=14)
                                                                          aaaGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACt                                                                              >  2:7821600/1‑151 (MQ=35)
                                                                           aaGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACtt                                                                             <  1:2422138/151‑1 (MQ=14)
                                                                            aGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTc                                                                            >  1:1265260/1‑151 (MQ=35)
                                                                                gCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTgg                                                                        <  2:7700572/151‑1 (MQ=14)
                                                                                 cATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTggg                                                                       >  1:6255890/1‑151 (MQ=35)
                                                                                   tCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTc                                                                     <  1:6619850/151‑1 (MQ=35)
                                                                                     gACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCTTCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCATCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCa                                                                   >  2:624248/1‑151 (MQ=25)
                                                                                     gACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCa                                                                   >  1:733321/1‑151 (MQ=14)
                                                                                     gACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCa                                                                   >  1:7453553/1‑151 (MQ=14)
                                                                                       cccGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCCGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAgg                                                                 <  2:4504798/151‑1 (MQ=14)
                                                                                          ggCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGa                                                              >  1:9998059/1‑151 (MQ=35)
                                                                                          ggCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGa                                                              >  1:3157865/1‑151 (MQ=14)
                                                                                               gggCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTg                                                         >  1:6169878/1‑151 (MQ=35)
                                                                                               gggCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTg                                                         >  1:4919311/1‑151 (MQ=35)
                                                                                                ggCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGc                                                        <  1:7005489/151‑1 (MQ=35)
                                                                                                ggCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGc                                                        <  2:4156338/151‑1 (MQ=14)
                                                                                                 gCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCGTGCg                                                       >  2:2616762/1‑151 (MQ=14)
                                                                                                    ttCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTa                                                    <  1:504821/151‑1 (MQ=14)
                                                                                                      ccTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAgg                                                  <  2:4924582/151‑1 (MQ=14)
                                                                                                         ggCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCTGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGa                                               <  2:10745823/151‑1 (MQ=25)
                                                                                                         ggCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGa                                               <  2:4772781/151‑1 (MQ=14)
                                                                                                            cACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGt                                            >  1:7990188/1‑151 (MQ=35)
                                                                                                            cACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTCGGTGAAGt                                            >  2:9045424/1‑151 (MQ=14)
                                                                                                            cACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGt                                            >  2:6891015/1‑151 (MQ=14)
                                                                                                             aCGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTg                                           <  1:5278235/151‑1 (MQ=35)
                                                                                                              cGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGc                                          >  1:1950698/1‑151 (MQ=14)
                                                                                                                gcgcTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCAc                                        >  2:3352642/1‑151 (MQ=14)
                                                                                                                 cgcTGTGGGTCATCGCGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACg                                       <  2:8591812/151‑1 (MQ=25)
                                                                                                                    tgtgGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGAGAAGTGCACGTCt                                    <  2:9564539/151‑1 (MQ=25)
                                                                                                                          tCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAgg                              <  2:7150711/151‑1 (MQ=35)
                                                                                                                               gAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAggtgggt                         <  2:4614841/151‑1 (MQ=35)
                                                                                                                               gAGGGCACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAggtgggt                         >  1:4614258/1‑151 (MQ=14)
                                                                                                                                   gCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAg                     <  2:10117724/151‑1 (MQ=35)
                                                                                                                                    cACGACGGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGCGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGt                    <  1:776299/151‑1 (MQ=14)
                                                                                                                                     acgacgGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCTACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTg                   <  2:5690433/151‑1 (MQ=25)
                                                                                                                                     acgacgGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTg                   <  2:5234550/151‑1 (MQ=14)
                                                                                                                                      cgacgGCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTgg                  <  2:10598995/151‑1 (MQ=14)
                                                                                                                                          ggCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGt              >  2:7048693/1‑151 (MQ=35)
                                                                                                                                          ggCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGt              <  1:8608817/151‑1 (MQ=14)
                                                                                                                                          ggCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGt              >  1:1610569/1‑151 (MQ=14)
                                                                                                                                           gCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGAGGGTCCAGGTGATCTTGCGGTAGGGGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGtt             >  1:8374531/1‑151 (MQ=25)
                                                                                                                                           gCCAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGtt             <  1:2136500/151‑1 (MQ=14)
                                                                                                                                            ccAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTg            <  2:2909126/151‑1 (MQ=35)
                                                                                                                                             cAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGc           >  2:8129041/1‑151 (MQ=35)
                                                                                                                                             cAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGc           >  2:8480538/1‑151 (MQ=35)
                                                                                                                                             cAGCGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGc           <  2:6527196/151‑1 (MQ=14)
                                                                                                                                                cGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGCCCg        <  1:259406/151‑1 (MQ=35)
                                                                                                                                                cGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGCCCg        <  2:3644463/151‑1 (MQ=14)
                                                                                                                                                cGAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGCCCg        >  1:3526538/1‑151 (MQ=14)
                                                                                                                                                 gAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGCCCgg       >  1:3825978/1‑151 (MQ=35)
                                                                                                                                                 gAGGCGGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGCCCgg       <  1:3614454/151‑1 (MQ=14)
                                                                                                                                                      |                                                                                                                                                      
GCGACCGATTCGATTGTCCTGCCTGGCGCGGGCGAGACAAGGGCGCGCAGGCTCGATCGGAAGCTGCCTACAACAAAGCGGCATCGACCCGGCGCGGGCTTTCCTGGCCACGGCGCTGTGGGTCATCGAGGGCACGACGGCCAGCGAGGCTGGCCGCGTGCGGCGGTCGCGGCTCAGCCGGCGACCGGGGAGCGCCAGTCGTCGGAGCCGGAGGTGCCGGAGACTTCGTGGGTCCAGGTGATCTTGCGGTAGGTGAAGTGCACGTCTTCCAGGTGGGTGAAGTGGGCGTTGCCCGGGTCCT  >  NC_008463/6205919‑6206219

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: