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

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

AAAAATCTGTTTTTATAACACGAATAAATTTTACAAAAGATCCAACTGTGGATAAATAACAGAGTTATACACAAGAAATTTGAGGAGGATCGGTAAGATCTTAACAACTTGTAGCATGATCCATCTCCTTCATTTACAACATAAAAAAGATCTTATCCAAGGTTTTAGTCGATCTTAATTATAGTAATAATAAACTTTATATATAGATAAAGATCTTATTACTTATTAACGTCGGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATGTTTTATACAGAAACCTATGATGTCATCGTAATCGGTGGCGGTCATGCAGGAACCGAAGCCGCCTTGGCACCTGCGCGCATGGGGTTAAAAACC  >  1/477532‑478001                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             aaaaaTCTGTTTTTATAACACGAATAAATTTTACAAAAGATCCAACTGTGGATAAATAACAGAGTTATACACAAGAAATTTGAGGAGGATCGGTAAGATCTTAACAACTTGTAGCATGATCCATCTCCTTCATTTACAACATAAAAAAGATCTTATCCAAGGTTTTCGTCGATCTTAATTATAGTAATAATAAACTTTATATATAGATAAAGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATc                                                                                                                                                                                                                             <  1:193102/251‑1 (MQ=255)              tataACACGAATAAATTTTACAAAAGATCCAACTGTGGATAAATAACAGAGTTATACACAAGAAATTTGAGGAGGATCGGTAAGATCTTAACAACTTGTAGCATGATCCATCTCCTTCATTTACAACATAAAAAAGATCTTATCCAAGGTTTTCGTCGATCTTAATTATAGTAATAATAAACTTTATATATAGATAAAGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTc                                                                                                                                                                                                                <  1:778382/251‑1 (MQ=255)                                      aGATCCAACTGTGGATAAATAACAGAGTTATACACAAGAAATTTGAGGAGGATCGGTAAGATCTTAACAACTTGTAGCATGATCCATCTCCTTCATTTACAACATAAAAAAGATCTTATCCAAGGTTTTCGTCGATCTTAATTATAGTAATAATAAACTTTATATATAGATAAAGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATc                                                                                                                                                                                                                             <  1:665302/214‑1 (MQ=255)                                                                                                                       aTCCATCTCCTTCATTTACAACATAAAAAAGATCTTATCCAAGGTTTTCGTCGATCTTAATTATAGTAATAATAAACTTTATATATAGATAAAGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATaaaaaa                                                                                                       >  1:324945/1‑251 (MQ=255)                                                                                                                                                                            aTCTTAATTATAGTAATAATAAACTTTATATATAGATAAAGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAAc                                                                                                                                                                                                       <  1:720056/102‑1 (MQ=255)                                                                                                                                                                                                  aCTTTATATATAGATAAAGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGaa                                                                                                                           >  1:585316/1‑156 (MQ=255)                                                                                                                                                                                                  aCTTTATATATAGATAAAGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATGTTTTATAc                                                                                       <  1:134639/192‑1 (MQ=255)                                                                                                                                                                                                  aCTTTATATATAGATAAAGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATGTTTTATACAGAAACCTatga                                                                           >  1:297865/1‑204 (MQ=255)                                                                                                                                                                                                   cTTTATATATAGATAAAGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATGTTTTATACAGAAACCTATGATGTCATCGTAATCGGTGGCGGTCATGCAGGAACCGAAGCCGCCTTgg                            >  1:243386/1‑250 (MQ=255)                                                                                                                                                                                                   cTTTATATATAGATAAAGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATGTTTTATACAGAAACCTATGATGTCATCGTAATCGGTGGCGGTCATGCAGGAACCGAAGCCGCCTTgg                            >  1:564156/1‑250 (MQ=255)                                                                                                                                                                                                       atatataGATAAAGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTAc                                                                                                                               <  1:363298/147‑1 (MQ=255)                                                                                                                                                                                                       atatataGATAAAGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATGTTTTATAc                                                                                       >  1:536016/1‑187 (MQ=255)                                                                                                                                                                                                       atatataGATAAAGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATGTTTTATACAGAAACCTATGATGTCATCGTAATCGGTGGCGGTCATGCAGGAACCGAAGCCGCCTTGGCAc                         >  1:269457/1‑249 (MQ=255)                                                                                                                                                                                                             aGATAAAGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATGTTTTATACAGAAACCTATGATGTCATCGTAAt                                                               >  1:169820/1‑205 (MQ=255)                                                                                                                                                                                                              gATAAAGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAAt                                                                                                                          >  1:701013/1‑145 (MQ=255)                                                                                                                                                                                                              gATAAAGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAag                                                                                                                 <  1:560274/154‑1 (MQ=255)                                                                                                                                                                                                              gATAAAGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATg                                                                                               <  1:291171/172‑1 (MQ=255)                                                                                                                                                                                                              gATAAAGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATGTTTTATACAGAAACCTatgatg                                                                         >  1:633592/1‑194 (MQ=255)                                                                                                                                                                                                                  aaGATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAAc                                                                                                                                                                                                       >  1:595223/1‑64 (MQ=255)                                                                                                                                                                                                                    gATCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATGTTTTATACAGAAACCTATGATGTCATCGTAATCGGTGGCGGTCATGCAGGAACCGAAgc                                    >  1:512330/1‑225 (MQ=255)                                                                                                                                                                                                                     aTCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATGTTTTATACAGAAACCTATGATGTCATCGTAATCGGTGGCGGTCATGCAGGTACCGAAGCCGCCTTGGCACCTGCGCGCATGGGGtt         >  1:861229/1‑251 (MQ=255)                                                                                                                                                                                                                      tCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAAc                                                                                                                                                                                                       <  1:685444/60‑1 (MQ=255)                                                                                                                                                                                                                      tCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCt                                                                                                                                                                                                 >  1:509614/1‑66 (MQ=255)                                                                                                                                                                                                                      tCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCg                                                                                                                                   >  1:13853/1‑128 (MQ=255)                                                                                                                                                                                                                      tCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAg                                                                                                                             <  1:60085/134‑1 (MQ=255)                                                                                                                                                                                                                      tCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATg                                                                                               >  1:894060/1‑164 (MQ=255)                                                                                                                                                                                                                      tCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATGTTTTATAc                                                                                       >  1:607276/1‑172 (MQ=255)                                                                                                                                                                                                                      tCTTATTACTTATTAACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATGTTTTATACAGAAACCTATGATGTCATCGTAATCGGTGGCGGTCATGCAGGAAc                                          <  1:291089/217‑1 (MQ=255)                                                                                                                                                                                                                             acttattaACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATGTTTTATACAGAAACCTATGATGTCATCGTAATCGGTGGCGGTCATGCAGGAACCGAAGCCGCCTTGGCACCTGCGCGCATGGGGTTaaaa     >  1:636476/1‑247 (MQ=255)                                                                                                                                                                                                                             acttattaACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATGTTTTATACAGAAACCTATGATGTCATCGTAATCGGTGGCGGTCATGCAGGAACCGAAGCCGCCTTGGCACCTGCGCGCATGGGGTTAAAAAcc  >  1:753122/1‑250 (MQ=255)                                                                                                                                                                                                                              cttattaACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATGTTTTATACAGAAACCTatgatg                                                                         >  1:720855/1‑178 (MQ=255)                                                                                                                                                                                                                                 attaACGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCcac                                                                                                                                                                                        >  1:64906/1‑64 (MQ=255)                                                                                                                                                                                                                                      cGTCCGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATGTTTTATAc                                                                                       <  1:472459/156‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             AAAAATCTGTTTTTATAACACGAATAAATTTTACAAAAGATCCAACTGTGGATAAATAACAGAGTTATACACAAGAAATTTGAGGAGGATCGGTAAGATCTTAACAACTTGTAGCATGATCCATCTCCTTCATTTACAACATAAAAAAGATCTTATCCAAGGTTTTAGTCGATCTTAATTATAGTAATAATAAACTTTATATATAGATAAAGATCTTATTACTTATTAACGTCGGATCCTTTCTTCAGATCCCGTCTTCTATTCCTCTCAAACTGACCTTTTATCCACACGGGATTTTCAGTAGAATAGCAGCCAATTTTTTTAAAGGATCAGATGATCCGTTACAGAATGCAAAAAAGAGATAAAAAATAATCATGTTTTATACAGAAACCTATGATGTCATCGTAATCGGTGGCGGTCATGCAGGAACCGAAGCCGCCTTGGCACCTGCGCGCATGGGGTTAAAAACC  >  1/477532‑478001

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