Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,410,419	G→A	L312L (CTG→TTG)	yfbS ←	putative transporter YfbS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,410,419	0	G	A	96.9%	91.6 / ‑3.5	32	L312L (CTG→TTG)	yfbS	putative transporter YfbS
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (18/13);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-01

GCATACCGATCAGTTTCCAGTTACCCACAACCAGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAGTACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGG  >  NC_000913/2410179‑2410644                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                  gCATACCGATCAGTTTCCAGTTACCCACAACCAGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGc                                                                                                                                                                                                                         <  1:103787/251‑1 (MQ=255)   aTACCGATCAGTTTCCAGTTACCCACAACCAGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTc                                                                                                                                                                                                                       >  2:137178/1‑251 (MQ=255)           cAGTTTCCAGTTACCCACAACCAGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAAtt                                                                                                                                                                                                               <  2:123299/251‑1 (MQ=255)            aGTTTCCAGTTACCCACAACCAGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTg                                                                                                                                                                                                              <  1:256200/251‑1 (MQ=255)            aGTTTCCAGTTACCCACAACCAGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTg                                                                                                                                                                                                              >  1:332027/1‑251 (MQ=255)                tCCAGTTACCCACAACCAGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTc                                                                                                                                                                                                           >  1:146959/1‑250 (MQ=255)                    attACCCACAACCAGGCTGATATCGCCCAGCCGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCCCTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTAc                                                                                                                                                                                                      >  2:22597/2‑251 (MQ=255)                       aCCCACAACCAGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACaa                                                                                                                                                                                                    >  1:225790/1‑250 (MQ=255)                                         atCGCCCAGCAGCAGGGGCTCATCCGCCAACGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACAt                                                                                                                                                                                 >  2:367727/1‑251 (MQ=255)                                                agcagcagAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGc                                                                                                                                                                          >  2:142035/1‑251 (MQ=255)                                                 gcagcagAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCa                                                                                                                                                                         <  1:245427/251‑1 (MQ=255)                                                       gagGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACata                                                                                                                                                                   <  2:375387/251‑1 (MQ=255)                                                        agGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAGTACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACatat                                                                                                                                                                  <  1:188875/251‑1 (MQ=255)                                                           cTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCATGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCa                                                                                                                                                                >  1:354844/1‑250 (MQ=255)                                                                acGCCAGCGAACCTTCCAGCGCCACACCATTGCGCTTTAGCCCCACCACATTCAGTCCGTATCGGGTACGCAAACCCATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAACTTGCCTGAGATCGACAGCAGCCGCAGCCATATCAATa                                                                                                                                                             >  2:214228/2‑248 (MQ=255)                                                                            cTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCAcgc                                                                                                                                              >  2:373407/1‑251 (MQ=255)                                                                                    gccacgccaTTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGaa                                                                                                                                       >  2:232911/1‑250 (MQ=255)                                                                                      cacgccaTTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAAtt                                                                                                                                    <  2:337095/251‑1 (MQ=255)                                                                                                            cccacATTCATTCCGTAGCTGGTACGAAAACCAATTTCGCGCACCTATTTTCCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGAGCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACTCGCGCGAAATTCAGAAACCCCATTAGCGTTCAc                                                                                                              <  2:390996/250‑1 (MQ=255)                                                                                                                   tCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAAcg                                                                                                       >  2:222488/1‑251 (MQ=255)                                                                                                                             gCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCATAAACCCCATTAACGTTCACGATAACGCGACGAAAAc                                                                                             <  1:123042/251‑1 (MQ=255)                                                                                                                                ggTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAAcgcc                                                                                          >  1:6750/1‑251 (MQ=255)                                                                                                                                           cAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAACTTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTc                                                                                >  2:342735/1‑250 (MQ=255)                                                                                                                                                                 aaTCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCAt                                                         <  1:361558/251‑1 (MQ=255)                                                                                                                                                                              tCAGGAATTAATGAAATCTCTGCCATGCCCCCCTCAAGGGCCTGGTAAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTa                                             >  1:258404/1‑250 (MQ=255)                                                                                                                                                                                           aaagctctGCCATTACCCCATCAAGGACCTGGTCAGAAAAATACTCGCCACGCAAGACCATCGGCTCCAGCAATGGCTCACTACAAAATTTCCGGAGAGCGACATAAGCCTCAGACATATCAATAAGCAAAACGTCACGCGCTCGAAATTCATAAACCCCACAAACGTTCACGATAACGCGACGAAAACGCCGCCAGCTTTCAACTCCGATGACGTTAGCGCCATAACGCGCACGTAATTTGAGACCATcc                               <  2:249388/247‑1 (MQ=255)                                                                                                                                                                                            aaTCTCTGCCATGCCCACATCAAGGGACTAGTCAGAAACATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACAACGATGACGTTAGCGCCATAACGCTCACGCAATTCGAGATCATcc                               >  2:172094/1‑250 (MQ=255)                                                                                                                                                                                              tctcTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATAACATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATACGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGc                            >  1:67950/1‑251 (MQ=255)                                                                                                                                                                                                 ctGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGcc                           <  1:137178/249‑1 (MQ=255)                                                                                                                                                                                                     cATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGAc                     >  1:181581/1‑251 (MQ=255)                                                                                                                                                                                                                    ggCCTGGTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATc      <  1:203477/251‑1 (MQ=255)                                                                                                                                                                                                                         ggTCAGAAAAATACTCGCCGCGCAATACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCggg  <  2:337043/250‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                  GCATACCGATCAGTTTCCAGTTACCCACAACCAGGATGATATCGCCCAGCAGCAGAGGCTCATCCGCCAGCGAACCTTCCAGCGCCACGCCATTGCGCTTTAGCCCCACCACATTCAGTCCGTAGCGGGTACGAAAACCAATTTCGCGCACCGATTTACCAATCAGTTCTGACTCAGGAATTAATGAAATCTCTGCCATGCCCACATCAAGGGCCTGGTCAGAAAAATACTCGCCGCGCAGTACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGG  >  NC_000913/2410179‑2410644

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