Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,392,453	C→T	intergenic (+18/‑19)	glyV → / → glyX	tRNA‑Gly/tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,392,453	0	C	T	85.7%	55.1 / 2.3	28	intergenic (+18/‑19)	glyV/glyX	tRNA‑Gly/tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (17/7);  total (19/9)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

CGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATA  >  NC_000913/4392304‑4392591                                                                                                                                                      |                                                                                                                                           cgcCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTaa                                                                                                                                          <  2:564363/151‑1 (MQ=255)      cccGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGca                                                                                                                                      <  2:401356/150‑1 (MQ=255)       ccGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGcac                                                                                                                                     >  1:889578/1‑150 (MQ=255)         gTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGcac                                                                                                                                     >  1:408284/1‑148 (MQ=255)              gacgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGaccacc                                                                                                                             >  1:219387/1‑151 (MQ=255)                cgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCa                                                                                                                           <  1:155782/151‑1 (MQ=255)                cgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGATGTAAGGCACAGaccacc                                                                                                                             >  2:143386/1‑149 (MQ=255)                 gCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCAccc                                                                                                                            >  2:901509/1‑149 (MQ=255)                 gCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCa                                                                                                                           >  1:68791/1‑150 (MQ=255)                  cAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGaccacc                                                                                                                             >  2:397831/1‑147 (MQ=255)                   aGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGc                                                                                                                        >  2:976526/1‑151 (MQ=255)                    gAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCg                                                                                                                       >  1:647747/1‑151 (MQ=255)                     aaaTGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCgg                                                                                                                      >  2:114167/1‑151 (MQ=255)                          cGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATa                                                                                                                 <  2:247589/151‑1 (MQ=255)                          cGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGCGCTGTAAGGCACAGAGCACCCAAGCGGGaa                                                                                                                   >  2:79899/1‑149 (MQ=14)                           gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGaa                                                                                                                   >  2:911167/1‑148 (MQ=255)                           gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAg                                                                                                                <  2:40267/151‑1 (MQ=255)                            aaaaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑CAAGCGCTGTAAGGCACAGACCACCCAAGCggg                                                                                                                     >  2:404751/1‑145 (MQ=14)                              aaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTc                                                                                                             <  2:635542/151‑1 (MQ=255)                               aTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTc                                                                                                             >  1:427421/1‑151 (MQ=255)                                    gAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTgg                                                                                                       >  1:447213/1‑151 (MQ=255)                                     aaaGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGCGCTGTAAGGCACAGACCACCCAAGCGCGAAAATCTTAGTTGGt                                                                                                      >  2:847765/1‑151 (MQ=11)                                      aaGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTa                                                                                                     <  1:114175/151‑1 (MQ=255)                                         cAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTgtt                                                                                                      >  2:563231/1‑145 (MQ=255)                                           aaaTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGc                                                                                                 >  1:491378/1‑150 (MQ=9)                                                                                                         gcgAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAAttt                                        >  2:974018/1‑146 (MQ=255)                                                                                                                       cggtaaccGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCTCCCAAGCGGGAATATCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGca                       <  2:292874/143‑1 (MQ=255)                                                                                                                                           gAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATa  <  2:5982/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                           CGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATA  >  NC_000913/4392304‑4392591

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