Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	191,720	G→A	13.1%	intergenic (+12/‑135)	tsf → / → pyrH	translation elongation factor EF‑Ts/uridylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	191,720	0	G	A	13.1%	104.5 / 6.1	38	intergenic (+12/‑135)	tsf/pyrH	translation elongation factor EF‑Ts/uridylate kinase
Reads supporting (aligned to +/- strand):  ref base G (17/16);  new base A (2/3);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTGTTGGTCAGCTGCTGAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATG  >  NC_000913/191576‑191867                                                                                                                                                 |                                                                                                                                                     aaCTGTTGGTCAGCTGCTGAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAgcc                                                                                                                                                  <  2:344699/149‑1 (MQ=255)         gTCAGCTGCTGAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAgg                                                                                                                                         <  2:37709/149‑1 (MQ=255)                 ctgAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGcttc                                                                                                                                 >  1:274523/1‑149 (MQ=255)                  tgAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGcttct                                                                                                                                >  1:133628/1‑149 (MQ=255)                   gAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGcttctt                                                                                                                               >  2:170926/1‑149 (MQ=255)                           aTAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGcc                                                                                                                       >  4:341936/1‑149 (MQ=255)                                gCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCt                                                                                                                  >  2:296046/1‑149 (MQ=255)                                gCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGAAGCCG‑CCCTCAGGCGGctcctt                                                                                                                               <  3:108156/136‑1 (MQ=255)                                              cTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTa                                                                                                    >  2:421673/1‑149 (MQ=255)                                               ttCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAg                                                                                                       >  3:255517/1‑145 (MQ=255)                                               ttCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAg                                                                                                       <  4:255517/145‑1 (MQ=255)                                                    ccGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCtt                                                                                              <  3:341936/149‑1 (MQ=255)                                                      gCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGAAGCCG‑CCCTCAGGCGGctcctt                                                                                                                               >  4:458941/1‑114 (MQ=38)                                                      gCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGAAGCCG‑CCCTCAGGCGGctcctt                                                                                                                               <  3:458941/114‑1 (MQ=38)                                                                     aaGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGaaa                                                                             >  3:351066/1‑149 (MQ=255)                                                                      aGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGaaaa                                                                            <  1:170926/149‑1 (MQ=255)                                                                      aGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGaaaa                                                                            <  3:276417/149‑1 (MQ=255)                                                                               gaAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCtttt                                                                                                                             <  3:12344/91‑1 (MQ=255)                                                                               gaAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCtttt                                                                                                                             >  4:12344/1‑91 (MQ=255)                                                                                   gTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATg                                                               <  2:133628/149‑1 (MQ=255)                                                                                        gactgactTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCt                                                          <  1:296046/149‑1 (MQ=255)                                                                                           tgactTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCtt                                                                                                                 >  4:324823/1‑91 (MQ=255)                                                                                           tgactTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCtt                                                                                                                 <  3:324823/91‑1 (MQ=255)                                                                                               tttGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTa                                                   >  4:182099/1‑149 (MQ=255)                                                                                                    agcagAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGc                                                                                       >  4:301489/1‑108 (MQ=255)                                                                                                    agcagAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGc                                                                                       <  3:301489/108‑1 (MQ=255)                                                                                                                     tGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTATCCCTTGTGGGGCTGCGCt                                                                                 <  2:267572/97‑1 (MQ=255)                                                                                                                     tGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCt                                                                                 >  1:267572/1‑97 (MQ=255)                                                                                                                           aGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAgg                       >  4:260841/1‑149 (MQ=255)                                                                                                                                tCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaaga                  >  1:324698/1‑149 (MQ=255)                                                                                                                                   tAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAAc               <  3:260841/149‑1 (MQ=255)                                                                                                                                        aTCAAAAAGGCGCCG‑CCTGAGGGCGGCTTCTTTTTGTTCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGtt                                   >  3:54813/1‑124 (MQ=255)                                                                                                                                        aTCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGtt                                   <  4:54813/124‑1 (MQ=255)                                                                                                                                         acAAAAAGAAGCCGCCCTCA‑GGCGGCTCCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCAtt                                            >  4:3958/2‑114 (MQ=37)                                                                                                                                         acAAAAAGAAGCCG‑CCCTCAGGCGGCTCCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCAtt                                            <  3:3958/113‑1 (MQ=37)                                                                                                                                            aaaaGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTAc       <  1:109030/148‑1 (MQ=255)                                                                                                                                            aaaaGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTAc       >  2:109030/1‑148 (MQ=255)                                                                                                                                               aGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAAt   <  4:1177/149‑1 (MQ=255)                                                                                                                                                ggAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAAt   <  1:256764/148‑1 (MQ=255)                                                                                                                                                ggAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAAt   >  2:256764/1‑148 (MQ=255)                                                                                                                                                ggAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATg  >  2:82646/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                     AACTGTTGGTCAGCTGCTGAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATG  >  NC_000913/191576‑191867

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