Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,906,7760AC29.6% 104.1 / 12.3 71G95G (GGT→GGG) rlmA23S rRNA m(1)G745 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (38/12);  new base C (0/21);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 6.48e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGT‑CCATAC‑TGGTATCGGAAAACGG‑CAAACGGTGG‑CTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGA  >  NC_000913/1906551‑1907024
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
gggCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAACTACGTATTATGGCGT‑CCATAC‑TGGTATCGGCCCACGG‑CAAACGGTGG‑CTGGAAGCGCCACAACAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTCTCGCTACCTTCGCAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGGGTACTACCCTTCACCACAGCCAATCTCCAGCAccg                                                                                                                                                                                                                                    >  2:419112/1‑248 (MQ=255)
gggCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACTTATTATGGCGT‑CCATAC‑TGGTATCGGAAAACGG‑CAAACGGTGG‑CTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAAGGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCAccg                                                                                                                                                                                                                                    <  1:419112/248‑1 (MQ=255)
      cACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGT‑CCATAC‑TGGTATCGGAAAACGG‑CAAACGGTGG‑CTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCt                                                                                                                                                                                                                           >  1:242510/1‑251 (MQ=255)
         tactCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGT‑CCATAC‑TGGTATCGGAAAACGG‑CAAACGGTGG‑CTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTat                                                                                                                                                                                                                        >  1:164505/1‑251 (MQ=255)
               tGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGT‑CCATAC‑TGGTATCGGAAAACGG‑CAAACGGTGG‑CTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTAtcatca                                                                                                                                                                                                                   >  1:246920/1‑250 (MQ=255)
                        ttctGCTTTACACGGCGAGTAAATACGTATTATGGCGT‑CCATAC‑TGGTATCGGAAAACGG‑CAAACGGTGG‑CTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTccc                                                                                                                                                                                                         >  1:144470/1‑251 (MQ=255)
                         gttccTTTAAACGCCGCGTAAATACGTATGATGG‑GTACCATCC‑TGGTGCGGGAAAACGG‑CAAATGGGGG‑CTGGAAGCTACACAAAAAGGGACCGCCGGTGACGGTTTCCCCGCCG‑TTTTTATCGCTACTTTCCAAAAATCCAGGCCAAACGGGTTTATTGGGGTGAACT‑CATCGGAAAATGAGTTTGGGTAATCCCCTTCACAACAGCTGATATCCCGCACCGCCGTGGCCTTATCATCAAGCCGTTcc                                                                                                                                                                                                          <  1:275146/245‑1 (MQ=255)
                            gttttACACGGCGCGTAAATACGTATTATGGCGT‑CCATCC‑TGGTATCGGAAATCGG‑CAAACGGTGG‑CTGGAAGCGACACAAAAAGGGACCTGCGGAGAGCTTTTCGCCGCCG‑CTTTTATCGCTACCTTGGAACCACCCAGACCAAACGTGGTGATTTCGGGCAAGG‑CATCGGCAATTGCGTGTGTGTAATCCCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCa                                                                                                                                                                                                      <  1:226011/248‑1 (MQ=255)
                               ttACACGGCGCGTAAATACGTATTATGGCGT‑CCATAC‑TGGTATCGGAAAACGG‑CAAACGGTGG‑CTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTCCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTg                                                                                                                                                                                                  >  1:340385/1‑251 (MQ=255)
                               ttACACGGCGCGTAAATACGTATTATGGCGT‑CCATAC‑GGGTATCGGAAACCGG‑CAAACGGTGG‑CTGGAAGCGACACAAAAAGTGCCCGGCGGATAGCTTTTCGCCGCCG‑CTTTTATCGCTCCCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATCCCCTTCCCCACAGCCAATATCCAGCGCCGCCGTTGCCTTATCAGCAAGCCGTTCCCTCAGtt                                                                                                                                                                                                   <  1:369814/250‑1 (MQ=255)
                               ttACACGGCGCGTAAATACGTATTAGGGCGT‑CCATAC‑TGGTATCGAAAAACGG‑CAACCGGTGG‑CTGAAAGCGACAAAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCAAAACATCCAGCCCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTTTAATACCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGtt                                                                                                                                                                                                   <  1:272685/250‑1 (MQ=255)
                                 acccGGCGCGAAAATACTTATTATGGCGT‑CCATAC‑TGGTACCGAAAAAGGG‑CAAACGGTGG‑CGGGAGGCGACACAAAAAGTGACCTGCGGGTGGCGGTGCGCCGCCG‑CTTTTACCGCTCCCTTCGAAAGATCCAGCCCAACGGTGGTGTGGTCGGGCAACT‑CCTGCGAAAATGCTTCTGTGTAATCCCCTTCCCCACAGCCATTACCCTGCACCGCCGTGGTCTTATCATCAAGCCGTTCCCTCAGTTggg                                                                                                                                                                                                <  1:456754/248‑1 (MQ=255)
                                 acacGGCGCGTAAATACGTATTATGGCGT‑CCATAC‑TGGTATCGGAAAACGG‑CAAACGGTGG‑CTGGAAGCGACACAAAAAGTGCCCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATCCCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTgg                                                                                                                                                                                                 <  1:14532/250‑1 (MQ=255)
                                     ggcgcataaagccgaattatGGCGT‑CCATAC‑GGGTCTCG‑AAAACGGGCAAGCGGGGG‑CGGGAAGCGACAAAAAAAGTACCCGGCGGAAAGCGTTTCGCCCCCG‑CTTTTACCGCTACCTTCGAAACATCCAGCCCAAACGGGGTGTTTTCGGGCACCG‑CATCGGCAAAGGCGTGTGGGTAATCCCCTTCACCCCAGCCAAGATCCAGCACCGCCGAGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGAc                                                                                                                                                                                            <  2:161620/236‑1 (MQ=255)
                                     ggCGCGTAAATGGGTATTAGGGCG‑‑CCATACGTGGTCTCGGAAA‑CGG‑CAAACGGGGGGCTGGAAGCAACAAAAAAAGTAACTTGCGAATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCAAAACAGCCAGACCAAACGTGGTGATTGGGGGCAACG‑CATCGGCAAATGCGTGTGGGTAATACCCTTCCCCACGGCCAATAGCCAGCCCCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGAc                                                                                                                                                                                            <  2:92909/251‑1 (MQ=255)
                                     ggCGCGTAAATACGTATTATGGCGT‑CCATAC‑TGGTATCGGAAAACGG‑CAAACGGTGG‑CTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGAc                                                                                                                                                                                            >  1:382265/1‑251 (MQ=255)
                                       cgcgTAAATACGTATTATGGCGT‑CCATAC‑TGGTATCGGAAAACGG‑CAAACGGTGG‑CTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACa                                                                                                                                                                                           >  1:247239/1‑250 (MQ=255)
                                            aaaTACGTATTATGGCGT‑CCATAC‑TGGTATCGGAAAACGG‑CAAACGGTGG‑CTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGc                                                                                                                                                                                      >  1:455804/1‑250 (MQ=255)
                                                  gTATTAGGGCGC‑CCATAC‑TGGTATCGGAAACCGG‑CAATCGGTGG‑CTGGAGGCGACAAAAAAAGTGCCCTGCGGATAGCTTTCCCCCGCCT‑CTTTTATCGCTCCCTTCGAAACATCCAGCCCAACCGTGGTGATTTCGGGCAACG‑CCTCGGCAAATGCGTGTGTTTAATCCCCTTCACCCCGGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGc                                                                                                                                                                               <  1:115630/251‑1 (MQ=255)
                                                      aatGGCGT‑CCATCC‑TGGTATCGGAAAAGGG‑CAAACGGTGG‑CTGGAAGCAAAACAAAAGGTGCCCTGCGGAAAGGTTTTCGCCCCCC‑CTTTTATCGCTCCCTCCGAAACATCCAGACCAAACGTGGTGATTTCGGGCACCG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGc                                                                                                                                                                            <  1:655/249‑1 (MQ=255)
                                                             t‑CCATAC‑TGTTAACGGAAAGCGG‑AAAACGGTGC‑CTGGAAGCGACACAAAAAGTGCCCGTCGGATGACGTCTTGCCCCCA‑TTTTTATCGCTACCTTCAAAACGCCCAAACCAACCGTGGTTATTCCGGGCACCG‑CACCGCCAAGTGCTTGTTTGTAATCCCCTTCCCCCCGCCCAATACCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGAt                                                                                                                                                                     <  1:64970/250‑1 (MQ=255)
                                                               ccATAC‑TGGTATCGGAAAACGG‑CAAACGGTGG‑CTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATa                                                                                                                                                                    >  2:116192/1‑250 (MQ=255)
                                                                            cGGAAAACGG‑CAAACGGTGG‑CTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCt                                                                                                                                                       >  2:379087/1‑251 (MQ=255)
                                                                            cGGAAAACGA‑CAAACGGGGG‑CAGAAAGTGACAAAAAAAGTAACATGCGGATAGCGTTTCGCCGCTG‑TTTTTACCGCGCCGTCCTAACCAACCAGACCAAGCGGGGGGATTAAGGACATGG‑CACGGGCAAATGCGTTTGTGGAATACCCTTCCCCACAGCCAAAATCCAGCGCCGCCGTGGCCTGATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAAGGTCCGGCATCt                                                                                                                                                       <  1:417604/251‑1 (MQ=255)
                                                                             ggAAAAGGG‑CAAACGGTGG‑CTGTAAGCGACAAAAAAAGTGTCCTGCGGATAGCGTGCCGGC‑CCGTCTTTTTTCGATCCTTTCGAACCCCCCAGACCAATGTTGGTGTTTTGGGGCCACG‑AATCGGCAAATGCGTGTGTGTAAACCCCTCCACCCCAGCCAATATCCAGCACCGCCGTGGCCTTATCAGCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAAGGTCCGGCAtata                                                                                                                                                      <  2:57157/251‑4 (MQ=255)
                                                                             ggAAAACGG‑CAAACGGTGG‑CTGGAAGCGACACAAAAAGTGACCTGCGGATAGAGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCACTATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCt                                                                                                                                                       >  2:251119/1‑250 (MQ=255)
                                                                               aaagcGG‑CAAACGGTGG‑CTGGAGGCGACAAAAAAAGTGCCCGGCGGATAGCGTTGCGCCGCCG‑CTTTTAGCGCTACCTTCAAAACATCCAGCCCAAACGTGGTGATTCCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAACACCCTTCACCACAGCCATTATCCAGAACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAg                                                                                                                                                    <  1:14013/247‑1 (MQ=255)
                                                                                            ggTGG‑CTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCgcggcgg                                                                                                                                        >  2:340490/1‑250 (MQ=255)
                                                                                               gg‑CTGGAAGCGACACACCAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTCTCGCTACCTTCGAAACATCCAGACCAAGCGTGGTGATTTCGGGCAACG‑CATAGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCCAGCCGTTCCCTCCGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCt                                                                                                                                     >  2:467772/1‑251 (MQ=255)
                                                                                                   gataaccGACACAAAAAGGGACCTTCGGAGTGCTTTTCCCCGCCG‑CTTTTCACGCTCCCTTCGAAACACCCAGACAAAACTGGGTGATTTCGGACAACG‑CTTCGGCAAATGCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCCCCGCCTTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCAGCTAATAATGCGCGGCGTGCTTGc                                                                                                                                  <  1:121782/245‑1 (MQ=255)
                                                                                                    ggAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGACG‑CTTTTATCACTACCTTCGAAACATCAAGACCAAACGTGGTGATTCCGGGCAAAG‑CATCGGCAGATGCGTGTGTGTAATACCCTTCACCAAAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGc                                                                                                                                  >  1:454678/1‑250 (MQ=255)
                                                                                                     gAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATC‑AGACCAAACGTGGTGAATTCGGGCAACG‑CATCAGCAAATGCGTGTGTGTAATACCCTTCACCCCAGCCAATATCCCGCACCGCCCGGGCCTTATCATCAAGACGCTCCCTCACTTTGGCGCCAATTTCCTCACGCCGCGGCTGATAATGTCCGGCCTCTAAGAATTCGGGGCGTTCTTGcgac                                                                                                                               >  2:96915/1‑248 (MQ=255)
                                                                                                     agagccacacAAAAAGTCACCGGGGGATGGCGTTGCGCCGCCT‑CTTTGATCGCTCCCTTCAAACCAGCCAGACCAAACGGGGTGATTTCGGGCAACG‑ACTCGGCAAAGGCGGGTGTGAAATCCCCTCCCCCACAGCCAATACCCAGCCCCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGcat                                                                                                                                <  2:344573/245‑1 (MQ=255)
                                                                                                      aaGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGcat                                                                                                                                >  1:344573/1‑250 (MQ=255)
                                                                                                       aGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGcatcc                                                                                                                              >  2:120353/1‑250 (MQ=255)
                                                                                                       aGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGcatc                                                                                                                               >  2:338025/1‑250 (MQ=255)
                                                                                                           ccaaaaaaatggcccTGGGGATAGCTTTCCGCCGCCT‑TTTTTATCGCTCCCTTCAAAACATCCAGCCAAAACGGGGTGTTTGCGGGAAACC‑CTCCGGAAAATGGGGGTGTTAAATCCCCTTCCCCCCAGCCAATACCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTc                                                                                                                          <  2:120838/238‑1 (MQ=255)
                                                                                                             acaAAAAGGCACCTGCGGATGCCTTTTCCCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAACCGTGGTTATTTCGGGCAACG‑CATCGGCAAAGGGGTGTGTGTAATACCCTTCCCCACAGCCAATACCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCcg                                                                                                                        <  1:433748/251‑1 (MQ=255)
                                                                                                               aaaaaGTGCCCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCGTGCATCATTTCcgcg                                                                                                                      <  1:116192/251‑1 (MQ=255)
                                                                                                                    gTGACATGCGGCTAGCGTTTCGCCGCCG‑CTTTTATCGATACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTc                                                                                                                 >  2:422817/1‑251 (MQ=255)
                                                                                                                               aTAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCAc                                                                                                      >  1:424110/1‑251 (MQ=255)
                                                                                                                                  cctgtgccacgcgg‑cTTTTATCGCTGCTTTCGAAACACCCACACCAATCGAGGTTGTTTCGGTCATCC‑CCTCGGCCAATGCGTGTGTTTAATACCCTTCCCCCCAGCCAATACCCTGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACgag                                                                                                   <  1:255675/238‑1 (MQ=255)
                                                                                                                                      ttCGCCCCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGAAAATGCGTGTGTGTAATACCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGAcc                                                                                                <  1:346678/250‑1 (MQ=255)
                                                                                                                                         cccgccg‑cTTTTATCGCTCCCTTCAAAACATCCAGACCAAACGTGGTGTTTTCGGGAAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGttt                                                                                            <  1:57130/250‑1 (MQ=255)
                                                                                                                                           cgccg‑cTTTTACCGCCCCCTACGAAACACCCACCCCAACCGTGGTAATGTCGGCCAACG‑CACCGGCAAATGCTGGTGTGAAATCCCCTCCCCCCCAGCCAATACCCAGCCCCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTGGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGTATCACGAGACCGTTTAt                                                                                          <  2:424110/251‑1 (MQ=255)
                                                                                                                                               g‑cTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTg                                                                                      >  2:20072/1‑251 (MQ=255)
                                                                                                                                               cgcgtttATGGCTACCTTCGAACCATCCAATCCAAGCGTGGTAATTTCGGGAAACG‑CGACGGCAAAGGGGTGTTTGTAAAACCCTGCACCCCAGCCATTTTCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGGGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCt                                                                                       <  1:467772/247‑1 (MQ=255)
                                                                                                                                                 cTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGCCCGTTTATGCTg                                                                                      >  2:113646/1‑250 (MQ=255)
                                                                                                                                                 cTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTg                                                                                      >  1:152242/1‑250 (MQ=255)
                                                                                                                                                    ttATCGCTACCTTCGAAACATCCAGCCCAACCGTGGGGATTTCGGGCAACG‑CATCGGCAAATGGGTGTGTGTAATCCCCTTCCCCACAGCCAATATCCAGCCCCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAAc                                                                                   <  1:190256/250‑1 (MQ=255)
                                                                                                                                                     aatgcgtaccttgaaacctcccAGACCAAACGTGGTGATTT‑GGGGAACGGCAGGGGAAAAGGGGTGGGTGAAATCCCCCTCCCCACAGCAAATACCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACg                                                                                  <  2:382265/230‑1 (MQ=255)
                                                                                                                                                      aTCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CGTAGGCAAATGCGTGTGTGTAATACCCTTCACCACAGACAATATCCAGCACCGCCGTGGCCTTATCGTCAAGCCGTTACCTCAGTTGGGCGCCAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGACTGCGCGGCGTGCTTGCATCATTTCCGCACTGTCGCCCGGATCAAGAGACCGTGTATGCTGAACggg                                                                                >  2:44337/1‑251 (MQ=255)
                                                                                                                                                       ccGCTCCCTTCGAAACACCCAGACCAAAGGTGGTGATTTGGGGCACCG‑CATCGGCAAATGCGTGTGTGTAATCCCCTTCCCCACAGCCAATACCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCATTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTGGCATCATTTCCTCGCGGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGgc                                                                               <  1:323272/250‑1 (MQ=255)
                                                                                                                                                              cTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTCAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGgcagc                                                                            >  1:229228/1‑247 (MQ=255)
                                                                                                                                                                tcgaaccgcccAGACCAACCGTGGTGATTGCGGGAACCG‑CATCGCAAAAGGCGTGTGTGAAATCCCCTCCCCCACGCCCAATACCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGAtt                                                                       <  2:117076/241‑1 (MQ=255)
                                                                                                                                                                    aaCATCCACACCAAACGGGGAGATTCCGGGAAACG‑CAGGGGCAAATGCGTGTGTGTACTCCCCTCCACCCCAGCAAATATCCAGCGCCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACAt                                                                  <  2:144470/251‑1 (MQ=255)
                                                                                                                                                                            gACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCcttct                                                          <  1:257551/251‑1 (MQ=255)
                                                                                                                                                                                          aTTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCaa                                             >  2:242253/1‑250 (MQ=255)
                                                                                                                                                                                           tttCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACACACAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGATGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAAc                                           >  1:226362/1‑251 (MQ=255)
                                                                                                                                                                                            ttCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACt                                          >  1:58311/1‑251 (MQ=255)
                                                                                                                                                                                              cGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTg                                         >  2:46358/1‑250 (MQ=255)
                                                                                                                                                                                               gggCAACG‑CATCGGCAAATGCGTGTTTGTAATACCCTTCACCACAGCCCATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGTGCGACACTTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGCCATACCCGTCTTTCGCCATGTCAAACTGAt                                       >  2:8790/1‑251 (MQ=255)
                                                                                                                                                                                                      g‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTg                                >  2:110756/1‑251 (MQ=255)
                                                                                                                                                                                                           cggaaaaTGCGTGGGTGTAATCCCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTgggg                             <  1:4326/246‑1 (MQ=255)
                                                                                                                                                                                                            ggCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCACGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGAc                           >  1:166588/1‑251 (MQ=255)
                                                                                                                                                                                                            ggCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCACACTGATGTCGCTGGGGAc                           >  2:295203/1‑251 (MQ=255)
                                                                                                                                                                                                              cAAATGCGTGTGTGTAATACCCTTCACCAAAGACAATATCCAGCACCGCAGTGGCCTTCTCATCCAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGAGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGCTTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGAc                           >  2:458707/1‑249 (MQ=255)
                                                                                                                                                                                                                  ggcgtgtgtaaaaaccccTCCACCACAGCAAATACCCAGCCCCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGatat                     <  2:247239/236‑1 (MQ=255)
                                                                                                                                                                                                                             aaTACCCTTCACCACAGCCAATCTCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTggg                                                                                                                                                                                                >  2:45161/1‑69 (MQ=255)
                                                                                                                                                                                                                                    ttCACCACAGCCAATATCCAGAACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCCCGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGa                                        >  2:463966/1‑214 (MQ=255)
                                                                                                                                                                                                                                     tCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATAATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCAcgcg   >  1:210757/1‑250 (MQ=255)
                                                                                                                                                                                                                                      caccacAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCa                                              >  2:66690/1‑206 (MQ=255)
                                                                                                                                                                                                                                      caccacAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGa  <  1:295203/250‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
GGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGT‑CCATAC‑TGGTATCGGAAAACGG‑CAAACGGTGG‑CTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCG‑CTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACG‑CATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGA  >  NC_000913/1906551‑1907024

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