Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	3,125,061	A→G	7.8%	F194S (TTT→TCT)	ytnA ←	putative amino acid permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	3,125,061	0	A	G	7.8%	146.6 / 4.7	51	F194S (TTT→TCT)	ytnA	putative amino acid permease
Reads supporting (aligned to +/- strand):  ref base A (23/24);  new base G (2/2);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

CTGTATCGATTGCTTTGGCTAACGATTTTTGCGGGTTTTTTACTTCGCCAGCCGTCACCCCGATCATTTCAATCCCTAAATAGGCGAACATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCAATAATCATCAATAAACCGACAGCAATCATGGATAAGATTGCAACAATTTTAATCAATGCGAACCAAAATTCCAGC  >  BS168_reference_genome/3124924‑3125198                                                                                                                                          |                                                                                                                                          cTGTATCGATTGCTTTGGCTAACGATTTTTGCGGGTTTTTTACTTCGCCAGCCGTCACCCCGATCATTTCAATCCCTAAATAGGCGAACATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGaaaa                                                                                                                                           <  2:20076/138‑1 (MQ=255)    tATCGATTGCTTTGGCTAACGATTTTTGCGGGTTTTTTACTTCGCCAGCCGTCACCCCGATCATTTCAATCCCTAAATAGGCGAACATGACCATTTGAAGAGACAATAAACCCCCTTTTAATCCATGAGGGAAAAAtcc                                                                                                                                       >  4:327497/1‑139 (MQ=255)    tATCGATTGCTTTGGCTAACGATTTTTGCGGGTTTTTTACTTCGCCAGCCGTCACCCCGATCATTTCAATCCCTAAATAGGCGAACATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAAtc                                                                                                                                        >  6:47221/1‑138 (MQ=255)    tATCGATTGCTTTGGCTAACGATTTTTGCGGGTTTTTTACTTCGCCAGCCGTCACCCCGATCATTTCAATACCTAAATAGGCGAACATGACCATTTGAAGAGACAATAAAACCCCTTTAAATCCATGAGGGAAAAAtc                                                                                                                                        >  5:259774/1‑138 (MQ=255)                tGGCTAACGATTTTTGCGGGTTTTTTACTTCGCCAGCCGTCACCCCGATCATTTCAATCCCTAAATAGGCGAACATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCca                                                                                                                           >  8:186244/1‑139 (MQ=255)                   cTAACGATTTTTGCGGGTTTTTTACTTCGCCAGCCGTCACCCCGATCATTTCAATCCCTAAATAGGCGAACATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACaa                                                                                                                        <  8:399023/139‑1 (MQ=255)                              tGCGGGTTTTTTACTTCGCCAGCCGTCACCCCGATCATTTCAATCCCTAAATACGCGAACATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGTAAAAATCCTCCGTTGTTCCACAAATTGCTGAtt                                                                                                              <  6:304427/138‑1 (MQ=255)                                       tttACTTCGCCAGCCGTGAACCCGATCATTTCAAACCCTAAATACGCGAACATGACCATTTGAAGACACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCa                                                                                                     <  6:173637/138‑1 (MQ=255)                                           cTTCGCCAGCCGTCACCCCGATCAATTAAATCCCTAAATAGGCGAACATGACCATTTGAAGACACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGc                                                                                                 <  2:20472/138‑1 (MQ=255)                                               gCCAGCCGTCACCCCGATCATATCAATCCCTAACAAGGCGAACATGACCATTTGAAGAGACAATACAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCa                                                                                            <  6:210775/139‑1 (MQ=255)                                                       tCACCCCGATCATTTCAATCCCTAAATAGGCGAACATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGAc                                                                                    <  3:145307/139‑1 (MQ=255)                                                       tCACCCCGATCATTTCAATCCCTAAATAGGCGAACATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCACATTGCCGAc                                                                                    <  8:189053/139‑1 (MQ=255)                                                       tCACCCCGATCATTTCAATCCCTAAATAGGAGAACATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGAc                                                                                    <  7:393479/139‑1 (MQ=255)                                                       gcACCCCGATCATTTCAATCCCTAAATAGGCGAACATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGAc                                                                                    <  6:102469/138‑1 (MQ=255)                                                               aTCATTTCAATCCCTAAATAGGCGAACATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGCCTCCTGCa                                                                             >  1:161857/1‑138 (MQ=255)                                                               aTCATTTCAATCCCTAAATAGGCGAACATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCaa                                                                            >  6:297546/1‑139 (MQ=255)                                                               aTCATTTCAATCCCTAAATAGGCGAACATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCaa                                                                            >  6:224163/1‑139 (MQ=255)                                                               aTCATTTCAATCCCTAAATAGGCGAACATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCa                                                                             >  2:187045/1‑138 (MQ=255)                                                                   tttCAATCCCTAAATAGGCGAACATGACCAGTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGGTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCaataat                                                                        <  6:17763/139‑1 (MQ=255)                                                                            cTAAATAGGCGAACATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCaataat                                                                        <  7:196116/130‑1 (MQ=255)                                                                            cTAAATAGGCGAACATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCaataat                                                                        >  8:196116/1‑130 (MQ=255)                                                                                    gCGAACATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCAATAATCATCAATAAACCGACAg                                                       >  8:163277/1‑139 (MQ=255)                                                                                       aaCATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAAt                                                                                                   >  7:219442/1‑92 (MQ=255)                                                                                       aaCATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAAt                                                                                                   <  8:219442/92‑1 (MQ=255)                                                                                         cATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCAATAATCATCAATAAACCGACAGCAATc                                                  <  3:81054/139‑1 (MQ=255)                                                                                         cATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCAATAATCATCAATAAACCGACAGCAATc                                                  <  5:47221/139‑1 (MQ=255)                                                                                         aaTGACCATTTGAAGAGACTATAAAACCCCTTTTAATCCATGAGGGAAAATTCCTCCGTTGTTCCACAAATTGCTGATTCCGTTTGCAATGCCCCCATTGCCGACTCCTGCAATAATCATCAATAAACCGACAGCAATc                                                  <  6:259774/138‑1 (MQ=255)                                                                                             aCCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCAATAATCATCAAt                                                                 <  4:336978/120‑1 (MQ=255)                                                                                             aCCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCAATAATCATCAAt                                                                 >  3:336978/1‑120 (MQ=255)                                                                                                 tttGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTcc                                                                                 >  5:32522/1‑100 (MQ=255)                                                                                                 tttGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTcc                                                                                 <  6:32522/100‑1 (MQ=255)                                                                                                 tttGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCAATAATCATCAATAACCCGACAGCAATCATGGAt                                            >  4:258578/1‑137 (MQ=255)                                                                                                      agagACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGcc                                                                                       >  5:3137/1‑89 (MQ=255)                                                                                                      agagACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGcc                                                                                       <  6:3137/89‑1 (MQ=255)                                                                                                         gaCAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACt                                                                                   <  5:280276/90‑1 (MQ=255)                                                                                                         gaCAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACt                                                                                   >  6:280276/1‑90 (MQ=255)                                                                                                         gaCAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCAATAATCATCAATAAACCGACAGCAATCATGGATAAGATTGcaa                                  <  3:125281/139‑1 (MQ=255)                                                                                                         gaCAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCAATAATCATCAATAAACCGACAGCAATCATGGATAAGATTGcaa                                  <  1:187045/139‑1 (MQ=255)                                                                                                         gaCAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCAATAATCATCAATAAACCGACAGCAATCATGGATAAGATTGcaa                                  <  8:94310/139‑1 (MQ=255)                                                                                                          aCAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCa                                                                                            >  5:400618/1‑80 (MQ=255)                                                                                                          aCAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCa                                                                                            <  6:400618/80‑1 (MQ=255)                                                                                                                      cTTTTAATCCATGAGGGAAAGATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGAc                                                                                    <  2:364795/76‑1 (MQ=255)                                                                                                                      cTTTTAATCCATGAGGGAAAGATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGAc                                                                                    <  2:359327/76‑1 (MQ=255)                                                                                                                      cTTTTAATCCATGAGGGAAAGATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGAc                                                                                    >  1:364795/1‑76 (MQ=255)                                                                                                                      cTTTTAATCCATGAGGGAAAGATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGAc                                                                                    >  1:359327/1‑76 (MQ=255)                                                                                                                      cTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCTATGCCCCCATTGCCGCCTCCTGCAATAATCATCAATAAACCGACAGCAATCATGGATAAGATTGCAACAATTTTAATCaa                     >  7:321775/1‑139 (MQ=255)                                                                                                                      cTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCAATAATCATCAATAAACCGACAGCAATCATGGATAAGATTGCAACAATTTTAATCaa                     >  2:261217/1‑139 (MQ=255)                                                                                                                      cTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCAATAATCATCAATAAACCGACAGCAATCATGGATAAGATTGCAACAATTTTAATCaa                     >  1:341393/1‑139 (MQ=255)                                                                                                                            aTCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGAAATAATCATCAATAAACCGACAGCAATCATGGATAAGATTGCAACAATTTTAATCAATGCGa                <  1:163227/138‑1 (MQ=255)                                                                                                                                aTGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCAATAATCATCAATAAACCGACAGCAATCATGGATAAGATTGCAACAATTTTAATCAATGCGAACCa            >  3:127198/1‑138 (MQ=255)                                                                                                                                 tGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCAATAATCATCAATAAACCGACAGCAATCATGGATAAGATTGCAACAATTTTAATCAATGCGAACCaa           >  5:304022/1‑138 (MQ=255)                                                                                                                                 tGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCAAGAATCATCAATAAACCGACAGCAATCATGGATAAGATTGCAACAATTTTAATCAATGCGAACCaa           >  1:201990/1‑138 (MQ=255)                                                                                                                                         aaaTCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCAATAATCATCAATAAACCGACAGCAATCATGGATAAGATTGCAACAATTTTAATCAATGCGAACCAAAATTCCAGc  <  7:163277/139‑1 (MQ=255)                                                                                                                                          |                                                                                                                                          CTGTATCGATTGCTTTGGCTAACGATTTTTGCGGGTTTTTTACTTCGCCAGCCGTCACCCCGATCATTTCAATCCCTAAATAGGCGAACATGACCATTTGAAGAGACAATAAAACCCCTTTTAATCCATGAGGGAAAAATCCTCCGTTGTTCCACAAATTGCTGATTCCGGTTGCAATGCCCCCATTGCCGACTCCTGCAATAATCATCAATAAACCGACAGCAATCATGGATAAGATTGCAACAATTTTAATCAATGCGAACCAAAATTCCAGC  >  BS168_reference_genome/3124924‑3125198

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