Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	973,577	C→T	S14F (TCC→TTC)	cmoM →	tRNA cmo(5)U34 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	973,577	0	C	T	100.0%	88.1 / NA	31	S14F (TCC→TTC)	cmoM	tRNA cmo(5)U34 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/14);  total (17/14)

CAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAGACGGGTTCATTCTCGCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGG‑AGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATCTTTCTGCGCAGATGATCGACCGCGCAAAACAGGCGGCAGAAGCAAAAGGTGTG  >  NC_000913/973332‑973809                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          cAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAGACGGGTTCATTCTCGCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTa                                                                                                                                                                                                                                      <  1:126713/251‑1 (MQ=255)                tCAGCGGAAGTGGCAGCATAATCTGGCCAATCACTTTTTTCAGTGTAAAAACCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAGACGGGCTCATTCTCGCCGAGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACAGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGGCGAATAGTTTTTCCGTAACATTTACGGcacc                                                                                                                                                                                                                        >  2:236375/1‑250 (MQ=255)                     ggAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAGACGGGTTCATTCTCGCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCaa                                                                                                                                                                                                                   >  2:40032/1‑250 (MQ=255)                     ggAAGGGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAACACGCATCCTTTCTGGTTCCTTTTTTAACAATATAGCAAGGGATCGCCAGCGATATTACACCAGACGGGTTCATTCTCGCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCCCCaa                                                                                                                                                                                                                   >  2:132630/1‑250 (MQ=255)                            gcagcaACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAGACGGGTTCATTCTCGCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGc                                                                                                                                                                                                           <  1:260715/251‑1 (MQ=255)                              agcaACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAGACGGGTTCATTCTCGCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCtt                                                                                                                                                                                                         >  1:21855/1‑251 (MQ=255)                               gcaACATATTTCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAGACGGGTTCATTCTCGCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTc                                                                                                                                                                                                        <  1:236375/251‑1 (MQ=255)                                        tGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCCTATAGCAAGGGATCGCCAGCGATATTACACCAGACGGGTTCATTCTCGCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCCACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGACAGCTTCGACAGGCtt                                                                                                                                                                                               >  1:49891/1‑250 (MQ=255)                                               cccTTTTTTCAGTGTAAAAAGCATCATTTCTGGGTCCTTTTTTAACCATATAGCCAGGAATCGCCAGCGATATTACACCAGACGGGTTCATTCTCGCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGGGGAGCAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCAACAAAGGGCAGATTCGACAGGCTaatatggg                                                                                                                                                                                        >  2:229343/3‑244 (MQ=255)                                                          gtgtAACAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAGACGGGTTCATTCTCGCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCg                                                                                                                                                                             >  1:348918/1‑251 (MQ=255)                                                                                     tttttAACCATATAGCAAGGGAGCGCCAGCGATATTACACCAGACGGGTTCATTCTCGCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGGGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCGCGATCGCGTGCTGGCGGAAATGGGCCCg                                                                                                                                                   <  1:179921/250‑1 (MQ=255)                                                                                          aaCCATATAGCAAGGGATCGCCAGCGATATTACACCAGACGGGTTCATTCTCGCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCaaaaa                                                                                                                                             <  1:343792/251‑1 (MQ=255)                                                                                                              ccAGCGATATTACACCAGACGGGTTCATTCTCGCCGTGGCTGTGACCAAATAGCGTCTTTATCAGAACACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCCAGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGGGTGCTGGAtggtgg                                                                                                                         >  2:300043/1‑251 (MQ=255)                                                                                                                         acacCAGACGGGTTCATTCTCGCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGG‑AGAAg                                                                                                             <  1:300043/251‑1 (MQ=255)                                                                                                                              aGACGGGTTCATTCTCGCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGG‑AGAAGGGCAg                                                                                                        <  1:40032/251‑1 (MQ=255)                                                                                                                                        aTTCTCGTCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGTAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCTCGCAAAAACTGCGTGTGCTGGATGCTGGAGGTGG‑AGAATGGCAGACCGCAATc                                                                                               <  1:277862/250‑1 (MQ=255)                                                                                                                                                        tgACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCAGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATAGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGG‑AGCAGGGCAGACCGCAATCAAAATGGCCGAGCgggg                                                                              >  1:201064/1‑248 (MQ=255)                                                                                                                                                               aTAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGGGCTGGATGCTGGCGGTGGTAGAAGG‑CACACCGCAATAAAAATggtcgggccgtgggatccg                                                                       >  2:124803/1‑234 (MQ=255)                                                                                                                                                                tAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGG‑AGAAGGGCAGACCGCAATCAAAAAGGCCGAGCGGGGGCATCAgg                                                                      >  1:222580/1‑251 (MQ=255)                                                                                                                                                                        ttATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGG‑AGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCAttt                                                                <  1:132630/249‑1 (MQ=255)                                                                                                                                                                              gAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGG‑AGGAAGGCAGACCGCAATAAAAATGGCCGAGCGTGGGCATCAGGTCCTTTTATGCGAt                                                        >  1:248927/1‑251 (MQ=255)                                                                                                                                                                                       ttAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGG‑GGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTAATTTTAAGCGATCTTTCTgcg                                               >  1:279072/1‑251 (MQ=255)                                                                                                                                                                                              aCCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGATGGATGCTGGCGGTGG‑GGAAGGGCAGACCGCAATCAAAATGGCCGCGCGTGGGCATCAGGTCATTTTATGCGATCTTTCTGCGCAgatg                                         >  1:334081/1‑250 (MQ=255)                                                                                                                                                                                                                tcAGGATCGCCATTTTGATGATATTTCTGTAAAGTTTTTCCGTTACATTTACGGCACCACCAAAGGGCAGCTTCGACCGGCTATTCTGTGGCCGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCTTGTGCTGGATGCTGGCGGTGG‑AGAAGGGCTGACCGCAATCAAAATGGCCGAGCCTGGGCATCAGGTCATTTTATGCGATCTTTCTGCGCATATGATCGACCGCGCAAAAc                         <  2:216181/247‑1 (MQ=255)                                                                                                                                                                                                                gCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGGGG‑AGAAGGGGAGAACGCAATCAAAATGGCCGAGCGGGGGCATCAGGGCATTTTATGGGATCTTTCTGCGCAgatgggcgaacgggcaaaaa                         >  1:216181/1‑232 (MQ=255)                                                                                                                                                                                                                  aGGATCGCAATTTAGATTATAGTGCGGAAAAGTTTTTCCGTAACATCTACGGCACCACCAAAGGGCAGCGACGACAGGCTAGTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGG‑AGAAGGGCAGACCGCAATCAAAAGGGCCGAGCGTGGGTATCAGGTCATTTGATGCGATCTTTCTGCGCAGATGATCGACCGCGCAAAACAggag                    <  2:167815/251‑3 (MQ=255)                                                                                                                                                                                                                          aatgttGATGATATGGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGAGCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGGGTGCTGGATGCTGGCGGTGG‑AGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATCTTTCTGCGCAGATGATCGACCGCGCAAAACAGGCGGCATAAGc            <  2:260183/247‑1 (MQ=255)                                                                                                                                                                                                                            ttttGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGG‑AGAAGGGCAGACCGCAATCAAAATGGCCGGGCGTGGGCATCAGGGCATTTTATGCGATCTTTCTGCGCAGATGATCGACCGCGCAAAACAGGCGGCAGAAGcc           >  1:208379/1‑249 (MQ=255)                                                                                                                                                                                                                             tttGATGATATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGG‑AGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATCTTTCTGCGCAGATGATCGACCGCGCCAAACAGGCGGCAGAAGCaa          >  1:159699/1‑250 (MQ=255)                                                                                                                                                                                                                                   gatATTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGG‑AGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATCTTTCTGCGCAGATGATCGACCGCGCAAAACAGGCGGCAGAAGCAAAAGgtgt   <  1:37532/251‑1 (MQ=255)                                                                                                                                                                                                                                    atatTGCGGAAAAGTTTTTCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGG‑AGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCAGTTTATGCGATCTTTCTGCGCAGATGATCGACCGCGCAAAACAGGCGGCAGAAGCAAAAGgtggg  <  2:132194/251‑3 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          CAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAGACGGGTTCATTCTCGCCGTGGCTGTGACAAAATAGCGTCTTTATCAGAAAACGTCTTAGCCAACCGGTGGAGAAAGCAATGCAGGATCGCAATTTTGATGATATTGCGGAAAAGTTTTCCCGTAACATTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGG‑AGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATCTTTCTGCGCAGATGATCGACCGCGCAAAACAGGCGGCAGAAGCAAAAGGTGTG  >  NC_000913/973332‑973809

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