Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,041,632	T→A	V107E (GTG→GAG)	ygfZ →	folate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,041,632	0	T	A	96.5%	83.5 / NA	29	V107E (GTG→GAG)	ygfZ	folate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (11/17);  total (11/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCT  >  NC_000913/3041492‑3041767                                                                                                                                             |                                                                                                                                        cgccCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCGTGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCAGCAAGGGCAGATACAATATGCGGTATTCTCTAAAgaga                                                                                                                                       >  1:851087/1‑143 (MQ=255)      tttGGGACTCCAAAGGTAAGATGGGGAGCGATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATc                                                                                                                                  <  2:258083/142‑1 (MQ=255)                     gTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGGGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAg                                                                                                                   >  1:765173/1‑143 (MQ=255)                        aaaTGTGGGGCGATTTACGTCTGTGCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCgt                                                                                                                <  2:896633/143‑1 (MQ=255)                            gttgAGCCATTTACGTCTGTTCCGCGACGGCGGTTGCTTTGCATTGGTTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGtgc                                                                                                            <  1:852595/140‑1 (MQ=255)                                  cAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTg                                                                                                      >  1:583117/1‑143 (MQ=255)                                             tGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGt                                                                                           >  1:468681/1‑143 (MQ=255)                                               ttCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGttt                                                                                         <  1:776960/143‑1 (MQ=255)                                                 cagcgACGGACATGGCTTTTGATGGATTGAACGGCGCCGCGAGAGTGAACCGCAGCTGACTGAACTGAAAAAATATCCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTc                                                                                       <  1:79703/141‑1 (MQ=255)                                                      aCGGCGAGGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGgcg                                                                                  <  1:358796/143‑1 (MQ=255)                                                             tgtcTTTGCATGGATTGAACGGCGCACCGTTCTTGAACCGCATCTGACTTAACTTAAAAAATATGCGGTAACCGCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGcgcgccg                                                                           <  1:177113/140‑1 (MQ=255)                                                                      atggatggAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCaa                                                                  <  1:269468/143‑1 (MQ=255)                                                                           ttGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAAtctc                                                             >  1:398387/1‑143 (MQ=255)                                                                            cgaccGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATctct                                                            <  2:35793/139‑1 (MQ=255)                                                                             gAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCtt                                                           <  2:471486/143‑1 (MQ=255)                                                                             gAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCtt                                                           <  2:583117/143‑1 (MQ=255)                                                                              aaCGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCCCttt                                                          >  2:407197/1‑143 (MQ=255)                                                                                           gcgtgAACCGCAGCTGACTGAACTGACAAAATATGCGGTACTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCtt                                             <  2:467334/143‑1 (MQ=255)                                                                                            ggagAACCAGAGCTGACGGAACTGAAAAAATATGCGGGATTCCCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTc                                            <  1:187497/140‑1 (MQ=255)                                                                                               gAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaa                                         >  1:732877/1‑143 (MQ=255)                                                                                                aacaggaGCTGACTTAGCTGATAAAATATTCGGTAGTGTCTAATGTGTCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTTCCGTTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaa                                        <  2:779589/137‑1 (MQ=255)                                                                                                   cGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGGTGGCAAATCTCTTTAGCGGACGGCCTTCgacagac                                     >  2:254269/1‑142 (MQ=255)                                                                                                    gCAGCTGACTGAACCGACAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGCGCTGCCTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaagaaa                                    <  1:668830/143‑1 (MQ=255)                                                                                                     cAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGaaaa                                   <  1:540514/143‑1 (MQ=255)                                                                                                      aGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAGTCTCTTTAGCGAACTGCCCTCcagagaaaga                                  >  1:774263/1‑141 (MQ=255)                                                                                                      aGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGACGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAACGCACTGCCTTCGAAAGacaaa                                  >  2:303359/1‑143 (MQ=255)                                                                                                         tgactgaAAGTAAAAAATAAGCGCTATTCTCTAAAGAGAACATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTTCCGGTTTTCAGGCGCGCGCCGCGCTGGAAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAg                               <  2:765971/143‑1 (MQ=255)                                                                                                                                    tctcTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTg    <  2:810669/143‑1 (MQ=255)                                                                                                                                      tctAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAGCGAAAAACAGGTAGTCAGACAAGGCCAGACCACTGTGcc  >  2:10402/1‑142 (MQ=255)                                                                                                                                      tctAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGCGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCt  >  1:46861/1‑143 (MQ=255)                                                                                                                                             |                                                                                                                                        CGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCT  >  NC_000913/3041492‑3041767

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑