Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,906,7690AC27.4% 111.8 / 11.5 62Y98D (TAT→GAT) rlmA23S rRNA m(1)G745 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (30/15);  new base C (0/17);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 9.27e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTT  >  NC_000913/1906520‑1907016
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        
gTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATAc                                                                                                                                                                                                                                                         <  1:232478/251‑1 (MQ=255)
   ccGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCt                                                                                                                                                                                                                                                      <  1:324615/251‑1 (MQ=255)
     ggCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTc                                                                                                                                                                                                                                                    >  2:142720/1‑251 (MQ=255)
      gCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTca                                                                                                                                                                                                                                                   <  1:101118/251‑1 (MQ=255)
           gCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACcaca                                                                                                                                                                                                                                              <  1:118622/251‑1 (MQ=255)
             aGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATACTTCTCCTTTACACGGCGCGTAACTACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCCAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGCAACATCCCGACCAACC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACACTTCACCACAGc                                                                                                                                                                                                                                            >  2:66737/1‑251 (MQ=255)
             aGTAATGACCCAGCCGCCGGGCTTCACTACTCGGGCTAATTCTTCTGCTTTACACGGCGCGTACCTACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGGGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGATACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGAGTTCGGGCAACGCAGCGGCAAATGCGTG‑‑GTGGTAATACCCTTCACCACAGc                                                                                                                                                                                                                                            >  2:48788/1‑251 (MQ=255)
              gTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAATCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGcc                                                                                                                                                                                                                                           <  1:367372/251‑1 (MQ=255)
                              cgGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCAcc‑gcc                                                                                                                                                                                                                          <  1:158978/251‑1 (MQ=255)
                                       ctactCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATACCCAGCACC‑GCCGTGGCCtt                                                                                                                                                                                                                  >  1:133891/1‑250 (MQ=255)
                                        tactCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATCCGTATTATGGCGTCCATACTGGTATCGGAAACCGGCACACGGTGGCTGGAAGCGACACACAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTat                                                                                                                                                                                                                >  1:258280/1‑251 (MQ=255)
                                            cGTGCAAATTCTTCTTCTTTACACGGCGCGTAAATACGTTTTATGGCGTCCATCCTGGTATCGGAAACCGGCAACCGGTGGCTGAAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTCCGCCGCCGCTTTTATCGCTCCCT‑TCGAAACACCCAGACCAAAC‑GTGG‑TGATTTCGGGCACCGCTTCGGCAAATGCGGG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTatca                                                                                                                                                                                                              <  1:380390/249‑1 (MQ=255)
                                             ggccaaaTTCTTCTGCTTGACCCGGCGCGAAAATGCGTATTATGGCGCCCAGACTGGTAGCGGAAAACGGCAAACGGGGGCTGGACGCAACAAAAAAAGTGACCC‑GCGGAGGGGGTTTCCCCGCCGCTTTTATCGCTACCT‑TCGAAACACCCAGCCCAAAC‑GGGT‑TGATTTCGGGCAACGCACCGGCAAAGGCGGG‑TGGG‑TAATACCCTCAACCCCAGCCAATATCCAGCACC‑GCCGTGGCCTTAtcatca                                                                                                                                                                                                           <  2:294625/246‑1 (MQ=255)
                                                  aaTTCTTCTGCTTTACCCGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAACCGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCg                                                                                                                                                                                                      <  1:129362/251‑1 (MQ=255)
                                                        tctGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCt                                                                                                                                                                                                >  2:17503/1‑251 (MQ=255)
                                                        tctGCTTCACACGGCGCGTAAATACGTATTAGGGCGTCCATACGGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGGGACCG‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAGC‑GGGT‑TGTTTTCGGGCAACGCATGGGCAAATGCGTT‑GGGG‑TAATCCCTTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTcc                                                                                                                                                                                                  <  1:39444/249‑1 (MQ=255)
                                                          cgcttaccccGGCGCGTAAATCCGTATTATGGCGCCCATACGGGTATCGAAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGCCCT‑GCGAATAGCGTTCCGCCGCCGCTTTTATCGCTCCCT‑TCAAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATCCCCTTCACCACAGCCAATACCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCa                                                                                                                                                                                              <  1:386232/242‑1 (MQ=255)
                                                                 cacGGCGCGTAAATACGTATTATGGCGTCCATACTGGTACCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGCCCG‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTCCCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAAGGCGTG‑TGTG‑TAATCCCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCAGCAAGCCGTTCCCTCAGTTGGGc                                                                                                                                                                                       <  1:452496/251‑1 (MQ=255)
                                                                         gTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTg                                                                                                                                                                               >  1:254200/1‑251 (MQ=255)
                                                                              aactttttaggcgccccagacgggTATCGAAAAACGGCAAAGGGTGGCTGAAAGCAACAAAAAAAGGGCCCG‑GCGGTTGGCTTTTCGCCGCCGTTTTTATCGCTCCCT‑ACAAACCATCCAGCCCAAAC‑GGGG‑TGATTTCGGGCAACGCATCGCCAAATG‑GTGGTGTG‑TAATCCCCTTCACCACGGCAAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCa                                                                                                                                                                          <  2:21144/229‑1 (MQ=255)
                                                                               aCGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCa                                                                                                                                                                          >  1:21144/1‑250 (MQ=255)
                                                                                 gTATTATGGCGACCACAATGGTATCGGAAAAAGGCAAACGGTGGCTGGAAGCGCCACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCCACCT‑TCGACACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCAAC‑GACGTGGCCTTATCATCATGCCGTTCCCTCAGTTGGGCGACCACTGCATCACGc                                                                                                                                                                       >  2:320302/1‑251 (MQ=255)
                                                                                      atGGCGTCCATACTGGTATCGGAAAACGGCAACCGGTGGCTGGAAGCGACAAAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGT‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATCCCCTTCCCCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCg                                                                                                                                                                   <  1:198859/250‑1 (MQ=255)
                                                                                       tGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGc                                                                                                                                                                 >  1:391461/1‑251 (MQ=255)
                                                                                            tCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGAt                                                                                                                                                             >  1:74138/1‑250 (MQ=255)
                                                                                              cATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCACCGGCTGATAAt                                                                                                                                                          >  2:8420/1‑251 (MQ=255)
                                                                                                        cggaaaaaCGGCAACGGGGGGCGGGAAGCAACA‑AAAAAGTGACCCCGCGGAGAGCGTTTCGCGCCCGCTTTTACCGCTACCT‑TCGAAACACCCAG‑CCAAACCGGGG‑TAATTCGGGGGAACGCAGCGGCAAATGCGGG‑GGTG‑TAATCCCCTTCACCACGGCCAAAATCCAGCCCC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATc                                                                                                                                                <  2:154962/247‑1 (MQ=255)
                                                                                                           aaaaaCGGCAAACGGGGGTTGGAGGCAACCAAAAAAGGGACCG‑GGGGAGAGCGTTCCGCGCCCGCGTTTACCGCGCCCT‑CCAAAACACCCAGAC‑AAAA‑GTGGGTGTTTGCGGACAACGCCTCGGCAAATGCGGG‑TGTG‑TAATCCCCTTCCCCACAGCCAATAGCCAGCTCC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCATTTGGGCGACAATTGCGTCACGCAGCGGCTGATAATGTCCGGCATCTaa                                                                                                                                             <  2:254200/250‑1 (MQ=255)
                                                                                                            aaaaCGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTa                                                                                                                                              >  1:406088/1‑249 (MQ=255)
                                                                                                               aCGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAatat                                                                                                                                                                                                                                      <  1:291310/159‑1 (MQ=255)
                                                                                                               aCGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAatat                                                                                                                                                                                                                                      >  2:291310/1‑159 (MQ=255)
                                                                                                                gggaaacgggtgtgtgaaagggACAAAAAAAGGGCCCG‑GCGGATAGGGTTTCGCCGCCGCTTTTATCGCTCCCT‑GCGAAAAACCAACACAAAAC‑GGGG‑GGATTTGGGGCACCGCATCGGAAAAGGCGTG‑GGTT‑AAATCCCCTCCCCCCCAGCCAATACCCACCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATg                                                                                                                                        <  2:524025/230‑1 (MQ=255)
                                                                                                                 ggCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTACCGCTACCT‑TCGAAACATCCAGACCAAAC‑GGGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATCCCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATgc                                                                                                                                       <  1:142720/251‑1 (MQ=255)
                                                                                                                   cAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCCAATAGCGATTCGCAGACGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GAGG‑GGATTTCGGGCAACCCATCGGCAAATCCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTCAACATCAAGCCGTTCCCTCAGTTGGGCAACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATgcgc                                                                                                                                     >  1:506986/1‑251 (MQ=255)
                                                                                                                    aaaCGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGcgcg                                                                                                                                    >  1:105319/1‑251 (MQ=255)
                                                                                                                        ggTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGt                                                                                                                                >  1:172691/1‑251 (MQ=255)
                                                                                                                             cTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCATAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTg                                                                                                                           >  1:192195/1‑251 (MQ=255)
                                                                                                                                  aGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGcatc                                                                                                                       >  2:88299/1‑250 (MQ=255)
                                                                                                                                        caaaaaaagtcccctgcgaaaggcttttCGCCGCCGTTTTTATCGCTCCCT‑TCAAAACACCCAGACAAAAC‑TGGG‑TGATTCCGGGCAACGCATCGGCAAATGGGTG‑TGTG‑TAATCCCCTTCCCCCCAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTcc                                                                                                                 <  2:391461/226‑1 (MQ=255)
                                                                                                                                              agttctccggcggagaggttttCGCCCCCGCTTTTCTCTCTCCCT‑CGGAAAAAGCCAGACAAAAC‑GGGG‑TGATGGCGGACACCCCACCGGCAAAGGCGTG‑GGGG‑TAAGCCCCTACCCCACAGCAAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTg                                                                                                           <  2:157909/232‑1 (MQ=255)
                                                                                                                                                           aTAGCGTTCCGCCGCCGCTTTTATCCCTACCT‑TCGAAACATCCAGACCAACC‑GTGG‑TGATTCCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCAc                                                                                              <  1:146521/251‑1 (MQ=255)
                                                                                                                                                                   ccgccgccgcTTTGACCTCTCCCT‑TCAAAACAGCCAGACAAAAC‑GTGT‑TGTTTTCGGGAAACGCATGGGCAAATGCGTG‑GTTG‑TAATACCCTCCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCg                                                                                       <  2:383222/249‑1 (MQ=255)
                                                                                                                                                                       cgcctgttttACCGCTACCT‑GCGAAACCACCGGCCCAAAC‑GTGG‑TAATTGCGGGCAACGCTTCGGCAAATGCGTG‑TGTG‑TAATCCCCTTCCCCCCAGCCAATACCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTAt                                                                                  <  2:172691/245‑1 (MQ=255)
                                                                                                                                                                        gccgcTTTTATCGCTCCCT‑TCGACACACCGACACAAAAA‑GTGG‑TGATATCGGGCCACGCATCGGCAAATGCGTA‑AGTG‑TAATACCCTTCACCACAGCCAATATACCGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCAGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACgag                                                                                           >  2:218220/1‑241 (MQ=255)
                                                                                                                                                                        gccgcTTTTAGCGCTACCT‑GCTAAACATCCAGACCAAAC‑GTGG‑TGATTCCGGGCAGCGAATGGGCAAATGCGTG‑TGTG‑TAATCCCCTTCACCCCAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGAGAATTGCCGGCATCTAAGAAGGCGCTGCTTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACgag                                                                                           <  1:218220/241‑1 (MQ=255)
                                                                                                                                                                         cccttgttATCGCTACCTCTC‑AACCATCCAGACCAACT‑GTGG‑TGATTTCGGGCAACCCATCTGCAAATGCGTT‑TGTG‑TAATCCCCTCCACCCCACCCAATATCCA‑CACCCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGTGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCTTCATTTCCGCGCTGTCGCCCGGATCACGAGCCCGTTTATGc                                                                                <  1:48788/245‑1 (MQ=255)
                                                                                                                                                                          cgcTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑AGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGACTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGACTGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCt                                                                               >  1:48299/1‑251 (MQ=255)
                                                                                                                                                                                    cgctaccctagccaCATCAAGCCCAAAC‑GTGG‑TGCTTTCGGGCAACGCATCGGCAAATACGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCCACCCGCCGCTCACTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCTTCTAAGAATGCGCGGCGTGCTTGACTCCTTTCCGCGCTGTCGACCGGAACACGAGACCGTTTATGCTGACCGGgc                                                                       >  2:210306/14‑250 (MQ=255)
                                                                                                                                                                                        ccT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTGATGCTGAACGGgcagc                                                                    >  2:132453/1‑248 (MQ=255)
                                                                                                                                                                                         cT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGgcagc                                                                    >  1:381457/1‑247 (MQ=255)
                                                                                                                                                                                                aaCATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGACCGGGCAGCAGATTGACAt                                                          >  1:41873/1‑251 (MQ=255)
                                                                                                                                                                                                    aacccagcccaaccggggggATTGCGGGAAACGCATCGGCAAAGGCGTG‑TGTG‑TAATACCCTTCCCCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATAc                                                        <  1:3112/232‑1 (MQ=255)
                                                                                                                                                                                                            aaaC‑GTGG‑AGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCg                                              >  1:402575/1‑251 (MQ=255)
                                                                                                                                                                                                                 gtgg‑tgATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCat                                          >  1:11522/1‑251 (MQ=255)
                                                                                                                                                                                                                   gg‑tgATTTCGGGAAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCctat                                        >  1:128397/1‑251 (MQ=255)
                                                                                                                                                                                                                    g‑tgATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTCCCCCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATc                                       <  1:344060/251‑1 (MQ=255)
                                                                                                                                                                                                                             gggaaacgtacgggaaaaTGGGTG‑TGTG‑AAATACCCTCCCCCCCACCCAAAATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGAt                               <  2:11522/236‑1 (MQ=255)
                                                                                                                                                                                                                                 catcgcagcggaaattgcttgtttg‑tAATCCCCTGCCCCACAGCCAATACCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCGGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTc                            <  1:66737/228‑1 (MQ=255)
                                                                                                                                                                                                                                         ggCAAATGCGTG‑TGTG‑TAATCCCCTTCCCCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGAt                                                                <  1:169730/206‑1 (MQ=255)
                                                                                                                                                                                                                                         ggCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGAt                                                                >  2:169730/1‑206 (MQ=255)
                                                                                                                                                                                                                                          gaaaaTGCGTG‑GGTG‑TAATCCCCTTCACCACAGCCAATATCCACCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACa                  <  1:8420/249‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaTACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGtttttt  >  1:485690/1‑251 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        
GTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCT‑GCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCT‑TCGAAACATCCAGACCAAAC‑GTGG‑TGATTTCGGGCAACGCATCGGCAAATGCGTG‑TGTG‑TAATACCCTTCACCACAGCCAATATCCAGCACC‑GCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTT  >  NC_000913/1906520‑1907016

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: ‑