Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,818,930	G→T	A19A (GCG→GCT)	yicG →	PF03458 family inner membrane protein YicG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,818,930	0	G	T	96.9%	104.5 / NA	32	A19A (GCG→GCT)	yicG	PF03458 family inner membrane protein YicG
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (20/11);  total (20/12)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑CTGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCCCTCCCCTTAAAAAATCTCAATCGTGACAATGCGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGCTCGATGCGCTCGGACTGGTGGTCTTT  >  NC_000913/3818693‑3819173                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    tttAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCgggc                                                                                                                                                                                                                                        >  1:120856/1‑249 (MQ=255)  ttAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCgggcg                                                                                                                                                                                                                                       >  2:379568/1‑248 (MQ=255)            ttttCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATg                                                                                                                                                                                                                             >  2:232032/1‑251 (MQ=255)               tCTTATGCATTTTCTCAGAAAGGAAGACCGAAAAATGCATCTACT‑CTACCGAACATAGCGATTTCCCCTAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACACATCGCTACCCTGCCAGACAGATTTTTAGGTAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGTCGGCCGGTCGTCGCCGCATGGAt                                                                                                                                                                                                                          >  2:424412/1‑251 (MQ=255)                 ttATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATAc                                                                                                                                                                                                                        >  2:219269/1‑251 (MQ=255)                     gCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACAttt                                                                                                                                                                                                                    >  1:45289/1‑251 (MQ=255)                                       agTCAGAAGAATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCg                                                                                                                                                                                                  >  2:126295/1‑251 (MQ=255)                                           agaagaATGCATCT‑ACGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCg                                                                                                                                                                                              >  1:88481/1‑251 (MQ=255)                                               gaaTGCATCT‑ACGCAACAGAAACTAGCGATTTCACATAACTACAAGATATCAATTTCCTCTCACCTTAAAAAATCTCAATCGAGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCAGGGCGACGCCGCATGGATACCTTTGGAGTGATTATTATTGCGACAGCCAc                                                                                                                                                                                          >  2:383755/1‑251 (MQ=255)                                                               aCAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAg                                                                                                                                                                            <  2:45289/250‑1 (MQ=255)                                                                   aaaaTAGCGATTTCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCa                                                                                                                                                                       <  1:379568/251‑1 (MQ=255)                                                                               tCACATATCTACAAGTTATCAATTTCCTCTCACCTGAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTGGGCATTACTGCCGAAGCCAGGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCaa                                                                                                                                                                                     <  2:140608/225‑1 (MQ=255)                                                                               tCACATAACTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCaa                                                                                                                                                                                     >  1:140608/1‑225 (MQ=255)                                                                                     aaCTACAAGTTATCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTctgct                                                                                                                                                      <  2:120856/250‑1 (MQ=255)                                                                                                aTCAATTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATcc                                                                                                                                          >  1:243762/1‑251 (MQ=255)                                                                                                    aTTTCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTc                                                                                                                                      >  1:436836/1‑251 (MQ=255)                                                                                                      ttCCTCTCACCTTAAAAAATCTCAATCGTGACAATACGCACAATTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGATGTTACACATTTTGTATATGGTTGGCATTACTGCCGAAGCCATGACGGGAGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGACACCGCAATTGTCGGCGGGTCGGTGCTGGTTATTCTTCTGGGCACCTATACGCTCgg                                                                                                                                    >  2:141833/1‑251 (MQ=255)                                                                                                          cctcacctTAAAAAATCGCAACCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTTCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCTCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGAGATTCTGCTGGGCCACTATCCGCTCGGCTgg                                                                                                                                <  2:415977/251‑1 (MQ=255)                                                                                                            cctcctTAAAAAATCTCAACCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTCCTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGGCGTCGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGt                                                                                                                              <  2:126550/248‑1 (MQ=255)                                                                                                             cacctTAAAAAATCTCAAGCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTgg                                                                                                                                <  1:219269/246‑1 (MQ=255)                                                                                                              accttccaaaaTCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCCCCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTc                                                                                                                             >  2:182965/8‑250 (MQ=255)                                                                                                                   aaaaaaTCTCAATCGTGACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAacac                                                                                                                       <  1:232032/251‑1 (MQ=255)                                                                                                                                cGTGACAATACGCACAAATCGCTACCCTGCCAGACCGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATAtgtgc                                                                                                          >  1:102221/1‑250 (MQ=255)                                                                                                                                   gACAATACGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATATGTGAtt                                                                                                        <  1:18705/250‑1 (MQ=255)                                                                                                                                           gcACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCt                                                                                               >  2:324878/1‑251 (MQ=255)                                                                                                                                                              cagacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGCATATGTGATTATCGTCGCTACCGCCGCAGTACTcaccc                                                                            >  2:88530/1‑250 (MQ=255)                                                                                                                                                                                    cacagggtgatACACATTTGGTATTGGGTTGGCATTCCGGCCGAAGCCATGGGGGGGGCTCTGGCGGCCGGGCGACGCCGCAGGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATAGGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGcc                                                      <  2:102221/241‑1 (MQ=255)                                                                                                                                                                                             ttACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTgcg                                             >  2:111470/1‑251 (MQ=255)                                                                                                                                                                                              tacacaaattGTATTTGAATGGCATTACTGCCGATGCCATGACGCGGGCTATGGAGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCAACGCAAATGGCGGAGGGTCAGTGCGCGATATTCTGCTGAGCAACTATCCGCTCGGCTGGGTCCACCACCCGGCATATGTGCTTATCGTCGCTACCGCCGCAGTACTCACCCCTATCGTTGATCCTGTAATGCCTTACCTgcgc                                            >  2:196929/9‑251 (MQ=255)                                                                                                                                                                                                                              gAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCAc                                                                                                                                               >  1:137523/1‑120 (MQ=255)                                                                                                                                                                                                                              gAAGCCATGACGGGGGCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCAc                                                                                                                                               <  2:137523/120‑1 (MQ=255)                                                                                                                                                                                                                                        ggggggCTCTGGCGGCCGGGCGACGCCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTAGCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGCTCGATGCGCTTGGACTGGTGGTCttt  <  1:182965/250‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    TTTAACCTATTTTTTCTTATGCATTTTCTCAGACAAGAAGTCAGAAGAATGCATCT‑CTGCTACAGAAAATAGCGATTTCACATAACTACAAGTTATCAATTTCCCCTCCCCTTAAAAAATCTCAATCGTGACAATGCGCACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGCTCGATGCGCTCGGACTGGTGGTCTTT  >  NC_000913/3818693‑3819173

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