Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,392,323	G→T	intergenic (+174/‑37)	orn → / → glyV	oligoribonuclease/tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,392,323	0	G	T	100.0%	93.2 / NA	27	intergenic (+174/‑37)	orn/glyV	oligoribonuclease/tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (13/14);  total (13/14)

CTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC  >  NC_000913/4392206‑4392458                                                                                                                      |                                                                                                                                        cTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAACTTTTTTGGGGGTTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCa                                                                                                                   <  2:759657/140‑1 (MQ=255) cTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCa                                                                                                                   <  2:416330/140‑1 (MQ=255) cTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCa                                                                                                                   <  1:45390/140‑1 (MQ=255)         aaCGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAag                                                                                                           <  1:430944/140‑1 (MQ=255)            gTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAagtag                                                                                                        <  2:210562/140‑1 (MQ=255)                tGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTAcg                                                                                                    >  1:355540/1‑140 (MQ=255)                 gCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTAcgc                                                                                                   <  1:90236/140‑1 (MQ=255)                               cACTTGAGATAAAAACGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGt                                                                                     >  1:266508/1‑140 (MQ=255)                                    gagaTAAAAACGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTa                                                                                >  1:764228/1‑140 (MQ=255)                                    gagaTAAAAACGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTa                                                                                >  2:377921/1‑140 (MQ=255)                                         aaaaaCGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCa                                                                           <  2:200440/140‑1 (MQ=255)                                            aaCGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGa                                                                        >  2:573718/1‑140 (MQ=255)                                             aCGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGAc                                                                       <  2:3452/140‑1 (MQ=255)                                             aCGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGAc                                                                       <  1:573718/140‑1 (MQ=255)                                                cAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCtt                                                                    >  1:125970/1‑140 (MQ=255)                                                   aaaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGcc                                                                 <  1:377921/140‑1 (MQ=255)                                                                                aaaGATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGt                                    >  1:688531/1‑140 (MQ=255)                                                                                 aaGATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGtt                                   >  1:526317/1‑140 (MQ=255)                                                                                   gATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCC‑AGGTCGGGGTCGCGAGTTc                                             >  1:528348/1‑127 (MQ=255)                                                                                   gATTTCTCGTATAATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCC‑AGGTCGGGGTCGCGAGTTc                                             <  2:528348/127‑1 (MQ=255)                                                                                              taATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTcc                                >  2:260746/1‑130 (MQ=255)                                                                                              taATGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTcc                                <  1:260746/130‑1 (MQ=255)                                                                                                aTGCGCCTCCCGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAAtt                    >  1:121093/1‑140 (MQ=255)                                                                                                        cccGTAACGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGt            >  2:22149/1‑140 (MQ=255)                                                                                                             aaCGACGCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGtgctgc       <  1:313115/140‑1 (MQ=255)                                                                                                              acgacgCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCa      <  2:8032/140‑1 (MQ=255)                                                                                                                  cgCATAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGc  >  1:352247/1‑140 (MQ=255)                                                                                                                      |                                                                                                                                        CTAAATTGAACGTCTTGCTGATAATTTCAGCACTTGAGATAAAAACGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC  >  NC_000913/4392206‑4392458

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