Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913732,7650CT28.2% 46.2 / 20.8 39T1061T (ACC→ACT) rhsCrhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base C (12/16);  new base T (3/8);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATAC  >  NC_000913/732518‑732993
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                      
cgcgGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGcc                                                                                                                                                                                                                                     <  1:82048/251‑1 (MQ=21)
   ggACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTCCGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCAcc                                                                                                                                                                                                                                  <  1:5248/251‑1 (MQ=21)
   ggACACAATATGCAGAGCCGCGGGCCGAAAGTCGCTATCTTTACGCCCCGCTGGGCCGCAGGGTGGCAAAAGGGGTATGGCGGGGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGTTTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCAcc                                                                                                                                                                                                                                  <  1:208125/251‑1 (MQ=255)
    gacccaagatgcggtgCCGCTGGTCGAAATTCGCGATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAGCGGGTATGGCGGCGGGAACGGGACCTGACGGGCTGGATGTCGCTGCCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCGTACCACGATACAGAACGAAAGATCCCGCATCCAGACGATTTATCAGCCGGGAAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCg                                                                                                                                                                                                                                 <  2:215630/236‑1 (MQ=11)
    gACACAAGGTGCAGAGCAGTTGGTGGAAAGTCGCTATCTTTACGACCTGCTGGGCCGTAGGGTGGAAAAACGTGTAGGGCTGCGTCAACGGGACCTGACGGTCTGGATGTCGCTGTCACGGAAACCGCAAG‑TGACTTGGTA‑CGGT‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTTTCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCg                                                                                                                                                                                                                                 <  2:215637/251‑1 (MQ=11)
        caATAGGCAGGGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTAGGGCGGCGGGAACGGGACCTGACGGGCGGGATGTCGCTGTCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGa                                                                                                                                                                                                                             <  2:140825/251‑1 (MQ=255)
              gCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGc                                                                                                                                                                                                                       <  1:231986/251‑1 (MQ=21)
                     cGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGCCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACg                                                                                                                                                                                                                <  1:59706/251‑1 (MQ=21)
                      gCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGc                                                                                                                                                                                                               >  1:199453/1‑251 (MQ=34)
                         ggTCGAAGGTGGTTATCTTTACGACCCGCTGGGCCGCGGGGGGGCAAAACGGGTAGGGCGGCGGGAACGGGACCGGACGGGCGGGATGTCGCTGTCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAgc                                                                                                                                                                                                            <  2:223476/251‑1 (MQ=11)
                                      tATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGc                                                                                                                                                                                              >  1:68939/1‑251 (MQ=18)
                                             aCGACCCGCGGGCCCGCAGGGTGACAAACGGGGCAAGGCGGCGTGAACGGGACCGGACGGGCTGTAGGTCGCTGTCACGTAAACCCCAAG‑TGACCGGTTA‑CGGC‑TGGAACGGCGACCGCCTGACCACGATACAGAACGACAGAGCCCGCATCCAGACGAGTGAGCAGCGGGGGAGCTGCACGCCACTCAGCGGTGTTGAAACTGCCACCTGTGAGCAGGCGAAAACGCAGCGCCGCATCC‑GGGCGGaggc                                                                                                                                                                                        <  2:124725/250‑4 (MQ=2)
                                               gACCCGCTGGGCCGCAGGGTGGCCAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGATGGATGTCGCTGTCCCGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGAGGGCGACCGCCTGACCCCGAAACAGAACGACAGAACCCGGATCCCGACGATTTATCAGCCGGGGAGCTTCACGCCAATTATCAGGGTTGACAATGCCAACGGGCAGCTCGCGGAAACGCAGCGCCGCAGCC‑TGGCggttggccc                                                                                                                                                                                     >  2:243916/1‑244 (MQ=11)
                                                cccagcggggCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGACCGGGACCTGACGGGCGGGATGTCGCTGTCACGGAAACCGCAAG‑TGACCTGGTC‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGACCCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATGCCCtt                                                                                                                                                                                    <  2:54614/244‑1 (MQ=11)
                                                aCCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGAACCTGACGGGCTGGATGTCGCTGTCACGGCAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGAACGCCTGACCAAGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGATGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCtt                                                                                                                                                                                    >  1:132318/1‑251 (MQ=255)
                                                  ccGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTc                                                                                                                                                                                   >  1:7641/1‑250 (MQ=33)
                                                     cTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGCACGGGACCTGCCGGGCTGGCTGTCGCTGTCACGGAAACCTCAAG‑TGACCTGGTA‑CGGC‑TGGGTCGGCGACCGCCTGCCCACGATACAGAAAGACAGCACCCCCCTCCCGACGATTTCTCAGCCTGGGAGCTTCTCGCCCCTCATCCGGGTTGTAACCTCCACCGGTTAGCTGTCGCACACTCAGCCCCGCCTCC‑TCGCGGATTCCCttttgcc                                                                                                                                                                               >  2:79962/1‑246 (MQ=255)
                                                       gggCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGTCAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTCAGCAGt                                                                                                                                                                             <  1:55757/251‑1 (MQ=18)
                                                         gCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATGCCCTTCAGCAGTcc                                                                                                                                                                           <  1:253788/251‑1 (MQ=18)
                                                         cccGCAGGGTGGCAAAGCGGGTAGGTCGGCTGGACCGGGACCTGCCGGGCTGGATGTCGCTGTCACGGAACCCGCAAG‑TGCCCTGTTA‑CGGC‑TGGGACGGGGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTCAGCAGTcc                                                                                                                                                                           <  2:60479/250‑1 (MQ=11)
                                                            gCAGGGTGGCAAACCGGGTATGGGGGCGGGAGCGGGACCGGCGGGGCGGGTTGTCGCTGTCACGGACACCGCAAG‑TGACC‑GGTAACGGC‑TGGGACGGCGACCGCCTGGCCACGATACAGACCGACAGACCCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTCAGCAGTCCGGt                                                                                                                                                                        <  2:12534/251‑1 (MQ=11)
                                                                   ggCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTCAGCAGTCCGGTGGCGAAg                                                                                                                                                                 <  1:148106/251‑1 (MQ=34)
                                                                    gaaaaccGGGTATGGCGGCGTGAACTGGACCTTGCGGCCTGGATGTCCCTGTCACGGAAACCGCAAG‑TGACCTGTTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGAGACATAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTTGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGa                                                                                                                                                                <  1:259554/245‑1 (MQ=11)
                                                                        aaCGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACggg                                                                                                                                                            >  1:211679/1‑250 (MQ=9)
                                                                                        tGAACGGGACCGTACTGGGTGAATGCCGCTGTCACCGAACCCGCGAC‑TACCCTGGTA‑CGGC‑TGGGACGGCAACCGCCTGGCCACGATACAGAACGACAGATCCCTCATCCAGACGATTTAGCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAATCCGCCACCGGTGAGCTGGCGAAAACGCAGCGGCGCTGCC‑TGGCGGATACCCTTCAGCAGGCCGGTGGCGAAGACGGTGGCAGTGTCgtg‑ta                                                                                                                                             <  2:22232/249‑2 (MQ=11)
                                                                                         gAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGATG‑TTccc                                                                                                                                          <  1:248262/251‑1 (MQ=14)
                                                                                                                gTCGCTGTCACGAAAACC‑CAAGGTGCCCTGGTA‑CGGC‑TGGGCCGGCGCCCGCCTGCCCACGATACAGAACGACAGACCCCGCATCCAGACGTTTTACCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATCCCCTTCAGCAGTCCGGTGGCGAAGACGGGGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCg                                                                                                                   <  1:190332/251‑1 (MQ=11)
                                                                                                                   cctgccACGGAACCCGCAAG‑TGACCTGGTA‑CGGC‑GGGGCCGGCGACCGCCTGCCCACGATACAGACCGCCAGAACCCGCACCCAGACGTTTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGAcc                                                                                                                <  1:201997/246‑1 (MQ=11)
                                                                                                                     tGTCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCgg                                                                                                              >  1:128565/1‑251 (MQ=31)
                                                                                                                           ggAAACCGCAAG‑TGCCGTGGT‑‑CGGCCTGGAACGGCGCCCGCCTGACCAGGATACAGAACGAAAAAACCCGCATCCAGACAATTTACCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAACCCGCCCCCGGTGAGCTGGGGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGTTT‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGaa                                                                                                        <  1:106725/251‑1 (MQ=11)
                                                                                                                                            cTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTg                                                                                        >  1:46562/1‑251 (MQ=31)
                                                                                                                                                            aCGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtgagtg                                                                            >  1:4249/1‑249 (MQ=31)
                                                                                                                                                                                aTACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGGG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGAGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCtggctggc                                                      >  2:12740/1‑251 (MQ=14)
                                                                                                                                                                                  aCAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCAt                                                    >  1:49995/1‑251 (MQ=31)
                                                                                                                                                                                           cAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGcc                                           >  2:101675/1‑251 (MQ=31)
                                                                                                                                                                                                                   ccccccGGGGACCTCCACGCCACATATCAGGGTGCAACCCGCCACCGGTGAGGGGGCGAAAGCCAAGCGCCTCAGCCGTGG‑GGCTCCCCTTCAGCGGTCGGGGGGCGAGGGCGTGGGCGGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATTCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGcc                   <  2:130971/247‑1 (MQ=2)
                                                                                                                                                                                                                       ccGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGcc                   >  1:181829/1‑247 (MQ=34)
                                                                                                                                                                                                                             gccgcccccccACTCAGCGGGGTTGAACCCGCCCCCGGTGAGCTGGCGAAAACGCACCGCCCCAGCC‑GGGCGGATACCCTTCACCAGTCGGGTGGCGAAGACGGTGGCAG‑GTGGGGGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGAccc         <  1:34205/242‑1 (MQ=11)
                                                                                                                                                                                                                                     gCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATac  >  2:222836/1‑250 (MQ=37)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                      
CGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAG‑TGACCTGGTA‑CGGC‑TGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCC‑TGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTG‑TTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATAC  >  NC_000913/732518‑732993

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