Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,247,134	A→C	Y467D (TAC→GAC)	lysP ←	lysine:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,247,134	0	A	C	96.0%	77.8 / NA	25	Y467D (TAC→GAC)	lysP	lysine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (15/9);  total (16/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATCAGACACAGAATGAATGCGAAGATCGGCCCCAGTGGGAAGAAACCTGAACGGTACGGCAGATCGTTAATGTCGTGTCCCTGCAATACGTAACCGCGACGGAAGCGATAGTGGCTAATGGCAATCCCCAGCCAGG  >  NC_000913/2246895‑2247379                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      cGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCaa                                                                                                                                                                                                                                             >  2:106626/1‑250 (MQ=255) cGGAACAGCCTCTGATTTGTCAGCATTTTTTGAGTGCGGCGTAACGATAATGCAAACGATAACAATTATCAAGCCAAATGTTAAATGTTATATCCGCATGATTTGTTAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTTATTTCCTTTAATCAGCTTGTCTCCGAACCaaa                                                                                                                                                                                                                                            <  1:85407/251‑1 (MQ=255)                      aGCCTATTTTGATTGCGGCGTAACGATAATGCAACCGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATAc                                                                                                                                                                                                                       >  2:296163/1‑251 (MQ=255)                                                      aaaCGATAACAATTATCAATCCAAATGGTAAATTTTATATTCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCGTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGAGACCAATAGACGTTGCCGCGACGCCGCCCCAGTCa                                                                                                                                                                                       <  2:120211/251‑1 (MQ=255)                                                                         tCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACg                                                                                                                                                                    >  2:77767/1‑251 (MQ=255)                                                                                   aaaTTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGtt                                                                                                                                                          <  1:106626/251‑1 (MQ=255)                                                                                     atgttATAGCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTGCTTATCGTTCTGCGGGAACTGCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGGTGCCGCTATGCCGCCCCAGTCAATAGTATCTGGCAGGAACGCTTCGTAGTTCt                                                                                                                                                        <  2:307259/248‑1 (MQ=255)                                                                                                      tttGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATc                                                                                                                                       >  1:187056/1‑251 (MQ=255)                                                                                                                          gAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCCTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATCAGACACAGAATGAATGCGa                                                                                                                    >  2:258862/1‑250 (MQ=255)                                                                                                                               aGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATCAGACACAGAATGAATGCGa                                                                                                                    >  2:180230/1‑245 (MQ=255)                                                                                                                                                gggCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCCGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATCAGACACAGAATGAATGCGAAGATCGGCCCCAGTGGgaagaa                                                                                              >  1:159967/1‑250 (MQ=255)                                                                                                                                                  gCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGCAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATCAGACACAGAATGAATGCGAAGATCGGCCCCAGTGGGAAGAAAcc                                                                                           >  2:296739/1‑251 (MQ=255)                                                                                                                                                      aGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGTATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCttt                                                                                                                                                                            <  2:143552/166‑1 (MQ=255)                                                                                                                                                      aGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCttt                                                                                                                                                                            >  1:143552/1‑166 (MQ=255)                                                                                                                                                                        cTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATCAGACACAGAATGAATGCGAAGATCGGCCCCAGTGGGAAGAAACCTGAACGGTACGGCAGATCGtt                                                                      >  2:170847/1‑250 (MQ=255)                                                                                                                                                                                  aTCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATCAGACACAGAATGAATGCGAAGATCGGCCCCAGTGGGAAGAAACCTGAACGGTACGGCAGATCGTTAATGTCGTGTc                                                           <  1:77767/251‑1 (MQ=255)                                                                                                                                                                                              gAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCACATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATCAGACACAGAATGAATGCGAAGATCGGCCCCAGTGGGACGAAACCTGAACGTTACGGCAGATCGTTAATGTCGTGTCCCTGCAATCCGt                                               >  2:303933/1‑251 (MQ=255)                                                                                                                                                                                                     aTTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATCAGACACAGAATGAATGCGAAGATCGGCCCCAGTGGGAAGAAACCTGAACGGTACGGCAGATCGTTAATGTCGTGTCCCTGCAATACGTAACcgcg                                        <  1:60619/251‑1 (MQ=255)                                                                                                                                                                                                        tCGCTGTAGAGTACGAAGTGAGTTCCTTTACTCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCACCCGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATCAGACACAGAATGAATACGAAGTTCGGCCCCAGTGGTAAGAAACCTGAACGGTACGGCAGATCGTTAATGTCGTGTCCCTGCAATACGTAACCGCGACg                                     >  1:240957/1‑251 (MQ=255)                                                                                                                                                                                                           cTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATCAGACACAGAATGAATGCGAAGATCGGCCCCAGTGGGAAGAAACCTGAACGGTACGGCAGATCGTTAATGTCGTGTCCCTGCAATACGTAACCGCGACGGa                                   >  1:149259/1‑250 (MQ=255)                                                                                                                                                                                                                                ccTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGTGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATCAGACACAGAATGAATGCGAAGATCGGCCCCAGTGGGAAGAAACCTGAACGGTACGGCAGATCGTTAATGTCGTGTCCCTGCAATACGTAACCGCGACGGAAGCGATAGTGGCTAATGGCAAt             <  2:325510/251‑1 (MQ=255)                                                                                                                                                                                                                                    tAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTc                                                                                                                                                                >  1:53165/1‑100 (MQ=255)                                                                                                                                                                                                                                    tAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTc                                                                                                                                                                <  2:53165/100‑1 (MQ=255)                                                                                                                                                                                                                                        cAGCTTGTCGCCGAACCAACTAATCAGGAACAGCGGGATACCAATCTACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAACGTGATAATCAGACACAGAATGAATGCGAAGATCGGCCCCAGTGGGAAGAAACCTGAACGGTACGGCAGATCGTTAATGTCGTGTCCCTGCAATACGTAACCGCGCCGGCAGCGATAGTGGCTAATGTCAATACccggcc     >  2:51792/1‑251 (MQ=255)                                                                                                                                                                                                                                           cTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGATACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTACGTAGTTCTGGCCCAAAGTGATAATCAGACACAGAATGAATGCGAAGATCGGCCCCAGTGGGAAGACACCTGAACGGTACGGCAGATCGTTAATGTAGTGTCCCTGCAATACGTAACCGCGACGGACGCGATAGTGGCTAATGGCACTCCCCAGCCAgg  >  1:53033/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      CGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATCAGACACAGAATGAATGCGAAGATCGGCCCCAGTGGGAAGAAACCTGAACGGTACGGCAGATCGTTAATGTCGTGTCCCTGCAATACGTAACCGCGACGGAAGCGATAGTGGCTAATGGCAATCCCCAGCCAGG  >  NC_000913/2246895‑2247379

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