Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 692,133 G→A intergenic (+194580/‑7210) unknown → / → unknown tRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1692,1330GA100.0% 65.0 / NA 26intergenic (+194580/‑7210)unknown/unknowntRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/14);  total (12/14)

CGGTAGGGCTAAATACAAAGACAGTGGAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGGATATGTGGACCCCTTTCAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAACTAACTGTTAGTTCGGCAGTAAGTAAACAAAGATTATTATAGTTTGTT  >  1/692046‑692359
                                                                                       |                                                                                                                                                                                                                                  
cGGTAGGGCTAAATACAAAGGTAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGgtg                                                                >  1:901509/1‑251 (MQ=255)
  gTAGGGCTAAATACAAAGGTAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTa                                                              >  1:239399/1‑251 (MQ=255)
      ggCTAAATACAAAGGTAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAAttt                                                          >  1:88441/1‑251 (MQ=255)
      ggCTAAATACAAAGGTAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAAt                                                            <  1:173696/249‑1 (MQ=255)
            aTACAAAGGTAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTAcc                                                      <  1:91189/249‑1 (MQ=255)
               caaaggtaGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGaa                                                  >  1:555476/8‑250 (MQ=255)
                       gtgtaggagTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAACTAACTg                                           >  1:654396/5‑249 (MQ=255)
                       gtgtaggagTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAACTAACTGtt                                         >  1:178769/5‑251 (MQ=255)
                       gtgtaggagTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAACTAACTGtt                                         >  1:74541/5‑251 (MQ=255)
                            ggagTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAACTAACTGTTAGtt                                     >  1:308845/1‑250 (MQ=255)
                             gagTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAACTAACTGTTAGtt                                     <  1:830410/249‑1 (MQ=255)
                              agTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAAttt                                                          <  1:224001/227‑1 (MQ=255)
                               gTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAAt                                                            <  1:527362/224‑1 (MQ=255)
                               gTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAACTAACTGTTAGTTCgg                                  >  1:444977/1‑250 (MQ=255)
                                tatataAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAACTAACTGTTAGTTCGGc                                 <  1:315518/250‑1 (MQ=255)
                                 atataAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAACTAACTGTTAGTTCGGCag                               <  1:299360/251‑1 (MQ=255)
                                 atacaaGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAACTAACTGTTAGTTCGGCag                               >  1:800039/5‑251 (MQ=255)
                                        gAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGaactaact                                            <  1:359192/231‑1 (MQ=255)
                                        gAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAACTAACTGTTAGTTCGGCagtaagt                          <  1:127330/249‑1 (MQ=255)
                                                aaaaGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAACTAACTGTTAGTTCGGCAGTAAGTAAACAAAGa                 <  1:5896/250‑1 (MQ=255)
                                                    ggAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAAttt                                                          <  1:525547/205‑1 (MQ=255)
                                                      aaGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTGGTACAGGGTGTAATTTTACCTGAACTAACTGTTAGTTCGGCAGTAAGTAAACAAAGATtattat          <  1:519128/251‑1 (MQ=255)
                                                              aTTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAACTAACTg                                           >  1:171817/1‑210 (MQ=255)
                                                              aTTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAACTAACTGTTAGTTCGGCAGTAAGTAAACAAAGATTATTATAgtttgtt  >  1:631322/1‑251 (MQ=255)
                                                                   ggATAATTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAACTAACTGTTAGTTCg                                   <  1:142111/213‑1 (MQ=255)
                                                                        aTTTATCGAAGAAAGAATATGTGGACCCCTTT‑AAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAACTAACTGTTAGTTCGGCAGTAAGTAAACAAAGatta              <  1:108958/229‑1 (MQ=255)
                                                                                       |                                                                                                                                                                                                                                  
CGGTAGGGCTAAATACAAAGACAGTGGAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGGATATGTGGACCCCTTTCAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAACTAACTGTTAGTTCGGCAGTAAGTAAACAAAGATTATTATAGTTTGTT  >  1/692046‑692359

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