Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	6,850	A→G	Y370Y (TAT→TAC)	yaaJ ←	putative transporter YaaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	6,850	0	A	G	94.7%	54.1 / NA	19	Y370Y (TAT→TAC)	yaaJ	putative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (13/5);  minor base T (0/1);  total (13/6)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGC  >  NC_000913/6608‑7078                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     gggCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGt                                                                                                                                                                                                                               >  1:263495/1‑250 (MQ=255)                                        gCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACcac                                                                                                                                                                                       <  1:33471/251‑1 (MQ=255)                                                              agagCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGAAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTAAGCAccc                                                                                                                                                                 >  2:169806/1‑251 (MQ=255)                                                                               gTCAAATGGGTAATAGCCAGGCAGGCCATTATGAGATCTGCCAGTTGCCACATCAGCGGAAGACGTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAacac                                                                                                                                                <  2:13825/251‑1 (MQ=255)                                                                                               ccATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTgg                                                                                                                                >  1:222525/1‑251 (MQ=255)                                                                                                         cATTATGATATCTGACAGTTGCCACATCAGCGGAAGACTAAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGCTTTCGGCGTAAATGGAGTTGGCAACGATGGCGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGa                                                                                                                       >  1:53750/1‑250 (MQ=255)                                                                                                            tATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATAc                                                                                                                   >  1:38053/1‑251 (MQ=255)                                                                                                              tGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGCACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCt                                                                                                                 >  1:114388/1‑251 (MQ=255)                                                                                                                   gcTGCCAGTTGCCACCTCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGAACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTcc                                                                                                              >  2:154318/2‑249 (MQ=255)                                                                                                                                                                           gAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAAt                                                                                                                     <  2:48520/186‑1 (MQ=255)                                                                                                                                                                           gAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAAt                                                                                                                     >  1:48520/1‑186 (MQ=255)                                                                                                                                                                                         gCAAACACAAGATCGTTTTGGGTTTTTTCAGGCAAAAAAAGAAGAGATTGTTTTCGGGTAAAATGTAGTGGACAACGATGGAGCTGAAGGCAAACAAAATAACCACGAGGTTAAAAACCCCAGCACCCCAGGAACCCATAAACCCCCCCATCGCCTTCGGGATAAGCTGAATCCCTTCCAGCGGCATGTGGTTTGGGCGTTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGgtg                                      <  1:169806/251‑1 (MQ=255)                                                                                                                                                                                             acacAAGATCGATTTCGGGGTGTACAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTCAATGTAGTTGGCACCGATGGAGCTGAAGGCAAACAGAATACACACGAGGGTAACAAACTCAGCACCCCAGGAGCACCTTAACACACGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGATGGTGCCGTTACCCGCCAATAATATCAGCATGTCGCTTGCCGTCCAGATGACCAggcggtcgc                                  >  2:244551/1‑244 (MQ=255)                                                                                                                                                                                                             gggTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCAAGGCGCTTGCCGTACAGATGACCCGGGTGTCGATAAAAATGCCAatcat                  >  2:162903/1‑251 (MQ=255)                                                                                                                                                                                                                             aaaaGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTCCAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTgc  >  1:11383/1‑251 (MQ=255)                                                                                                                                                                                                                                   agagATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGcccccc                                                                                                                                                                >  2:46659/1‑87 (MQ=255)                                                                                                                                                                                                                                   agagATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCAcccc                                                                                                                                                                <  1:46659/87‑1 (MQ=255)                                                                                                                                                                                                                                    gagaTTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGAGCCGTTACCCGCCAg                                                                              >  1:175523/1‑168 (MQ=255)                                                                                                                                                                                                                                    gagaTTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGGTGTGCCGTAACCCGCCAg                                                                              <  2:175523/168‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     GGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGC  >  NC_000913/6608‑7078

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