Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913730,7480TC30.8% 90.0 / 10.9 54I389T (ATC→ACC) rhsCrhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base T (35/1);  new base C (0/16);  total (35/19)
Fisher's exact test for biased strand distribution p-value = 1.64e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCAC  >  NC_000913/730526‑730995
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        
gCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTc                                                                                                                                                                                                                             >  1:573538/1‑251 (MQ=255)
              ggcgCACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGACATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑AACATCACCGACAGCCTT‑AACCGCCGTGAAGTCCTGCACACGCAggg                                                                                                                                                                                                               >  2:360127/1‑251 (MQ=255)
              ggcgAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAggg                                                                                                                                                                                                               >  1:411910/1‑251 (MQ=255)
                 gAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGGGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTg                                                                                                                                                                                                             >  1:446754/1‑250 (MQ=255)
                          gcggTTTAGGACCGCAGCAATACACGGTTGCGCAGTTTTACTTACGATAATAAATACCGGGGCCGGATGGTGGCGCCCCGTCCCCGGGGCCGACCGGAAATC‑TTTTACCGTTACGA‑CAGCGACGGGGGGGGGACAGAACACCTAACCCCGGCAGGTTTAATTTCCCCGTATCAGTATGAGAAGGACCGCTTC‑CCCACCACCGACAGCCTA‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCt                                                                                                                                                                                                   <  1:445855/251‑1 (MQ=255)
                                tATGACCGCAGCAATACACAGGAGCGAAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGAACGGACATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGAATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGACCACGCAGGGCGACGGCGGGCTGAAGCg                                                                                                                                                                                             >  2:513181/1‑251 (MQ=17)
                                 aTGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCgg                                                                                                                                                                                            >  2:524459/1‑251 (MQ=255)
                                        cagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGgtggtg                                                                                                                                                                                      >  2:360441/1‑250 (MQ=255)
                                            aaTACACCGGAGCGCAGCTTTACTTACGATGATACATAACGGGGCCAGATGCTGGCGCACCGTCACACGGGCCGACCGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCCTC‑CCCATCACCGACAGCCTG‑ACCCGCCGTGAAGTCCTGCACACGCCGGGCGAAGGCGGGCTGAAGCGGGTGGTGaaaa                                                                                                                                                                                  >  2:58280/1‑250 (MQ=14)
                                                     gTGCGAAGCTTTACCTACGATGATAAATACCGGCGCCGGATGGTGGCGCACCGTCACACGGGACGACCGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGCGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCCTC‑ACCATCACCGACAGCCTG‑AAACGCCGGGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAACAAGGAACCcgc                                                                                                                                                                        >  2:85032/1‑251 (MQ=14)
                                                       gcgcAGCTTTACTTACGATGATAAATACCGGGGACGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑ACCCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCgg                                                                                                                                                                      >  1:177903/1‑251 (MQ=255)
                                                                  cttacTATGATAAATACAGGGGCCGGACGGTAGCTCAAAGTCAAACAGGCAGAACGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGAGGGTAACAGAACAGATAAAACCGGCAGGCTTAAGGTACACGTATCAGTGTGAGACGGACCGCATC‑AGCATAACCGACAGCCTG‑AACCGCCGTGAAGACCTGCACACGCAGGGTGAAGGCGGTCTGAAGAGGGTGGTGAAACAGGAACACGAGGCCGGCAGCGTc                                                                                                                                                           >  2:518108/1‑251 (MQ=255)
                                                                     aggagaaaaaaacccGGGGCGGGAGGGGGGCGCACCGTCAAAGGGGCCGACGGGAAATC‑TGT‑ACCGGTAAGAACGGCGAGGGGGGGGGGAAAGAACAGAAAAACCCGGAAGGCTTAAGTAACACGTAAAAGTAGGAGAAGAACCGAATC‑CCAAACACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACg                                                                                                                                                        <  1:72616/238‑1 (MQ=9)
                                                                      cGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACCCGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGc                                                                                                                                                       >  1:250446/1‑251 (MQ=255)
                                                                      cGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGc                                                                                                                                                       >  1:98341/1‑251 (MQ=255)
                                                                                tACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGGGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGCCAGCGGCACGCAGAGTCAGtt                                                                                                                                             >  1:452451/1‑251 (MQ=255)
                                                                                cccgggggCGGGATGGGGGCGCCCCGTCACCGGGGCCGACCGAAAATT‑TGTTCCGGTTACGA‑CAGCGAGGGGCGGGTGAAAAAACAGCAAACCCCGGCGGGCTAAAGTAACCCGTATCAGTAGGAGAAGGACCGCACC‑CCAACCCCCGACAGCCGG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGtt                                                                                                                                             <  1:62458/247‑1 (MQ=9)
                                                                                 ccgggggcggggggggggCGCGCCGTCACGGGGGCCGACCGTAAATC‑TGTACCCGTTACGA‑CAGGGAGGGGGGGGTGAAAGAAAAGCAAACCCCGGCGGGCTTAAGTTACACGTATCAGTAGGAGAAGGACCGCATC‑ACCACCCCCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGtt                                                                                                                                             <  1:277004/234‑1 (MQ=14)
                                                                                    ggggCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGAc                                                                                                                                         >  1:551969/1‑251 (MQ=255)
                                                                                     gggcgggggggggccgcACCGTCAAGCGGGCAGACCGAAAATT‑GGTTCCCGTGACGA‑CGGCGAGGGGGGGGAGAAAGAACAGCAAACCCCGCCAGGCTTAAGTAACACGTATAAGTATGAGAAGGACCCCACC‑CCAACCACCGACAGCCTA‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCTTCACGCAGAGTCAGTTTGAcg                                                                                                                                        <  2:40526/237‑1 (MQ=9)
                                                                                          ggATGGTGGCGCACCGTCACACGGGCCGACCGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTg                                                                                                                                   >  1:488888/1‑251 (MQ=255)
                                                                                              ggtggCGCACCGTCACACGGGCCGACCGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGc                                                                                                                                >  2:378483/1‑250 (MQ=255)
                                                                                                tggCGCACCGTCACACGGGCCGACCGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGGGGGCAgg                                                                                                                             >  2:294728/1‑251 (MQ=255)
                                                                                                 ggccccccgacccaggggCCGCCCGGAAATC‑TGTTACCTTTACGA‑CAGCGACGGGGGGGTGACAGAACAGCAAAACCCGGCAGTCTTAATTTACACGAATCATAAGGAGAAGGCCCGCACC‑CCCACCCCCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAgg                                                                                                                             <  1:483042/235‑1 (MQ=14)
                                                                                                              gggggCCGAC‑GGAAATCCTGTCCCCGTGACGA‑CGGCGAGGGGGGGGTGAAAAAAACGCTAACCCCGGCGGGTTTAAGTTCCGCGTATAAGTATGAGAAGGCCCGCACC‑CCCACCCCCGACAGCCGG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGAc                                                                                                               <  2:219104/249‑1 (MQ=9)
                                                                                                               cGGGCCGACCGGAAATC‑TGTTACCGTTACGC‑CCGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACg                                                                                                              >  1:474462/1‑251 (MQ=255)
                                                                                                               cGGGCCGACCGGAAAAC‑TGATAACGATACGA‑CAGCGACGGGCGGGCGACAGAACAGCTAAACCCGGAAGGCTTAAATTAAACGTATCAGTATGAGAAGGCCCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGCGAAGGCCTGCACACGCAGGGTGAAGGCGGGCTGCAGCGGGTGGGGAAAAAGGAACACGCGGACGGAAGCGTAACGCAGAGTCAGTTTGAAGCGGGGGGACGGCCCAGGGCACAGACg                                                                                                              >  2:518932/1‑251 (MQ=255)
                                                                                                                    cGACCGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGCGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGc                                                                                                         >  1:28774/1‑251 (MQ=255)
                                                                                                                        cgaaaaTC‑TGTTCCCGTAACGA‑CAGCGACGGGGGGGTGACAGAACAGCAAACCCCGGAGGGCTTAATTTCAACGTAACAGTAGGAAAAGGCCCG‑ATCCACCACCCCCGACGGCCGG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCgca                                                                                                     <  1:383081/248‑1 (MQ=14)
                                                                                                                           aaaTC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCgcaggc                                                                                                  >  1:67233/1‑251 (MQ=255)
                                                                                                                                 tGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCGCCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCGGcaggcaggcc                                                                                             >  1:94367/1‑249 (MQ=255)
                                                                                                                                     ccccgtacggacgggggcgggcgggAGAAAGAACAGCAAACCCCGGCAGGCTAAAGTTACACGTATCAGTAGGAGAGGCACCGCACC‑CCCACCACCGACGGCGTA‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGacaac                                                                                          <  2:79027/234‑1 (MQ=14)
                                                                                                                                        ggtaccgtcgggggcgggggggTGAAAAAACAGCAAAACCCGCCAGGCTAAATTAACACTGACGAGTATGAGAGGGCCCGCACC‑ACCACCACCGACA‑CCTGTAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGAGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGa                                                                                       <  2:307686/237‑1 (MQ=14)
                                                                                                                                         ttACGA‑CAGCGACGGGCGGGTGACAGAACGGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAAt                                                                                     >  1:176265/1‑251 (MQ=255)
                                                                                                                                            cGA‑CAGCGAGGGGGGGGTAAAAGAACAGCTAAACCCGGCAGGCTAAAGTTCCCCGTATCAGTAGGAGAAGGACCGAACC‑ACCACCACCGACAGCTTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACa                                                                                  <  1:13949/251‑1 (MQ=14)
                                                                                                                                                cagcgaggggggtgggACAAAACAGCTAACCCCGGAAGGCTAAAGTTACGCATTTCAGTATGGAAAGGACCGCAGC‑GCCAACCCCGACAGCGTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTc                                                                               <  1:544491/236‑1 (MQ=14)
                                                                                                                                                 aGCGAAGGGCGGGTGACAGAACAGCTAAAACCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGAACGCATC‑AACATAACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGAGGTGACCAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCAGCAGGCAGGAACACAGAATACAGTcc                                                                              >  2:543974/1‑251 (MQ=14)
                                                                                                                                                  gcgaggggggggtaaaaGAACAGCAAACCCCGGCAGGTTAAATTTACACGTACCAGTAGGAGAAGGCCCCCATC‑CCCATCCCCGACAGCCGA‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCg                                                                             <  1:425532/235‑1 (MQ=14)
                                                                                                                                                       gggcgggTCAAAGAACAGCAAAACCCGGCAGGTTTAAGTTACACGTAACAGTATGAGAAGGACCGCATC‑CCCACCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCAAAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGAtgt                                                                        <  1:126746/251‑1 (MQ=14)
                                                                                                                                                           gggTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGCGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATgtgggg                                                                    >  2:80214/1‑249 (MQ=255)
                                                                                                                                                           gggTGAAAAAACAGCAAAACCCGGCAGGCTAAATTAACACGTAACAGTATGAAAAGGCCCGCACC‑ACCACCACCGAAAGCCGG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATgtggtg                                                                    <  1:360127/251‑1 (MQ=14)
                                                                                                                                                            ggggaaaGAACAGCTAACCCCGGCAGGCTTAAGTTCCACGTATCAGTATGAAAAGGCCCCCATC‑ACCACCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGa                                                                   <  1:524459/245‑1 (MQ=14)
                                                                                                                                                               gACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGAGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACgg                                                                >  1:581050/1‑251 (MQ=255)
                                                                                                                                                                    aaCAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCtc                                                           >  1:264062/1‑251 (MQ=255)
                                                                                                                                                                            aaCCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCAcgcg                                                   >  1:19663/1‑251 (MQ=255)
                                                                                                                                                                               ccGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCCGGACAACAGAAGACAGTCCGGATGTTGTGACGGGCCTCATCCCGCGCAt                                                >  2:259156/1‑251 (MQ=255)
                                                                                                                                                                                     ggcgaaaGTAACACGGATCAGTATGAAAGGGACCGCATC‑GCCTCCACGGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGGCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCGTCACGCGCATcaccac                                          <  1:482706/246‑1 (MQ=14)
                                                                                                                                                                                       cTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGc                                        >  1:280088/1‑251 (MQ=255)
                                                                                                                                                                                                 acGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGc                                >  1:192949/1‑249 (MQ=255)
                                                                                                                                                                                                     agctggaggagaGGGACCCCAAC‑ACCACCACCGACAGCAGG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGc                           <  1:379757/242‑1 (MQ=14)
                                                                                                                                                                                                       cAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATc                        >  1:324825/1‑251 (MQ=255)
                                                                                                                                                                                                         gtaggagaAGGCCCGCATC‑CCCACCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCgg                      <  1:21524/247‑1 (MQ=255)
                                                                                                                                                                                                                  ggACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGCAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGGGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCTGGGCATCGGCGTTTTACt             >  2:53875/1‑251 (MQ=255)
                                                                                                                                                                                                                       gCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAac        >  1:122863/1‑251 (MQ=255)
                                                                                                                                                                                                                        aatc‑accaaccccGACAGCCGA‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGGCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAacc       <  2:385733/250‑1 (MQ=255)
                                                                                                                                                                                                                                atcaccGACAGCCTG‑AACCGCCGTGAAGTCCTGCACAAGCAGGGCGAAGGCGGGCTGCAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGAGACGGGCATCATCAAGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACcaccac  >  2:578967/1‑249 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        
GCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATC‑TGTTACCGTTACGA‑CAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATC‑ACCATCACCGACAGCCTG‑AACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCAC  >  NC_000913/730526‑730995

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