Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,162,672	A→G	L200P (CTG→CCG)	trmA ←	tRNA m(5)U54 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,162,672	0	A	G	100.0%	140.5 / NA	40	L200P (CTG→CCG)	trmA	tRNA m(5)U54 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/18);  total (22/18)

ACCGACGGCTTAGCGATTTCGGTGGCTAATACCCGATCAAAATTACGCGCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCC  >  NC_000913/4162537‑4162806                                                                                                                                        |                                                                                                                                       aCCGACGGCTTAGCGATTTCGGTGGCTAATACCCGATCAAAATTACGCGCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCAt                                                                                                                                    >  1:163488/1‑140 (MQ=255)         cTTAGCGATTTCGGTGGCTAATACCCGATCAAAATTACGCGCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAatat                                                                                                                            >  1:599467/1‑140 (MQ=255)             gCGATTTCGGTGGCTAATACCCGATCAAAATTACGCGCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCAt                                                                                                                        >  2:109037/1‑140 (MQ=255)                      gTGGCTAATACCCGATCAAAATTACGCGCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGtt                                                                                                               <  2:599467/140‑1 (MQ=255)                       tGGCTAATACCCGATCAAAATTACGCGCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGt                                                                                                                >  1:415831/1‑138 (MQ=255)                       tGGCTAATACCCGATCAAAATTACGCGCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGt                                                                                                                <  2:415831/138‑1 (MQ=255)                         gCTAATACCCGATCAAAATTACGCGCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCgg                                                                                                            <  1:73266/140‑1 (MQ=255)                         gCTAATACCCGATCAAAATTACGCGCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCgg                                                                                                            <  2:482271/140‑1 (MQ=255)                                   gATCAAAATTACGCGCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGtt                                                                                                               <  2:92608/127‑1 (MQ=255)                                   gATCAAAATTACGCGCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGtt                                                                                                               >  1:92608/1‑127 (MQ=255)                                   gATCAAAATTACGCGCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAg                                                                                                  <  1:154895/140‑1 (MQ=255)                                   gATCAAAATTACGCGCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAg                                                                                                  >  1:454654/1‑140 (MQ=255)                                    aTCAAAATTACGCGCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGc                                                                                                 >  2:302024/1‑140 (MQ=255)                                              cgcgCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTAc                                                                                       <  2:275442/140‑1 (MQ=255)                                               gcgcCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTAcc                                                                                      <  2:34534/140‑1 (MQ=255)                                                cgcCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTT‑AGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTg                                                                                    >  1:626706/1‑140 (MQ=255)                                                  ccAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGa                                                                                   <  1:47811/140‑1 (MQ=255)                                                           aTGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGatc                                                                          <  2:454654/140‑1 (MQ=255)                                                                     aCCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTcc                                                                >  2:501633/1‑140 (MQ=255)                                                                                      gCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCgg                                                      >  1:802493/1‑133 (MQ=255)                                                                                      gCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCgg                                                      <  2:802493/133‑1 (MQ=255)                                                                                       cTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGtt                                              >  1:102897/1‑140 (MQ=255)                                                                                       cTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGtt                                              >  1:150910/1‑140 (MQ=255)                                                                                          cAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCAt                                           <  1:109037/140‑1 (MQ=255)                                                                                          cAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCAt                                           >  2:167599/1‑140 (MQ=255)                                                                                          cAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCAt                                           <  1:167599/140‑1 (MQ=255)                                                                                            gTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCg                                         <  2:163488/140‑1 (MQ=255)                                                                                                aTCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGt                                     <  2:102897/140‑1 (MQ=255)                                                                                                aTCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGt                                     <  1:777371/140‑1 (MQ=255)                                                                                                 tCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTa                                    >  2:482055/1‑140 (MQ=255)                                                                                                    ccTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATc                                 >  1:268810/1‑140 (MQ=255)                                                                                                    ccTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATc                                 >  1:145422/1‑140 (MQ=255)                                                                                                               tttGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCt                      <  2:150910/140‑1 (MQ=255)                                                                                                                 tGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGt                                     >  1:713306/1‑123 (MQ=255)                                                                                                                 tGGTTACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGt                                     <  2:713306/123‑1 (MQ=255)                                                                                                                     tACGTCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGAttt                >  1:576914/1‑140 (MQ=255)                                                                                                                         tCCAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGt            >  2:357695/1‑140 (MQ=255)                                                                                                                           cAGCGCCCATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTATCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGttt          >  1:469579/1‑140 (MQ=255)                                                                                                                                  cATTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGGCGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGc   >  2:556277/1‑140 (MQ=255)                                                                                                                                   aTTCCGGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGcc  >  2:122003/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       ACCGACGGCTTAGCGATTTCGGTGGCTAATACCCGATCAAAATTACGCGCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCC  >  NC_000913/4162537‑4162806

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