Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	493,370	G→T	intergenic (+129953/‑78649)	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	493,370	0	G	T	100.0%	105.9 / NA	36	intergenic (+129953/‑78649)	unknown/unknown	tRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (14/22);  total (14/22)

TATTTCATCTTGCAATACAATTCTTTGAATAGGATATGAATTTTTCAATAAACTAATTGAATTTACCTCAAAATCTCTTGCTAAATATATTTCTGTTAGAGCTGCTTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATGTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAACTTCTTTAGCAAATTCATCTTTATTAAAAGTTTTTAATTGGTTTATATTGATCACATTTAAAGTATCTTCAGCCTCTTTAATATTTCTACTATCCATAGAGCTATAAAAAACGGTTATATTCCCAAAGCCCTGATAATGGCGTACCTTATGCCACAACTCATT  >  1/493204‑493593                                                                                                                                                                       |                                                                                                                                                                                                                                tatttCATCTTGCAATACAATTCTTTGAATAGGATATGAATTTTTCAATAAACTAATTGAATTTACCTCAAAATCTCTTGCTAAATATATTTCTGTTAGAGCTGCTTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTta                                                                                                                                                                                  <  1:756252/214‑1 (MQ=255)                      tCTTTGAATAGGATATGAATTTTTCAATAAACTAATTGAATTTACCTCAAAATCTCTTGCTAAATATATTTCTGTTAGAGCTGCTTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGAc                                                                                                                                                                                                                           >  1:490985/1‑152 (MQ=255)                                ggATATGAATTTTTCAATAAACTAATTGAATTTACCTCAAAATCTCTTGCTAAATATATTTCTGTTAGAGCTGCTTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGAc                                                                                                                                                                                                                           >  1:329351/1‑142 (MQ=255)                                ggATATGAATTTTTCAATAAACTAATTGAATTTACCTCAAAATCTCTTGCTAAATATATTTCTGTTAGAGCTGCTTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCt                                                                                                                                                                                           <  1:211705/174‑1 (MQ=255)                                 gATATGAATTTTTCAATAAACTAATTGAATTTACCTCAAAATCTCTTGCTAAATATATTTCTGTTAGAGCTGCTTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGa                                                                                                                                                                                                            >  1:328534/1‑156 (MQ=255)                                             tCAATAAACTAATTGAATTTACCTCAAAATCTCTTGCTAAATATATTTCTGTTAGAGCTGCTTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATg                                                                                                                                                                                                 <  1:776770/155‑1 (MQ=255)                                                       aaTTGAATTTACCTCAAAATCTCTTGCTAAATATATTTCTGTTAGAGCTGCTTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAAt                                                                                                                                                                                             <  1:539977/149‑1 (MQ=255)                                                                  ccTCAAAATCTCTTGCTAAATATATTTCTGTTAGAGCTGCTTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCt                                                                                                                                                                                           >  1:109332/1‑140 (MQ=255)                                                                  ccTCAAAATCTCTTGCTAAATATATTTCTGTTAGAGCTGCTTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAACTTCTTTAGCAAATTCATCTTTATTAAAAGTTTTTAATTGGTTTATATTGATCACATTTAAAGTATCTTCAGCCTCTTTAat                                                                                   <  1:10893/244‑1 (MQ=255)                                                                                gCTAAATATATTTCTGTTAGAGCTGCTTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAg                                                                                                                                                                                         <  1:321351/128‑1 (MQ=255)                                                                                gCTAAATATATTTCTGTTAGAGCTGCTTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAACTTCTTTAGCAAATTCATCTTTATTAAAAGTTTTTAATTGGTTTATATTGATCACATTTAAAGTATCTTCAGCCTCTTTAATACTTCTACTATCCATagag                                                                >  1:389903/1‑249 (MQ=255)                                                                                   aaaTATATTTCTGTTAGAGCTGCTTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGc                                                                                                                                                                                        <  1:103608/126‑1 (MQ=255)                                                                                   aaaTATATTTCTGTTAGAGCTGCTTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAAc                                                                                                                                                                    <  1:683589/146‑1 (MQ=255)                                                                                     atatatTTCTGTTAGAGCTGCTTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAACTTCTTTAGCAAATTCATCTTTATTAAAAGTTTTTAATTGGTTTatat                                                                                                                     <  1:725719/191‑1 (MQ=255)                                                                                         atTTCTGTTAGAGCTGCTTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAACTTCTTTAGCAAATTCATCTTTATTAAAAGTTTTTAATTGGTTTATATTGATCACATTTAAAGTATCTTCAGCCTCTTTAATACTTctacta                                                                         <  1:832168/231‑1 (MQ=255)                                                                                             cTGTTAGAGCTGCTTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAACTTCTTTAGCAAATTCATCTTTATTAAAAGTTTTTAATTGGTTTATATTGATCACATTTAAAg                                                                                                      >  1:739141/1‑198 (MQ=255)                                                                                                        gctTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAggg                                                                                                                                                                                                             <  1:819355/84‑1 (MQ=255)                                                                                                        gctTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGa                                                                                                                                                                                                            >  1:407439/1‑85 (MQ=255)                                                                                                        gctTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCt                                                                                                                                                                                           <  1:155938/102‑1 (MQ=255)                                                                                                        gctTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGc                                                                                                                                                                                        <  1:505543/105‑1 (MQ=255)                                                                                                        gctTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAACTTCTTTAGCAAATTCATCTTTATTAAAAGTTTTTAATTGGTTTatat                                                                                                                     >  1:22530/1‑172 (MQ=255)                                                                                                        gctTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAACTTCTTTAGCAAATTCATCTTTATTAAAAGTTTTTAATTGGTTTatat                                                                                                                     <  1:394549/172‑1 (MQ=255)                                                                                                        gctTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAACTTCTTTAGCAAATTCATCTTTATTAAAAGTTTTTAATTGGTTTatat                                                                                                                     <  1:19918/172‑1 (MQ=255)                                                                                                        gctTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAACTTCTTTAGCAAATTCATCTTTATTAAAAGTTTTTAATTGGTTTATATTGATCACATTTAAAGTATCTTCAGCCTCTTTAATACTTCTACTATCCATAGAGCTATAAAAAACGGTTATATTCCCaaa                                      >  1:895954/1‑251 (MQ=255)                                                                                                        gctTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAACTTCTTTAGCAAATTCATCTTTATTAAAAGTTTTTAATTGGTTTATATTGATCACATTTAAAGTATCTTCAGCCTCTTTAATACTTCTACTATCCATAGAGCTATAAAAAACGGTTATATTCCCaaa                                      >  1:450807/1‑251 (MQ=255)                                                                                                        gctTTAGACTAATCTAATTTAAACCCAACTTCTT‑‑‑‑AAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGt                                                                                                                                                                                                <  1:279750/93‑1 (MQ=255)                                                                                                             aGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGaa                                                                                                                                                                                                           <  1:41221/81‑1 (MQ=255)                                                                                                             aGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAata                                                                                                                                                                            <  1:593617/112‑1 (MQ=255)                                                                                                              gACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCa                                                                                                                                                                                                                   <  1:467689/72‑1 (MQ=255)                                                                                                              gACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAACTTCTTTAGCGAATTCATCTTTATTAAAAg                                                                                                                                       <  1:707358/148‑1 (MQ=255)                                                                                                                  aatctaatTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAACTTCTTTAGCAAATTCATCTTTATTAAAAGTTTTTAATTGGTTTATATTGATc                                                                                                                >  1:238341/1‑167 (MQ=255)                                                                                                                  aatctaatTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAACTTCTTTAGCAAATTCATCTTTATTAAAAGTTTTTAATTGGTTTATATTGATCACATTTAAAGTATCTTCAGCCTCTTTAATACTTCTACTATCCATAGAGCTATAAAAAACGGTTATATTCCCAAAGCCCTGATa                             >  1:726123/1‑250 (MQ=255)                                                                                                                   atctaatTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAAc                                                                                                                                                                    >  1:184787/1‑114 (MQ=255)                                                                                                                                 cAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAACTTCTTTAGCAAATTCATCTTTATTAAAAGTTTTTAATTGGTTTATATTGATCACATTTAAAGTATCTTCAGCCTCTTTAATACTTCTACTATCCATAGAGCTATAAAAAACGGTTATATTCCCAAAGCCCTGATAATGGCGTACCTTATGc             <  1:522455/251‑1 (MQ=255)                                                                                                                                       cttACGAAAGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAACTTCTTTAGCAAATTCATCTTTATTAAAAGTTTTTAATTGGTTTATATTGATCACATTTAAAGTATCTTCAGCCTCTTTAATACTTCTACTATCCATAGAGCTATAAAAAACGGTTATATTCCCAAAGCCCTGATAATGGCGTACCTTATGCCACaa        >  1:106563/1‑250 (MQ=255)                                                                                                                                             aaaGAATTCGAGTAAATTCTCGTAATTTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAACTTCTTTAGCAAATTCATCTTTATTAAAAGTTTTTAATTGGTTTATATTGATCACATTTAAAGTATCTTCAGCCTCTTTAATACTTCTACTATCCATAGAGCTATAAAAAACGGTTATATTCCCAAAGCCCTGATAATGGCGTACCTTATGCCACAACTCAtt  <  1:501545/250‑1 (MQ=255)                                                                                                                                                                       |                                                                                                                                                                                                                                TATTTCATCTTGCAATACAATTCTTTGAATAGGATATGAATTTTTCAATAAACTAATTGAATTTACCTCAAAATCTCTTGCTAAATATATTTCTGTTAGAGCTGCTTTAGACTAATCTAATTTAAACCCAACTTCTTACGAAAGAATTCGAGTAAATTCTCGTAATGTTCGACTATTCCCATCAGGGAATGGGTGAATGTAATCTAGCTTTTTATAAATATCTACAACTTCTTTAGCAAATTCATCTTTATTAAAAGTTTTTAATTGGTTTATATTGATCACATTTAAAGTATCTTCAGCCTCTTTAATATTTCTACTATCCATAGAGCTATAAAAAACGGTTATATTCCCAAAGCCCTGATAATGGCGTACCTTATGCCACAACTCATT  >  1/493204‑493593

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