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

CGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGC  >  NC_000913/1906527‑1907008
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
cGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTcac                                                                                                                                                                                                                                         <  1:367899/251‑1 (MQ=255)
      aggaaTGCCCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTGCTGCTTTACACGGCGCGTAAATACGTAGTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGc                                                                                                                                                                                                                                   <  2:248437/248‑1 (MQ=255)
                 aGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGc                                                                                                                                                                                                                        >  1:285317/1‑251 (MQ=255)
                  gccgccGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCa                                                                                                                                                                                                                       <  1:426942/251‑1 (MQ=255)
                    cgccgGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCAcc                                                                                                                                                                                                                     >  1:240237/1‑251 (MQ=255)
                       cgGGCTTCACTACTCGTGCTCATTCTTCTGCTTTACACGGCGCGTAAATCCGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGACGCGACACACAAAGTGCCCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCACACGTGGTGCTTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACCTCCAATATCCAGCAccgcc                                                                                                                                                                                                                  >  2:33866/1‑251 (MQ=255)
                       cgGGCTTCACTACTCGTGCAAATGCTTCTGCTTTACACGGCGCGTAAATACTTATTATGGCGCCCAAACTGTTATGGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGGGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTTCCTT‑CGAAACAGCCAGACCAAACGGGGTGATTTCGGGCACCGC‑ATCGGCAAATGCGTGTGTGAAAGACCCTTCACCACAGCCAATAGCCAGCAccgcc                                                                                                                                                                                                                  <  2:48469/251‑1 (MQ=255)
                               actactCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATGCGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAACCGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCTTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTTTGTTTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCtt                                                                                                                                                                                                          <  1:289146/251‑1 (MQ=255)
                                           aaTTCTTCTGCTTTACACGGCGCGTAAATCCGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCg                                                                                                                                                                                              >  1:147070/1‑251 (MQ=255)
                                             tcttcctGCTTTACACGCGGCGAAAATACGTATTAGGGCGTCCAACCGGGTATGGAAAAACGGAAACCGGTGGCGGGAGGCAACAAAAAAAGTACCCTGGGGATAGGGTTGCC‑CCGCCGTTTTTATCGCTACCTT‑CGAAACATCCAGACCAACCGTGGTGTTTGCGGGCAACGC‑TTCGGCAAATGCGTGTGTTTAATCCCCTTCCCCACAGCCAATATCAAGCACCGCCGTGGCCTTATCATCAAGCCGt                                                                                                                                                                                             <  2:128692/245‑1 (MQ=255)
                                                 tctGCTTGACACGGCGCGAAAAGACGTAGTATGGCGTCCACACTGGTATCGAAAAGCGGCAAAGGGTGGCTGGAAGCCACAAAAAAAGTGCCCTGCGGATAGCGTTCCC‑CCCCCGCTTTTCTCCCTACCTC‑CGAAACACCCATACCAACCGTGGGGATTTCGGGCAACGC‑ATGGGCAAAGTCTGGTGGGTAACCCTCTTCCCCGCAGCCAAGATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTc                                                                                                                                                                                           <  2:254391/248‑1 (MQ=255)
                                                       gtccccGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTg                                                                                                                                                                                  >  2:274838/6‑251 (MQ=255)
                                                              gcgcGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGAc                                                                                                                                                                            >  2:377422/1‑250 (MQ=255)
                                                                gcgTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAAt                                                                                                                                                                         >  2:112149/1‑251 (MQ=255)
                                                                gcgTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAAt                                                                                                                                                                         >  1:278737/1‑251 (MQ=255)
                                                                          gTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGc                                                                                                                                                               >  1:58056/1‑251 (MQ=255)
                                                                               atGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATCTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCgg                                                                                                                                                          >  2:407120/1‑251 (MQ=255)
                                                                               agggCGCCCATACGGGTAGCGAAAAACGGCAAACGGTGGCTGAAAGCGACAAAAAAAGTGCCCGGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACAAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAAGGCGTGGGTGTAATCCCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCg                                                                                                                                                           <  2:418635/248‑1 (MQ=255)
                                                                                          actggttggggAAAACGCCAACGGGGGGCGGGAAGCAACAAAAAAAGTGACTGGCGGATGGCGTTCCG‑CCGCCTTTTTTACCGCTACCTT‑GGAAACACCCAGACCAAACGGGGTGATTCCGGGCAACGC‑ACCGGCAAATGCGTGTGTTTAATCCCCTCCCCCACACCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTc                                                                                                                                               <  2:138145/242‑1 (MQ=255)
                                                                                             ggTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCgg                                                                                                                                            >  2:144101/1‑251 (MQ=255)
                                                                                                aTCGGAAAACGGCAAACGGTGGCGGGAAGCGACACAAAAAGTTACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTCCCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTTTAATCCCCTTCACCACAGCCAATACCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCAt                                                                                                                                         <  1:285982/251‑1 (MQ=255)
                                                                                                       aaCGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGa                                                                                                                                   >  1:38595/1‑250 (MQ=255)
                                                                                                        acggaaaccGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGC‑TTTCGGCCGCCGCTTTTATCGCTACCTT‑CGAAACACCCAGACCAAACGTGTTGATTTCGGGAAACGC‑TTCGGAAAATGCGTGGGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGTCATCTAAGAAt                                                                                                                                 <  1:370554/243‑1 (MQ=255)
                                                                                                        aCGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGCGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATACTGTCCGGCATCTAAGAAt                                                                                                                                 >  2:64360/1‑251 (MQ=255)
                                                                                                                 ggTGGCTGGAAGCGACACACAAAGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGCAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTACGAATGCGCGGCGt                                                                                                                        >  2:206078/1‑251 (MQ=255)
                                                                                                                  gTGGCTGGAAGCGCCCAAAAAAGGGACCTGCGGATAGCTTTCCG‑CCGCCGCTTTTATCGCTCCCTT‑AGAAACACCCAGACCAAACGGGGTGATTTCGGGCAGCGC‑ATCGGCAAATGCGTGTGTGTAATCCCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTg                                                                                                                       <  1:1459/251‑1 (MQ=255)
                                                                                                                      cgggAAGCGACACAAAAAGTGCCCTGGGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGGGGTGATTTCGGTAAACGA‑TTGGGCAAATGCGTGTGTGTAATACCTTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTg                                                                                                                   <  1:33866/249‑1 (MQ=255)
                                                                                                                         aaaGCGACACAAAAAGTGCCCTGCGGAGAGCTTTCCG‑CCGCCGTTTTTATCGCTACCTT‑CGAAACACCCAGCCCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAAGGCGTGTGTGTAATCCCCTTCCCCACAGCCAATACCCAGCCCCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGcat                                                                                                                <  2:101974/250‑1 (MQ=255)
                                                                                                                           aGCGACACACACAGTGACCTGCGGATAGAGTTTCG‑CCACCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGcatct                                                                                                              >  1:159961/1‑250 (MQ=255)
                                                                                                                              gACACACAAAGTGACCTGCGGATAGCGTTTCG‑ACGCCGCTTTTATCGCTACCTT‑CGAACCATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCAttt                                                                                                           >  2:48267/1‑251 (MQ=255)
                                                                                                                              aacacAAAAAGTGACCTGCGGAGAGCGTTTCG‑CCGCCGCTTTTATCTCTCCCTC‑CAAAACAACCAGACCAACCGGGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGGGTGTGTAAACCCCTTCACCACAGCAAATATCCACCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCAttt                                                                                                           <  1:48267/250‑1 (MQ=255)
                                                                                                                                caaaaaaaGTGACCTGCGGATAGCGTTTCG‑CCGCCTCTTTTATCGCTACCTT‑CGAAGCATCCATACCAACCGGGTTGATTTCGGGCAACGC‑ATCGGCAAATGCTTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTcc                                                                                                         <  1:59038/248‑1 (MQ=255)
                                                                                                                                     aaaggacccTGGGGATAGGGTTCCG‑CCGCCGATTTTACCGCTCCCTT‑CAAAACATCAAGACCAAACGGGGTGTTTTCGGGCAACGC‑ACGGGCAAAGGCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCt                                                                                                    <  2:278737/244‑1 (MQ=255)
                                                                                                                                       aGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAACCATCCAGACCACACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATGCCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCCAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGt                                                                                                  >  1:202385/1‑251 (MQ=255)
                                                                                                                                        gggACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACACCCAGACCAAACGTGGGGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATCCCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTc                                                                                                 <  1:251465/249‑1 (MQ=255)
                                                                                                                                          gACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCCCCATT‑CGAAACCTCCAGACCCAACGTGGTGATTTCGGGCCACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCAAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACACTTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGCATGCGCGGCGTGCTTGACTCATTTCCGCGCTGTCGc                                                                                               >  2:278907/1‑251 (MQ=255)
                                                                                                                                           cccTGCGGATAGCGTGTCC‑CCCCCGTTTTTATCGCT‑CCTTGCAAAACATCCAGACCAAACGTGGTGATTGCGGGCAACGC‑ATCGGCAAAGGCGTGTGTGTAATACCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGcc                                                                                              <  2:388971/250‑1 (MQ=255)
                                                                                                                                             cTGCGGATAGCTTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGCCCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATgcgc                                                                                                                             <  1:378456/218‑1 (MQ=255)
                                                                                                                                             cTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATgcgc                                                                                                                             >  2:378456/1‑218 (MQ=255)
                                                                                                                                                   aTAGCGTTCCG‑CCGCCGCTTTTACCGCTACCTT‑CGAAACACCCAGACCAATCGTGGTGTTTTCGGGCAACGC‑ATCGGAAAAGGCGTGTTTTTAATCCCCTTCACCACAGCCAATACCCAGCTCCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCTCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCAc                                                                                      <  1:455584/251‑1 (MQ=255)
                                                                                                                                                    tAGCGTTTCG‑CCGCCGATTTTATCGCTACCTT‑CGAAACATACAACCCCAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGCGTGTGTAATACCCTTCACCACAGCGAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCACTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGTTAATGTCAGGCATCTAAGACTGCGCGGCGTGCTTGCATAATTTCCGCGCTGTGGCCCGGATCACg                                                                                     >  1:467854/1‑251 (MQ=255)
                                                                                                                                                    aggcctttcccccgccgcTTTTCTCGCTCCCTG‑CGAAACGTCCAGACCAACCGTGGTGTTTCCGGGCAACGC‑ACGGGCAAAGGCGTGGTTGTAATTCCTTTCCCCCCGGCCAAGATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTGCCGCGCTGTCGCCCGGATCAc                                                                                      <  2:285317/240‑1 (MQ=255)
                                                                                                                                                        gTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGTCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACATTTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACTAGAc                                                                                 >  1:14658/1‑251 (MQ=255)
                                                                                                                                                             g‑ccgccgcTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGttt                                                                            >  1:326333/1‑251 (MQ=255)
                                                                                                                                                                  ccgcTTTTATCGATACCTT‑CGAAACATCCAGACCAACCGTACTGATTTCGGGCAACGC‑CTCGGCAAATGCGTGTGTGTAATACCCTTCACCATAGCGAATATCCAGAACTGCCGTGGAATTATAATAAAGCCGTTCCCTCAGTTGGGCGACAGTTGCATCACGCAGCGGCGGATAATGTCCGGCAACTAAGAATGCGCGGCGAGCTTGCCTAATTTCCGCGCTGTCGACCGGCTCACGAGACCGTTTAt                                                                          >  2:465777/1‑249 (MQ=255)
                                                                                                                                                                        ttATCGCTACCTT‑CGAAACATCCAGACCAACCGGGGTGATTTCGGGCAACGC‑ATCGGAAAAGGCGGGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAAc                                                                   <  1:482776/250‑1 (MQ=255)
                                                                                                                                                                        ttATCGCTACCTT‑CAAAACATCCAGACCAACCGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAAc                                                                   <  1:32814/250‑1 (MQ=255)
                                                                                                                                                                           ccGCTACTTT‑CGAAACATCCAGACCAAAGGTGGTGATTTCGGGAAACTC‑ACCGGCAAACGCGGGTGTGTAATCCCCTTCACCACAGCCAATATCCACCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGgc                                                               <  1:194484/250‑1 (MQ=255)
                                                                                                                                                                              gcccccttagaacatacaagcccAATCGTGGTG‑TTTTGGGCAACCCCATCGGAAAATGGGTTGGTGTAAACCCCTTCCCCACAGCCAATTTCCGGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGgcag                                                             <  2:8485/230‑1 (MQ=255)
                                                                                                                                                                                       gAAACATCCAGACCAAACGTGGTGCTTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCAc                                                                                                                                                                 >  2:306087/1‑142 (MQ=255)
                                                                                                                                                                                       gAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCAc                                                                                                                                                                 <  1:306087/142‑1 (MQ=255)
                                                                                                                                                                                                aGACCAAACTTGGTGATTTCGGGCAACGC‑ATCGACAAATGCGTGTGTGTAATCCCCCTCCCCCCAGCCAACATCCAGCCCCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCcttc                                           <  1:913/251‑1 (MQ=255)
                                                                                                                                                                                                    cAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATAACCTTCTTTc                                       >  2:18763/1‑251 (MQ=255)
                                                                                                                                                                                                    cAAACGAGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGAGACAATTGCATCACGCAGCGGCTGATACTGTCCGGCATCTTAGAATGGGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGACTCACGAGCCCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTGCTTTc                                       >  2:18757/1‑251 (MQ=255)
                                                                                                                                                                                                     aaaCGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCg                                      >  1:134273/1‑251 (MQ=255)
                                                                                                                                                                                                              gATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCatataa                              >  2:335003/1‑248 (MQ=255)
                                                                                                                                                                                                                        aaccacctCGGCAATTGCGGGTGTGAAATCCCCCCCCCCACAGCCAAGCTCCAGCCCCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTGCGCCATATCAAACTGATGTc                    <  2:240237/244‑1 (MQ=255)
                                                                                                                                                                                                                                cGGCAAATGCGGGTGTGTAATCCCCTCCACCACAGCCAATATCCAGCCCCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGa            <  1:340184/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aatgcggggtggtAATCCCCTTCCCCCCAGCCAATACCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTGCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGat       <  1:486883/240‑1 (MQ=255)
                                                                                                                                                                                                                                     aaTGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGAc           >  1:6003/1‑247 (MQ=255)
                                                                                                                                                                                                                                      aTGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGat       >  1:149579/1‑250 (MQ=255)
                                                                                                                                                                                                                                         cGTGTGTGTAATACACTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGatat     >  1:308253/1‑249 (MQ=255)
                                                                                                                                                                                                                                          gtgtgtgtAAGCCCCTCCCCCACAGCCAATACCCAGCCCCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTGGCATCATTTCCGCGCTGTCGCCCGGAGCACGAGACCGTTTATGCTGAACGGGCAGCAGAGTGACATACCCTTCTGTCGCCATATCAAACTGATGTCGCTGGGGACAGAGATAGc  <  2:202385/251‑1 (MQ=255)
                                                                                                                                                                                                                                               tgtAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCttt                                        >  1:159177/1‑208 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
CGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCG‑CCGCCGCTTTTATCGCTACCTT‑CGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGC‑ATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGC  >  NC_000913/1906527‑1907008

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