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→GGGrlmA23S 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.

GGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGA  >  NC_000913/1906551‑1907024
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
gggCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAACTACGTATTATGGCGTCCATACTGGTATCGGCCCACGGCAAACGGTGGCTGGAAGCGCCACAACAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTCTCGCTACCTTCGCAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGGGTACTACCCTTCACCACAGCCAATCTCCAGCAccg                                                                                                                                                                                                                                    >  2:419112/1‑248 (MQ=255)
gggCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACTTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAGGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCAccg                                                                                                                                                                                                                                    <  1:419112/248‑1 (MQ=255)
      cACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCt                                                                                                                                                                                                                           >  1:242510/1‑251 (MQ=255)
         tactCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTat                                                                                                                                                                                                                        >  1:164505/1‑251 (MQ=255)
               tGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTAtcatca                                                                                                                                                                                                                   >  1:246920/1‑250 (MQ=255)
                        ttctGCTTTACACGGCGAGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTccc                                                                                                                                                                                                         >  1:144470/1‑251 (MQ=255)
                         gttccTTTAAACGCCGCGTAAATACGTATGATGGGTACCATCCTGGTGCGGGAAAACGGCAAATGGGGGCTGGAAGCTACACAAAAAGGGACCGCCGGTGACGGTTTCCCCGCCGTTTTTATCGCTACTTTCCAAAAATCCAGGCCAAACGGGTTTATTGGGGTGAACTCATCGGAAAATGAGTTTGGGTAATCCCCTTCACAACAGCTGATATCCCGCACCGCCGTGGCCTTATCATCAAGCCGTTcc                                                                                                                                                                                                          <  1:275146/245‑1 (MQ=255)
                            gttttACACGGCGCGTAAATACGTATTATGGCGTCCATCCTGGTATCGGAAATCGGCAAACGGTGGCTGGAAGCGACACAAAAAGGGACCTGCGGAGAGCTTTTCGCCGCCGCTTTTATCGCTACCTTGGAACCACCCAGACCAAACGTGGTGATTTCGGGCAAGGCATCGGCAATTGCGTGTGTGTAATCCCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCa                                                                                                                                                                                                      <  1:226011/248‑1 (MQ=255)
                               ttACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTCCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTg                                                                                                                                                                                                  >  1:340385/1‑251 (MQ=255)
                               ttACACGGCGCGTAAATACGTATTATGGCGTCCATACGGGTATCGGAAACCGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGCCCGGCGGATAGCTTTTCGCCGCCGCTTTTATCGCTCCCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATCCCCTTCCCCACAGCCAATATCCAGCGCCGCCGTTGCCTTATCAGCAAGCCGTTCCCTCAGtt                                                                                                                                                                                                   <  1:369814/250‑1 (MQ=255)
                               ttACACGGCGCGTAAATACGTATTAGGGCGTCCATACTGGTATCGAAAAACGGCAACCGGTGGCTGAAAGCGACAAAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCAAAACATCCAGCCCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTTTAATACCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGtt                                                                                                                                                                                                   <  1:272685/250‑1 (MQ=255)
                                 acccGGCGCGAAAATACTTATTATGGCGTCCATACTGGTACCGAAAAAGGGCAAACGGTGGCGGGAGGCGACACAAAAAGTGACCTGCGGGTGGCGGTGCGCCGCCGCTTTTACCGCTCCCTTCGAAAGATCCAGCCCAACGGTGGTGTGGTCGGGCAACTCCTGCGAAAATGCTTCTGTGTAATCCCCTTCCCCACAGCCATTACCCTGCACCGCCGTGGTCTTATCATCAAGCCGTTCCCTCAGTTggg                                                                                                                                                                                                <  1:456754/248‑1 (MQ=255)
                                 acacGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGCCCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATCCCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTgg                                                                                                                                                                                                 <  1:14532/250‑1 (MQ=255)
                                     ggcgcataaagccgaattatGGCGTCCATACGGGTCTCGAAAACGGGCAAGCGGGGGCGGGAAGCGACAAAAAAAGTACCCGGCGGAAAGCGTTTCGCCCCCGCTTTTACCGCTACCTTCGAAACATCCAGCCCAAACGGGGTGTTTTCGGGCACCGCATCGGCAAAGGCGTGTGGGTAATCCCCTTCACCCCAGCCAAGATCCAGCACCGCCGAGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGAc                                                                                                                                                                                            <  2:161620/236‑1 (MQ=255)
                                     ggCGCGTAAATGGGTATTAGGGCGCCATACGTGGTCTCGGAAACGGCAAACGGGGGGCTGGAAGCAACAAAAAAAGTAACTTGCGAATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCAAAACAGCCAGACCAAACGTGGTGATTGGGGGCAACGCATCGGCAAATGCGTGTGGGTAATACCCTTCCCCACGGCCAATAGCCAGCCCCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGAc                                                                                                                                                                                            <  2:92909/251‑1 (MQ=255)
                                     ggCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGAc                                                                                                                                                                                            >  1:382265/1‑251 (MQ=255)
                                       cgcgTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACa                                                                                                                                                                                           >  1:247239/1‑250 (MQ=255)
                                            aaaTACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGc                                                                                                                                                                                      >  1:455804/1‑250 (MQ=255)
                                                  gTATTAGGGCGCCCATACTGGTATCGGAAACCGGCAATCGGTGGCTGGAGGCGACAAAAAAAGTGCCCTGCGGATAGCTTTCCCCCGCCTCTTTTATCGCTCCCTTCGAAACATCCAGCCCAACCGTGGTGATTTCGGGCAACGCCTCGGCAAATGCGTGTGTTTAATCCCCTTCACCCCGGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGc                                                                                                                                                                               <  1:115630/251‑1 (MQ=255)
                                                      aatGGCGTCCATCCTGGTATCGGAAAAGGGCAAACGGTGGCTGGAAGCAAAACAAAAGGTGCCCTGCGGAAAGGTTTTCGCCCCCCCTTTTATCGCTCCCTCCGAAACATCCAGACCAAACGTGGTGATTTCGGGCACCGCATCGGCAAATGCGTGTGTGTAATACCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGc                                                                                                                                                                            <  1:655/249‑1 (MQ=255)
                                                             tCCATACTGTTAACGGAAAGCGGAAAACGGTGCCTGGAAGCGACACAAAAAGTGCCCGTCGGATGACGTCTTGCCCCCATTTTTATCGCTACCTTCAAAACGCCCAAACCAACCGTGGTTATTCCGGGCACCGCACCGCCAAGTGCTTGTTTGTAATCCCCTTCCCCCCGCCCAATACCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGAt                                                                                                                                                                     <  1:64970/250‑1 (MQ=255)
                                                               ccATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATa                                                                                                                                                                    >  2:116192/1‑250 (MQ=255)
                                                                            cGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCt                                                                                                                                                       >  2:379087/1‑251 (MQ=255)
                                                                            cGGAAAACGACAAACGGGGGCAGAAAGTGACAAAAAAAGTAACATGCGGATAGCGTTTCGCCGCTGTTTTTACCGCGCCGTCCTAACCAACCAGACCAAGCGGGGGGATTAAGGACATGGCACGGGCAAATGCGTTTGTGGAATACCCTTCCCCACAGCCAAAATCCAGCGCCGCCGTGGCCTGATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAAGGTCCGGCATCt                                                                                                                                                       <  1:417604/251‑1 (MQ=255)
                                                                             ggAAAAGGGCAAACGGTGGCTGTAAGCGACAAAAAAAGTGTCCTGCGGATAGCGTGCCGGCCCGTCTTTTTTCGATCCTTTCGAACCCCCCAGACCAATGTTGGTGTTTTGGGGCCACGAATCGGCAAATGCGTGTGTGTAAACCCCTCCACCCCAGCCAATATCCAGCACCGCCGTGGCCTTATCAGCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAAGGTCCGGCAtata                                                                                                                                                      <  2:57157/251‑4 (MQ=255)
                                                                             ggAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGAGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCACTATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCt                                                                                                                                                       >  2:251119/1‑250 (MQ=255)
                                                                               aaagcGGCAAACGGTGGCTGGAGGCGACAAAAAAAGTGCCCGGCGGATAGCGTTGCGCCGCCGCTTTTAGCGCTACCTTCAAAACATCCAGCCCAAACGTGGTGATTCCGGGCAACGCATCGGCAAATGCGTGTGTGTAACACCCTTCACCACAGCCATTATCCAGAACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAg                                                                                                                                                    <  1:14013/247‑1 (MQ=255)
                                                                                            ggTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCgcggcgg                                                                                                                                        >  2:340490/1‑250 (MQ=255)
                                                                                               ggCTGGAAGCGACACACCAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTCTCGCTACCTTCGAAACATCCAGACCAAGCGTGGTGATTTCGGGCAACGCATAGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCCAGCCGTTCCCTCCGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCt                                                                                                                                     >  2:467772/1‑251 (MQ=255)
                                                                                                   gataaccGACACAAAAAGGGACCTTCGGAGTGCTTTTCCCCGCCGCTTTTCACGCTCCCTTCGAAACACCCAGACAAAACTGGGTGATTTCGGACAACGCTTCGGCAAATGCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCCCCGCCTTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCAGCTAATAATGCGCGGCGTGCTTGc                                                                                                                                  <  1:121782/245‑1 (MQ=255)
                                                                                                    ggAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGACGCTTTTATCACTACCTTCGAAACATCAAGACCAAACGTGGTGATTCCGGGCAAAGCATCGGCAGATGCGTGTGTGTAATACCCTTCACCAAAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGc                                                                                                                                  >  1:454678/1‑250 (MQ=255)
                                                                                                     gAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCAGACCAAACGTGGTGAATTCGGGCAACGCATCAGCAAATGCGTGTGTGTAATACCCTTCACCCCAGCCAATATCCCGCACCGCCCGGGCCTTATCATCAAGACGCTCCCTCACTTTGGCGCCAATTTCCTCACGCCGCGGCTGATAATGTCCGGCCTCTAAGAATTCGGGGCGTTCTTGcgac                                                                                                                               >  2:96915/1‑248 (MQ=255)
                                                                                                     agagccacacAAAAAGTCACCGGGGGATGGCGTTGCGCCGCCTCTTTGATCGCTCCCTTCAAACCAGCCAGACCAAACGGGGTGATTTCGGGCAACGACTCGGCAAAGGCGGGTGTGAAATCCCCTCCCCCACAGCCAATACCCAGCCCCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGcat                                                                                                                                <  2:344573/245‑1 (MQ=255)
                                                                                                      aaGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGcat                                                                                                                                >  1:344573/1‑250 (MQ=255)
                                                                                                       aGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGcatcc                                                                                                                              >  2:120353/1‑250 (MQ=255)
                                                                                                       aGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGcatc                                                                                                                               >  2:338025/1‑250 (MQ=255)
                                                                                                           ccaaaaaaatggcccTGGGGATAGCTTTCCGCCGCCTTTTTTATCGCTCCCTTCAAAACATCCAGCCAAAACGGGGTGTTTGCGGGAAACCCTCCGGAAAATGGGGGTGTTAAATCCCCTTCCCCCCAGCCAATACCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTc                                                                                                                          <  2:120838/238‑1 (MQ=255)
                                                                                                             acaAAAAGGCACCTGCGGATGCCTTTTCCCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAACCGTGGTTATTTCGGGCAACGCATCGGCAAAGGGGTGTGTGTAATACCCTTCCCCACAGCCAATACCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCcg                                                                                                                        <  1:433748/251‑1 (MQ=255)
                                                                                                               aaaaaGTGCCCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCGTGCATCATTTCcgcg                                                                                                                      <  1:116192/251‑1 (MQ=255)
                                                                                                                    gTGACATGCGGCTAGCGTTTCGCCGCCGCTTTTATCGATACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTc                                                                                                                 >  2:422817/1‑251 (MQ=255)
                                                                                                                               aTAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCAc                                                                                                      >  1:424110/1‑251 (MQ=255)
                                                                                                                                  cctgtgccacgcggcTTTTATCGCTGCTTTCGAAACACCCACACCAATCGAGGTTGTTTCGGTCATCCCCTCGGCCAATGCGTGTGTTTAATACCCTTCCCCCCAGCCAATACCCTGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACgag                                                                                                   <  1:255675/238‑1 (MQ=255)
                                                                                                                                      ttCGCCCCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGAAAATGCGTGTGTGTAATACCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGAcc                                                                                                <  1:346678/250‑1 (MQ=255)
                                                                                                                                         cccgccgcTTTTATCGCTCCCTTCAAAACATCCAGACCAAACGTGGTGTTTTCGGGAAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGttt                                                                                            <  1:57130/250‑1 (MQ=255)
                                                                                                                                           cgccgcTTTTACCGCCCCCTACGAAACACCCACCCCAACCGTGGTAATGTCGGCCAACGCACCGGCAAATGCTGGTGTGAAATCCCCTCCCCCCCAGCCAATACCCAGCCCCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTGGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGTATCACGAGACCGTTTAt                                                                                          <  2:424110/251‑1 (MQ=255)
                                                                                                                                               gcTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTg                                                                                      >  2:20072/1‑251 (MQ=255)
                                                                                                                                               cgcgtttATGGCTACCTTCGAACCATCCAATCCAAGCGTGGTAATTTCGGGAAACGCGACGGCAAAGGGGTGTTTGTAAAACCCTGCACCCCAGCCATTTTCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGGGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCt                                                                                       <  1:467772/247‑1 (MQ=255)
                                                                                                                                                 cTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGCCCGTTTATGCTg                                                                                      >  2:113646/1‑250 (MQ=255)
                                                                                                                                                 cTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTg                                                                                      >  1:152242/1‑250 (MQ=255)
                                                                                                                                                    ttATCGCTACCTTCGAAACATCCAGCCCAACCGTGGGGATTTCGGGCAACGCATCGGCAAATGGGTGTGTGTAATCCCCTTCCCCACAGCCAATATCCAGCCCCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAAc                                                                                   <  1:190256/250‑1 (MQ=255)
                                                                                                                                                     aatgcgtaccttgaaacctcccAGACCAAACGTGGTGATTTGGGGAACGGCAGGGGAAAAGGGGTGGGTGAAATCCCCCTCCCCACAGCAAATACCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACg                                                                                  <  2:382265/230‑1 (MQ=255)
                                                                                                                                                      aTCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCGTAGGCAAATGCGTGTGTGTAATACCCTTCACCACAGACAATATCCAGCACCGCCGTGGCCTTATCGTCAAGCCGTTACCTCAGTTGGGCGCCAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGACTGCGCGGCGTGCTTGCATCATTTCCGCACTGTCGCCCGGATCAAGAGACCGTGTATGCTGAACggg                                                                                >  2:44337/1‑251 (MQ=255)
                                                                                                                                                       ccGCTCCCTTCGAAACACCCAGACCAAAGGTGGTGATTTGGGGCACCGCATCGGCAAATGCGTGTGTGTAATCCCCTTCCCCACAGCCAATACCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCATTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTGGCATCATTTCCTCGCGGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGgc                                                                               <  1:323272/250‑1 (MQ=255)
                                                                                                                                                              cTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTCAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGgcagc                                                                            >  1:229228/1‑247 (MQ=255)
                                                                                                                                                                tcgaaccgcccAGACCAACCGTGGTGATTGCGGGAACCGCATCGCAAAAGGCGTGTGTGAAATCCCCTCCCCCACGCCCAATACCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGAtt                                                                       <  2:117076/241‑1 (MQ=255)
                                                                                                                                                                    aaCATCCACACCAAACGGGGAGATTCCGGGAAACGCAGGGGCAAATGCGTGTGTGTACTCCCCTCCACCCCAGCAAATATCCAGCGCCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACAt                                                                  <  2:144470/251‑1 (MQ=255)
                                                                                                                                                                            gACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCcttct                                                          <  1:257551/251‑1 (MQ=255)
                                                                                                                                                                                          aTTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCaa                                             >  2:242253/1‑250 (MQ=255)
                                                                                                                                                                                           tttCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACACACAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGATGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAAc                                           >  1:226362/1‑251 (MQ=255)
                                                                                                                                                                                            ttCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACt                                          >  1:58311/1‑251 (MQ=255)
                                                                                                                                                                                              cGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTg                                         >  2:46358/1‑250 (MQ=255)
                                                                                                                                                                                               gggCAACGCATCGGCAAATGCGTGTTTGTAATACCCTTCACCACAGCCCATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGTGCGACACTTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGCCATACCCGTCTTTCGCCATGTCAAACTGAt                                       >  2:8790/1‑251 (MQ=255)
                                                                                                                                                                                                      gCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTg                                >  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)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
GGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGA  >  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: