Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	571,925	A→G	intergenic (+208508/‑94)	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	571,925	0	A	G	100.0%	62.7 / NA	22	intergenic (+208508/‑94)	unknown/unknown	tRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/5);  total (17/5)

GTAAGTTGTATTTTTTCTAAGCTATCGACAAAAAAATGTAAAAATCCGCTTGATCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAAGATTTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATAATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTTTGAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCTT  >  1/571686‑572070                                                                                                                                                                                                                                                |                                                                                                                                                  gtaagtTGTATTTTTTCTAAGCTATCGACAAAAAAATGTAAAAATCCGCTTGATCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGAtt                                                                                                                                        >  1:819090/1‑251 (MQ=255)        gTATTTTTTCTAAGCTATCGACAAAAAAATGTAAAAATCCGCTTGATCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGAt                                                                                                                                         >  1:433639/1‑243 (MQ=255)             ttttCTAAGCTATCGACAAAAAAATGTAAAAATCCGCTTGATCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAAc                                                                                                                                    >  1:772768/1‑243 (MQ=255)              tttCTAAGCTATCGACAAAAAAATGTAAAAATCCGCTTGATCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGt                                                                                                                           >  1:259269/1‑251 (MQ=255)                                      gTAAAAATCCGCTTGATCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAg                                                                                                                                           <  1:644177/211‑1 (MQ=255)                                                gCTTGATCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAgg                                                                                                                            >  1:536855/1‑216 (MQ=255)                                                     aTCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGAt                                                                                                                                         >  1:831553/1‑198 (MQ=255)                                                     aTCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGatt                                                                                    >  1:633800/1‑251 (MQ=255)                                                     aTCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGatt                                                                                    >  1:642263/1‑251 (MQ=255)                                                        cTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTGCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTataat                                                                                         <  1:450090/243‑1 (MQ=255)                                                             aTAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTg                                                                                                                          >  1:516909/1‑205 (MQ=255)                                                                   aTCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTGCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTa                                                                      <  1:760200/251‑1 (MQ=255)                                                                                 cATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAg                                                         >  1:760756/1‑250 (MQ=255)                                                                                         cTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTataat                                                                                         <  1:39879/210‑1 (MQ=255)                                                                                               gggCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTTTGAATTGAag                                          >  1:575543/1‑251 (MQ=255)                                                                                                 gCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTataat                                                                                         >  1:708114/1‑202 (MQ=255)                                                                                                                         aaaTATTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAgg                                                                                                                            >  1:835565/1‑143 (MQ=255)                                                                                                                             atTTTCTTGACGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTTTGAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTgg            >  1:735344/1‑251 (MQ=255)                                                                                                                                       cGGCAAATAATATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTTTGAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCtt  >  1:235485/1‑251 (MQ=255)                                                                                                                                            aataatATCTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTataat                                                                                         >  1:341244/1‑159 (MQ=255)                                                                                                                                                  atcttCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTataat                                                                                         <  1:335128/150‑1 (MQ=255)                                                                                                                                                                              tAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTttatttatatta                                                                                                 >  1:370777/1‑117 (MQ=21)                                                                                                                                                                                                                                                |                                                                                                                                                  GTAAGTTGTATTTTTTCTAAGCTATCGACAAAAAAATGTAAAAATCCGCTTGATCCTATCATAGAGATCCCTATAATACGCATCCGCACTGCAAGGGCAGTAACTAAAGGAAAATAAACAAAATATTTTCTTGACGGCAAATAAGATTTTCATCTATAATTTGCGGCCTTGCTTAATGCAAATGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATAATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCGTTTATTTATATTATTTATAATAGATTAATTAAGTCAGCTATTGAGGGATTAAGCTTTTTGAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCTT  >  1/571686‑572070

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: ‑