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.

GCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCAC  >  NC_000913/730526‑730995
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        
gCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTc                                                                                                                                                                                                                             >  1:573538/1‑251 (MQ=255)
              ggcgCACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGACATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCAACATCACCGACAGCCTTAACCGCCGTGAAGTCCTGCACACGCAggg                                                                                                                                                                                                               >  2:360127/1‑251 (MQ=255)
              ggcgAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAggg                                                                                                                                                                                                               >  1:411910/1‑251 (MQ=255)
                 gAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGGGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTg                                                                                                                                                                                                             >  1:446754/1‑250 (MQ=255)
                          gcggTTTAGGACCGCAGCAATACACGGTTGCGCAGTTTTACTTACGATAATAAATACCGGGGCCGGATGGTGGCGCCCCGTCCCCGGGGCCGACCGGAAATCTTTTACCGTTACGACAGCGACGGGGGGGGGACAGAACACCTAACCCCGGCAGGTTTAATTTCCCCGTATCAGTATGAGAAGGACCGCTTCCCCACCACCGACAGCCTAAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCt                                                                                                                                                                                                   <  1:445855/251‑1 (MQ=255)
                                tATGACCGCAGCAATACACAGGAGCGAAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGAACGGACATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGAATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGACCACGCAGGGCGACGGCGGGCTGAAGCg                                                                                                                                                                                             >  2:513181/1‑251 (MQ=17)
                                 aTGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCgg                                                                                                                                                                                            >  2:524459/1‑251 (MQ=255)
                                        cagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGgtggtg                                                                                                                                                                                      >  2:360441/1‑250 (MQ=255)
                                            aaTACACCGGAGCGCAGCTTTACTTACGATGATACATAACGGGGCCAGATGCTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCCTCCCCATCACCGACAGCCTGACCCGCCGTGAAGTCCTGCACACGCCGGGCGAAGGCGGGCTGAAGCGGGTGGTGaaaa                                                                                                                                                                                  >  2:58280/1‑250 (MQ=14)
                                                     gTGCGAAGCTTTACCTACGATGATAAATACCGGCGCCGGATGGTGGCGCACCGTCACACGGGACGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGCGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCCTCACCATCACCGACAGCCTGAAACGCCGGGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAACAAGGAACCcgc                                                                                                                                                                        >  2:85032/1‑251 (MQ=14)
                                                       gcgcAGCTTTACTTACGATGATAAATACCGGGGACGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGACCCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCgg                                                                                                                                                                      >  1:177903/1‑251 (MQ=255)
                                                                  cttacTATGATAAATACAGGGGCCGGACGGTAGCTCAAAGTCAAACAGGCAGAACGGAAATCTGTTACCGTTACGACAGCGACGGGAGGGTAACAGAACAGATAAAACCGGCAGGCTTAAGGTACACGTATCAGTGTGAGACGGACCGCATCAGCATAACCGACAGCCTGAACCGCCGTGAAGACCTGCACACGCAGGGTGAAGGCGGTCTGAAGAGGGTGGTGAAACAGGAACACGAGGCCGGCAGCGTc                                                                                                                                                           >  2:518108/1‑251 (MQ=255)
                                                                     aggagaaaaaaacccGGGGCGGGAGGGGGGCGCACCGTCAAAGGGGCCGACGGGAAATCTGTACCGGTAAGAACGGCGAGGGGGGGGGGAAAGAACAGAAAAACCCGGAAGGCTTAAGTAACACGTAAAAGTAGGAGAAGAACCGAATCCCAAACACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACg                                                                                                                                                        <  1:72616/238‑1 (MQ=9)
                                                                      cGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACCCGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGc                                                                                                                                                       >  1:250446/1‑251 (MQ=255)
                                                                      cGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGc                                                                                                                                                       >  1:98341/1‑251 (MQ=255)
                                                                                tACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGGGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGCCAGCGGCACGCAGAGTCAGtt                                                                                                                                             >  1:452451/1‑251 (MQ=255)
                                                                                cccgggggCGGGATGGGGGCGCCCCGTCACCGGGGCCGACCGAAAATTTGTTCCGGTTACGACAGCGAGGGGCGGGTGAAAAAACAGCAAACCCCGGCGGGCTAAAGTAACCCGTATCAGTAGGAGAAGGACCGCACCCCAACCCCCGACAGCCGGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGtt                                                                                                                                             <  1:62458/247‑1 (MQ=9)
                                                                                 ccgggggcggggggggggCGCGCCGTCACGGGGGCCGACCGTAAATCTGTACCCGTTACGACAGGGAGGGGGGGGTGAAAGAAAAGCAAACCCCGGCGGGCTTAAGTTACACGTATCAGTAGGAGAAGGACCGCATCACCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGtt                                                                                                                                             <  1:277004/234‑1 (MQ=14)
                                                                                    ggggCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGAc                                                                                                                                         >  1:551969/1‑251 (MQ=255)
                                                                                     gggcgggggggggccgcACCGTCAAGCGGGCAGACCGAAAATTGGTTCCCGTGACGACGGCGAGGGGGGGGAGAAAGAACAGCAAACCCCGCCAGGCTTAAGTAACACGTATAAGTATGAGAAGGACCCCACCCCAACCACCGACAGCCTAAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCTTCACGCAGAGTCAGTTTGAcg                                                                                                                                        <  2:40526/237‑1 (MQ=9)
                                                                                          ggATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTg                                                                                                                                   >  1:488888/1‑251 (MQ=255)
                                                                                              ggtggCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGc                                                                                                                                >  2:378483/1‑250 (MQ=255)
                                                                                                tggCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGGGGGCAgg                                                                                                                             >  2:294728/1‑251 (MQ=255)
                                                                                                 ggccccccgacccaggggCCGCCCGGAAATCTGTTACCTTTACGACAGCGACGGGGGGGTGACAGAACAGCAAAACCCGGCAGTCTTAATTTACACGAATCATAAGGAGAAGGCCCGCACCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAgg                                                                                                                             <  1:483042/235‑1 (MQ=14)
                                                                                                              gggggCCGACGGAAATCCTGTCCCCGTGACGACGGCGAGGGGGGGGTGAAAAAAACGCTAACCCCGGCGGGTTTAAGTTCCGCGTATAAGTATGAGAAGGCCCGCACCCCCACCCCCGACAGCCGGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGAc                                                                                                               <  2:219104/249‑1 (MQ=9)
                                                                                                               cGGGCCGACCGGAAATCTGTTACCGTTACGCCCGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACg                                                                                                              >  1:474462/1‑251 (MQ=255)
                                                                                                               cGGGCCGACCGGAAAACTGATAACGATACGACAGCGACGGGCGGGCGACAGAACAGCTAAACCCGGAAGGCTTAAATTAAACGTATCAGTATGAGAAGGCCCGCATCACCATCACCGACAGCCTGAACCGCCGCGAAGGCCTGCACACGCAGGGTGAAGGCGGGCTGCAGCGGGTGGGGAAAAAGGAACACGCGGACGGAAGCGTAACGCAGAGTCAGTTTGAAGCGGGGGGACGGCCCAGGGCACAGACg                                                                                                              >  2:518932/1‑251 (MQ=255)
                                                                                                                    cGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGCGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGc                                                                                                         >  1:28774/1‑251 (MQ=255)
                                                                                                                        cgaaaaTCTGTTCCCGTAACGACAGCGACGGGGGGGTGACAGAACAGCAAACCCCGGAGGGCTTAATTTCAACGTAACAGTAGGAAAAGGCCCGATCCACCACCCCCGACGGCCGGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCgca                                                                                                     <  1:383081/248‑1 (MQ=14)
                                                                                                                           aaaTCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCgcaggc                                                                                                  >  1:67233/1‑251 (MQ=255)
                                                                                                                                 tGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCGCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCGGcaggcaggcc                                                                                             >  1:94367/1‑249 (MQ=255)
                                                                                                                                     ccccgtacggacgggggcgggcgggAGAAAGAACAGCAAACCCCGGCAGGCTAAAGTTACACGTATCAGTAGGAGAGGCACCGCACCCCCACCACCGACGGCGTAAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGacaac                                                                                          <  2:79027/234‑1 (MQ=14)
                                                                                                                                        ggtaccgtcgggggcgggggggTGAAAAAACAGCAAAACCCGCCAGGCTAAATTAACACTGACGAGTATGAGAGGGCCCGCACCACCACCACCGACACCTGTAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGAGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGa                                                                                       <  2:307686/237‑1 (MQ=14)
                                                                                                                                         ttACGACAGCGACGGGCGGGTGACAGAACGGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAAt                                                                                     >  1:176265/1‑251 (MQ=255)
                                                                                                                                            cGACAGCGAGGGGGGGGTAAAAGAACAGCTAAACCCGGCAGGCTAAAGTTCCCCGTATCAGTAGGAGAAGGACCGAACCACCACCACCGACAGCTTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACa                                                                                  <  1:13949/251‑1 (MQ=14)
                                                                                                                                                cagcgaggggggtgggACAAAACAGCTAACCCCGGAAGGCTAAAGTTACGCATTTCAGTATGGAAAGGACCGCAGCGCCAACCCCGACAGCGTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTc                                                                               <  1:544491/236‑1 (MQ=14)
                                                                                                                                                 aGCGAAGGGCGGGTGACAGAACAGCTAAAACCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGAACGCATCAACATAACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGAGGTGACCAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCAGCAGGCAGGAACACAGAATACAGTcc                                                                              >  2:543974/1‑251 (MQ=14)
                                                                                                                                                  gcgaggggggggtaaaaGAACAGCAAACCCCGGCAGGTTAAATTTACACGTACCAGTAGGAGAAGGCCCCCATCCCCATCCCCGACAGCCGAAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCg                                                                             <  1:425532/235‑1 (MQ=14)
                                                                                                                                                       gggcgggTCAAAGAACAGCAAAACCCGGCAGGTTTAAGTTACACGTAACAGTATGAGAAGGACCGCATCCCCACCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCAAAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGAtgt                                                                        <  1:126746/251‑1 (MQ=14)
                                                                                                                                                           gggTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGCGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATgtgggg                                                                    >  2:80214/1‑249 (MQ=255)
                                                                                                                                                           gggTGAAAAAACAGCAAAACCCGGCAGGCTAAATTAACACGTAACAGTATGAAAAGGCCCGCACCACCACCACCGAAAGCCGGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATgtggtg                                                                    <  1:360127/251‑1 (MQ=14)
                                                                                                                                                            ggggaaaGAACAGCTAACCCCGGCAGGCTTAAGTTCCACGTATCAGTATGAAAAGGCCCCCATCACCACCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGa                                                                   <  1:524459/245‑1 (MQ=14)
                                                                                                                                                               gACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGAGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACgg                                                                >  1:581050/1‑251 (MQ=255)
                                                                                                                                                                    aaCAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCtc                                                           >  1:264062/1‑251 (MQ=255)
                                                                                                                                                                            aaCCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCAcgcg                                                   >  1:19663/1‑251 (MQ=255)
                                                                                                                                                                               ccGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCCGGACAACAGAAGACAGTCCGGATGTTGTGACGGGCCTCATCCCGCGCAt                                                >  2:259156/1‑251 (MQ=255)
                                                                                                                                                                                     ggcgaaaGTAACACGGATCAGTATGAAAGGGACCGCATCGCCTCCACGGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGGCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCGTCACGCGCATcaccac                                          <  1:482706/246‑1 (MQ=14)
                                                                                                                                                                                       cTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGc                                        >  1:280088/1‑251 (MQ=255)
                                                                                                                                                                                                 acGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGc                                >  1:192949/1‑249 (MQ=255)
                                                                                                                                                                                                     agctggaggagaGGGACCCCAACACCACCACCGACAGCAGGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGc                           <  1:379757/242‑1 (MQ=14)
                                                                                                                                                                                                       cAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATc                        >  1:324825/1‑251 (MQ=255)
                                                                                                                                                                                                         gtaggagaAGGCCCGCATCCCCACCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCgg                      <  1:21524/247‑1 (MQ=255)
                                                                                                                                                                                                                  ggACCGCATCACCATCACCGACAGCCTGAACCGCCGTGCAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGGGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCTGGGCATCGGCGTTTTACt             >  2:53875/1‑251 (MQ=255)
                                                                                                                                                                                                                       gCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAac        >  1:122863/1‑251 (MQ=255)
                                                                                                                                                                                                                        aatcaccaaccccGACAGCCGAAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGGCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAacc       <  2:385733/250‑1 (MQ=255)
                                                                                                                                                                                                                                atcaccGACAGCCTGAACCGCCGTGAAGTCCTGCACAAGCAGGGCGAAGGCGGGCTGCAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGAGACGGGCATCATCAAGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACcaccac  >  2:578967/1‑249 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        
GCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCAC  >  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: