Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	2	124,839	C→G	intergenic (‑141/+302)	unknown ← / ← unknown	tRNA‑5‑carboxymethylaminomethyl‑2‑ thiouridine(34) synthesis protein MnmG/Flavoprotein MioC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	2	124,839	0	C	G	100.0%	92.0 / NA	35	intergenic (‑141/+302)	unknown/unknown	tRNA‑5‑carboxymethylaminomethyl‑2‑ thiouridine(34) synthesis protein MnmG/Flavoprotein MioC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (14/21);  total (14/21)

GGTTTTTAACCCCATGCGCGCAGGTGCCAAGGCGGCTTCGGTTCCTGCATGACCGCCACCGATTACGATGACATCATAGGTTTCTGTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCCGACGTTAATAAGTAATAAGATCTTTATCTATATATAAAGTTTATTATTACTATAATTAAGATCGACGAAAACCTTGGATAAGATCTTTTTTATGTTGTAAATGAAGGAGATGGATCATGCTACAAGTTGTTAAGATCTTACCGATCCTCCTCAAATTTCTTGTGTATAACTCTGTTATTTATCCACAGTTGGATCTTTTGTAAAATTTATTCGTGTTATAAAAACAGATTTTT  >  2/124603‑125072                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          ggTTTTTAACCCCATGCGCGCAGGTGCCAAGGCGGCTTCGGTTCCTGCATGACCGCCACCGATTACGATGACATCATAGGTTTCTGTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTtaataagt                                                                                                                                                                                                                              <  1:753122/250‑1 (MQ=255)    ttttAACCCCATGCGCGCAGGTGCCAAGGCGGCTTCGGTTCCTGCATGACCGCCACCGATTACGATGACATCATAGGTTTCTGTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTtaataagt                                                                                                                                                                                                                              <  1:636476/247‑1 (MQ=255)        aaCCCCATGCGCGCAGGTGCCAAGGCGGCTTCGGTACCTGCATGACCGCCACCGATTACGATGACATCATAGGTTTCTGTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGAt                                                                                                                                                                                                                      <  1:861229/251‑1 (MQ=255)                        gTGCCAAGGCGGCTTCGGTTCCTGCATGACCGCCACCGATTACGATGACATCATAGGTTTCTGTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCTTTATCtatatat                                                                                                                                                                                                        <  1:269457/249‑1 (MQ=255)                           ccAAGGCGGCTTCGGTTCCTGCATGACCGCCACCGATTACGATGACATCATAGGTTTCTGTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCTTTATCTATATATAAAg                                                                                                                                                                                                    <  1:564156/250‑1 (MQ=255)                           ccAAGGCGGCTTCGGTTCCTGCATGACCGCCACCGATTACGATGACATCATAGGTTTCTGTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCTTTATCTATATATAAAg                                                                                                                                                                                                    <  1:243386/250‑1 (MQ=255)                                   gcTTCGGTTCCTGCATGACCGCCACCGATTACGATGACATCATAGGTTTCTGTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATc                                                                                                                                                                                                                     <  1:512330/225‑1 (MQ=255)                                         gTTCCTGCATGACCGCCACCGATTACGATGACATCATAGGTTTCTGTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGa                                                                                                                                                                                                                       >  1:291089/1‑217 (MQ=255)                                                              aTTACGATGACATCATAGGTTTCTGTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCTTTATCt                                                                                                                                                                                                              <  1:169820/205‑1 (MQ=255)                                                                       aCATCATAGGTTTCTGTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGAc                                                                                                                                                                                                                                        >  1:377336/1‑170 (MQ=255)                                                                        catcatAGGTTTCTGTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTtaataag                                                                                                                                                                                                                               <  1:720855/178‑1 (MQ=255)                                                                        catcatAGGTTTCTGTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCTTTATc                                                                                                                                                                                                               <  1:633592/194‑1 (MQ=255)                                                                          tcatAGGTTTCTGTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCTTTATCTATATATAAAGt                                                                                                                                                                                                   <  1:297865/204‑1 (MQ=255)                                                                                      gTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACg                                                                                                                                                                                                                                       >  1:472459/1‑156 (MQ=255)                                                                                      gTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGa                                                                                                                                                                                                                       <  1:607276/172‑1 (MQ=255)                                                                                      gTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCTTTATCtatatat                                                                                                                                                                                                        <  1:536016/187‑1 (MQ=255)                                                                                      gTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCTTTATCTATATATAAAGt                                                                                                                                                                                                   >  1:134639/1‑192 (MQ=255)                                                                                              cATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGa                                                                                                                                                                                                                       <  1:894060/164‑1 (MQ=255)                                                                                              cATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCTTTATc                                                                                                                                                                                                               >  1:291171/1‑172 (MQ=255)                                                                                                      ttttttATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCTTTATCTATATATAAAGTTTATTATTACTATAATTAAGATCGACGAAAACCTTGGATAAGATCTTTTTTATGTTGTAAATGAAGGAGATGGAt                                                                                                                        <  1:324945/251‑1 (MQ=255)                                                                                                                ctTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCTTTATc                                                                                                                                                                                                               >  1:560274/1‑154 (MQ=255)                                                                                                                         aTTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCTTTATc                                                                                                                                                                                                               <  1:701013/145‑1 (MQ=255)                                                                                                                          ttCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCTTTATCTATATATAAAGt                                                                                                                                                                                                   <  1:585316/156‑1 (MQ=255)                                                                                                                            cTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGAc                                                                                                                                                                                                                                        >  1:645429/1‑117 (MQ=255)                                                                                                                            cTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGa                                                                                                                                                                                                                       >  1:60085/1‑134 (MQ=255)                                                                                                                              gTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCTTTATCtatatat                                                                                                                                                                                                        >  1:363298/1‑147 (MQ=255)                                                                                                                                  cGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGa                                                                                                                                                                                                                       <  1:13853/128‑1 (MQ=255)                                                                                                                                                                                       gtgGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTtaat                                                                                                                                                                                                                                  <  1:64906/64‑1 (MQ=255)                                                                                                                                                                                                aGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGa                                                                                                                                                                                                                       <  1:509614/66‑1 (MQ=255)                                                                                                                                                                                                      gTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGa                                                                                                                                                                                                                       >  1:685444/1‑60 (MQ=255)                                                                                                                                                                                                      gTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCtt                                                                                                                                                                                                                   <  1:595223/64‑1 (MQ=255)                                                                                                                                                                                                      gTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCTTTATCTATATATAAAGTTTATTATTACTATAATTAAGAt                                                                                                                                                                             >  1:720056/1‑102 (MQ=255)                                                                                                                                                                                                               gAATAGAAGACGGGATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCTTTATCTATATATAAAGTTTATTATTACTATAATTAAGATCGACGAAAACCTTGGATAAGATCTTTTTTATGTTGTAAATGAAGGAGATGGATCATGCTACAAGTTGTTAAGATCTTACCGATCCTCCTCAAATTTCTTGTGTATAACTCTGTTATTTATCCACAGTTGGATCTTTTGTAAAATTTATTCGTGTtata               >  1:778382/1‑251 (MQ=255)                                                                                                                                                                                                                            gATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCTTTATCTATATATAAAGTTTATTATTACTATAATTAAGATCGACGAAAACCTTGGATAAGATCTTTTTTATGTTGTAAATGAAGGAGATGGATCATGCTACAAGTTGTTAAGATCTTACCGATCCTCCTCAAATTTCTTGTGTATAACTCTGTTATTTATCCACAGTTGGATCt                                       >  1:665302/1‑214 (MQ=255)                                                                                                                                                                                                                            gATCTGAAGAAAGGATCGGACGTTAATAAGTAATAAGATCTTTATCTATATATAAAGTTTATTATTACTATAATTAAGATCGACGAAAACCTTGGATAAGATCTTTTTTATGTTGTAAATGAAGGAGATGGATCATGCTACAAGTTGTTAAGATCTTACCGATCCTCCTCAAATTTCTTGTGTATAACTCTGTTATTTATCCACAGTTGGATCTTTTGTAAAATTTATTCGTGTTATAAAAACAGAttttt  >  1:193102/1‑251 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          GGTTTTTAACCCCATGCGCGCAGGTGCCAAGGCGGCTTCGGTTCCTGCATGACCGCCACCGATTACGATGACATCATAGGTTTCTGTATAAAACATGATTATTTTTTATCTCTTTTTTGCATTCTGTAACGGATCATCTGATCCTTTAAAAAAATTGGCTGCTATTCTACTGAAAATCCCGTGTGGATAAAAGGTCAGTTTGAGAGGAATAGAAGACGGGATCTGAAGAAAGGATCCGACGTTAATAAGTAATAAGATCTTTATCTATATATAAAGTTTATTATTACTATAATTAAGATCGACGAAAACCTTGGATAAGATCTTTTTTATGTTGTAAATGAAGGAGATGGATCATGCTACAAGTTGTTAAGATCTTACCGATCCTCCTCAAATTTCTTGTGTATAACTCTGTTATTTATCCACAGTTGGATCTTTTGTAAAATTTATTCGTGTTATAAAAACAGATTTTT  >  2/124603‑125072

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