Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	218,835	(A)8→7	intergenic (‑60/+52)	proS ← / ← trmO	proline‑‑tRNA ligase/tRNA m(6)t(6)A37 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	218,828	0	A	.	100.0%	154.0 / NA	36	intergenic (‑53/+59)	proS/trmO	proline‑‑tRNA ligase/tRNA m(6)t(6)A37 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (21/15);  total (21/15)

CGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGA  >  NC_000913/218693‑218963                                                                                                                                        |                                                                                                                                        cgcAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTTAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACaaaaa                                                                                                                                     >  2:168393/1‑140 (MQ=255)  gcagcaTCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACaaaaaa                                                                                                                                    <  2:553329/140‑1 (MQ=255)  gcagcaTCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACaaaaaa                                                                                                                                    <  2:726726/140‑1 (MQ=255)       aTCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACaaaaaaagtga                                                                                                                               <  1:116419/140‑5 (MQ=255)                tGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTa                                                                                                                     >  1:129235/1‑140 (MQ=255)                tGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTa                                                                                                                     >  1:725712/1‑140 (MQ=255)                      aTCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTg                                                                                                               >  1:18416/1‑140 (MQ=255)                        cACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCa                                                                                                             <  1:968822/140‑1 (MQ=255)                                  tCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGt                                                                                                   >  1:1005954/1‑140 (MQ=255)                                          tgtCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCa                                                                                                 <  2:750053/134‑1 (MQ=255)                                          tgtCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCa                                                                                                 >  1:750053/1‑134 (MQ=255)                                           gtCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAgaga                                                                                          <  2:3148/140‑1 (MQ=255)                                                        tGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAAttt                                                                             >  2:511400/1‑140 (MQ=255)                                                          gagaGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTaa                                                                           >  2:277450/1‑140 (MQ=255)                                                          gagaGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTaa                                                                           >  1:786757/1‑140 (MQ=255)                                                                   tATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTcc                                                                  <  2:129235/140‑1 (MQ=255)                                                                         cTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCa                                                            <  1:511400/140‑1 (MQ=255)                                                                         cTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCa                                                            <  2:725712/140‑1 (MQ=255)                                                                          tAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCaa                                                           >  2:1027352/1‑140 (MQ=255)                                                                                  aTGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCa                                                   >  2:938832/1‑140 (MQ=255)                                                                                      tgttACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAAc                                               <  2:842593/140‑1 (MQ=255)                                                                                       gttACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAAcc                                              >  2:990530/1‑140 (MQ=255)                                                                                           cGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCg                                          <  2:18416/140‑1 (MQ=255)                                                                                                ccAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAGAACCTGCGTCGGt                                     >  1:933476/1‑140 (MQ=255)                                                                                                ccAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGt                                     >  2:371724/1‑140 (MQ=255)                                                                                                  cgTTGGAAGGAAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGCGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGa                                   >  2:161804/2‑140 (MQ=255)                                                                                                  aGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGa                                   >  1:322238/1‑140 (MQ=255)                                                                                                            tAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCg                         <  1:938832/140‑1 (MQ=255)                                                                                                            tAGAACAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCg                         <  1:277450/140‑1 (MQ=255)                                                                                                                 cAGGCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAGAACCTGCGTCGGTGACGCGCCAGCGAACGt                    <  2:933476/140‑1 (MQ=255)                                                                                                                    gCTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTaa                 >  1:201911/1‑140 (MQ=255)                                                                                                                     cTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTaaa                >  2:926824/1‑140 (MQ=255)                                                                                                                     cTCAAGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTaaa                <  1:371724/140‑1 (MQ=255)                                                                                                                         aGGCGAGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCa            >  2:1059865/1‑140 (MQ=255)                                                                                                                              aGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCg                         <  2:512379/122‑1 (MQ=255)                                                                                                                              aGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCg                         >  1:512379/1‑122 (MQ=255)                                                                                                                              aGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAg       >  1:500863/1‑140 (MQ=255)                                                                                                                              aGCCTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAg       >  1:429789/1‑140 (MQ=255)                                                                                                                                 cTGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCa    <  2:322238/140‑1 (MQ=255)                                                                                                                                  tGGGAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAg   >  2:232991/1‑140 (MQ=255)                                                                                                                                   gggAC‑AAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGa  <  1:1059865/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        CGCAGCATCAGCTGATGGCTGATCACCTCGGCGTCGGCAGGTGTCTCCTTGAGAGTGGAGAGCAGGTATTGGCTAGTACGCATGTTGTTACGGTTCCAGTTGGAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGA  >  NC_000913/218693‑218963

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