Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	443,388	C→G	intergenic (+79971/‑128631)	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
*	1	443,388	0	C	G	100.0%	90.1 / NA	35	intergenic (+79971/‑128631)	unknown/unknown	tRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (13/22);  total (13/22)

AGGAGTAAGATTTAAAACGCCAGAGCAATCTAGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTACTTATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTCCGTTAAATTTGAGGGTTGTTAGTTGTTGGTGATGTGTAAAATTTACGTAGTGTTACTTAATTTAAAACCATATAAATTTAGAAAATTTATTTCTAGTTTTATATGGTTATTTTTCTAAAATCTCGCGGTCGGCTACACGCAAAACTGCTT  >  1/443141‑443538                                                                                                                                                                                                                                                         |                                                                                                                                                       aagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGt                                                                                                                                                     <  1:442944‑M2/221‑1 (MQ=255)       taaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAAtt                                                                                                                                               <  1:436222‑M2/227‑1 (MQ=255)                 gagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATAC‑CCGTCGTTTTTGTTTTTGTTTATATGATGATTTTGGGTTAAATttgagggcgat                                                                                                                                     >  1:421794‑M2/15‑247 (MQ=255)                      atctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTtaaatttgagggcgattagtc                                                                                                                                <  1:681563‑M2/242‑10 (MQ=255)                      atctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTtaaatttgagggcgattagtc                                                                                                                                <  1:824457‑M2/242‑10 (MQ=255)                      atctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAAt                                                                                                                                                <  1:290052‑M2/226‑1 (MQ=255)                      atctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATAC‑CCCTCGGTTTTGTTTTTGTTTTTATTATGGTTCTGCGGTtaattttagggcgattagt                                                                                                                                 >  1:279353‑M2/10‑242 (MQ=255)                            cggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTctgcgttaaatttgagggcgattagtcgttatt                                                                                                                          <  1:859078‑M2/248‑16 (MQ=255)                            cggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTctgcgttaaatttgagggcgattagtcgttatt                                                                                                                          <  1:321418‑M2/248‑16 (MQ=255)                              gtggCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAgggcgattagtcgttattc                                                                                                                         <  1:603904/247‑17 (MQ=255)                               tggCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAgggcgattagtcgttattca                                                                                                                        <  1:767941/248‑18 (MQ=255)                                    gTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAgggcgattagtcgttatt                                                                                                                          <  1:462008/243‑16 (MQ=255)                                    gTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGAGGGAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAgggcgattagtcgttatt                                                                                                                          >  1:432932/1‑228 (MQ=255)                                              gCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAggg                                                                                                                                         <  1:691838/218‑1 (MQ=255)                                              gCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTTAgggcgattagtcgttatt                                                                                                                          >  1:647138/1‑218 (MQ=255)                                                cTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGt                                                                                                                                                     <  1:13219/204‑1 (MQ=255)                                                cTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTTTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAg                                                                                                      >  1:793873/1‑251 (MQ=255)                                                        cTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                                                                                              <  1:126024/251‑5 (MQ=255)                                                        cTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                                                                                              <  1:622206/251‑5 (MQ=255)                                                             tgtTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttatta                                                                                         <  1:783339/251‑10 (MQ=255)                                                             tgtTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttatta                                                                                         <  1:821071/251‑10 (MQ=255)                                                                   aGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTTTTTATGTTTATATGATGATTCTGCGTTAAATTTTAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAgtgt                                                                                                   >  1:322279/1‑235 (MQ=255)                                                                   aGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTTAgggcgattagtcgttatt                                                                                                                          >  1:348518/1‑197 (MQ=255)                                                                     caaacaaaGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaattattc                                                                                 <  1:543917/251‑18 (MQ=255)                                                                      aaacaaaGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAA‑CCGTCGCTTTTGTTTTTGTTTATATGATGATTCTGCGTTAAATTTGAgggcgatttgtcgttatt                                                                                                                          >  1:649193/1‑194 (MQ=255)                                                                      aaacaaaGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATAC‑CCATCGTTTTTGTTTATGTTTATATGATGATTTTGCGTTAAATTTTAGGGGGGTTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                                                                                              >  1:158898/1‑233 (MQ=255)                                                                         caaaGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaattattcatac                                                                             <  1:859192/251‑22 (MQ=255)                                                                                     gctgAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaatta                                                                                    <  1:482745/232‑15 (MQ=255)                                                                                        gAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTTTGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                                                                                              >  1:826309/1‑215 (MQ=255)                                                                                                           tGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgtt                                                                                               >  1:775368/1‑196 (MQ=255)                                                                                                                                    gCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttatta                                                                                         >  1:440256/1‑171 (MQ=255)                                                                                                                                                       aaaTAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                                                                                              <  1:258670/156‑5 (MQ=255)                                                                                                                                                         aTAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaat                                                                                      <  1:284319/162‑13 (MQ=255)                                                                                                                                                                    gctCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaat                                                                                      <  1:366774/149‑13 (MQ=255)                                                                                                                                                                                                     aaaCGGCGGGTTAATACTCCGCCGTTTTTCTTTTCATGTACATGATGATTCTGCTTTAAATTTGAGGGTAATTAGTCGCTGATGATGTGTAAAATTGACGTAGCGCTATTTGGCGTTAAACCACATAAGTTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACgtccggtaatgt  <  1:82718/203‑12 (MQ=255)                                                                                                                                                                                                                          tcgtTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                                                                                              <  1:72476/90‑5 (MQ=255)                                                                                                                                                                                                                                               tgatgaTTCTGCTTTAAATTTGAGGGTAATTAGTCGCTGATGATGTGTAAAATTGACGTAGCGCTATTTGGCGTTAAACCACATAAGTTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAacac              >  1:770647/1‑149 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                       AGGAGTAAGATTTAAAACGCCAGAGCAATCTAGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTACTTATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTCCGTTAAATTTGAGGGTTGTTAGTTGTTGGTGATGTGTAAAATTTACGTAGTGTTACTTAATTTAAAACCATATAAATTTAGAAAATTTATTTCTAGTTTTATATGGTTATTTTTCTAAAATCTCGCGGTCGGCTACACGCAAAACTGCTT  >  1/443141‑443538

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