Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,041,398	T→G	L29R (CTT→CGT)	ygfZ →	folate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,041,398	0	T	G	96.8%	78.6 / NA	32	L29R (CTT→CGT)	ygfZ	folate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (18/12);  total (19/13)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCT  >  NC_000913/3041261‑3041536                                                                                                                                          |                                                                                                                                           gCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGATGCTGATGACGCTTGATGACTGGGCGCGTGcca                                                                                                                                       >  1:113665/1‑143 (MQ=255) gCTCCTTTACCATTAGTCATTAATGCGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACGTTTTCCTCCCCGTAATCGAACTGCGTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGcca                                                                                                                                       <  1:490139/143‑1 (MQ=255)       ttACCATTAGTCATGAATCCGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTa                                                                                                                                 <  2:94067/143‑1 (MQ=255)                 tCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTACCCGTCAGCCGACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGGTGACTGGGCGCGTGGCGCGGTGACTGGGGCGGa                                                                                                                       >  1:272764/1‑143 (MQ=255)                   aTTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCCGGTGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACa                                                                                                                     >  1:522422/1‑143 (MQ=255)                          cGCCGTTAAGCACCTCCGGGCTCTATTATGGCTTTTACACCTTTCCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGaaaa                                                                                                              <  2:524497/143‑1 (MQ=255)                              gTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCATACGGCGTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAatat                                                                                                          >  1:490248/1‑143 (MQ=255)                                          aGGACTCTATTATGGCTTTTACACCTTTTCCTCCCAGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGACGCGTGCCACCATTGCTGACGCGGCACGCGTCAAAAATATGCACGGTCaa                                                                                              >  2:429903/1‑142 (MQ=255)                                            gACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGCCAGCGACAACTATACGCAGGGTCACGt                                                                                            >  2:393346/1‑143 (MQ=255)                                              ctctATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGa                                                                                          <  1:537264/143‑1 (MQ=255)                                               tctATTATGGCTTTTACACCTTTTCCTCCCTGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGCCGCTTGCTGACTGTGCGCGTGCCACCATTACTGGCGCGGACAGCGCAAAATATATGCAGGGTCAGGGGAc                                                                                         >  1:272637/1‑143 (MQ=255)                                               tctATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGGTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGGGGAAAAGATATGCGGGGTCAGGTGAc                                                                                         >  1:354394/1‑143 (MQ=255)                                                      tGGCTTTTACGCCTGTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGAGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGAt                                                                                  <  1:438268/143‑1 (MQ=255)                                                        gCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGt                                                                                >  2:54951/1‑143 (MQ=255)                                                        gCTTTTACACCTCTCCCTCCCCGTCAGCCTACGGCGTCCGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGt                                                                                <  1:192327/143‑1 (MQ=255)                                                          ttttACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCa                                                                              >  1:196544/1‑143 (MQ=255)                                                                   ttttCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAggagt                                                                     >  1:17357/1‑139 (MQ=255)                                                                          ccccGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGGCGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCagaaga                                                              <  2:196554/143‑1 (MQ=255)                                                                           cccGTCAGCTTACGTCTTCTGCCCGTTTACCGCTGCCGCTGATGACGCTTGATGCCTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCGGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGAt                                                             <  2:173955/143‑1 (MQ=255)                                                                                   ccTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGCTGACCGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGCCCGGTCACCGCAGAGGTCAGCCAACGGGCAGAAGAGGAGCAGCt                                                     >  2:437804/1‑143 (MQ=255)                                                                                     tACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTAGTGGATGGGACAGCTAAACCCACAGGTTCAGTGGTGTGACAGCAGATGTCAGCCAGATGGCAGCAGATCAGAACCTGc                                                   >  2:277963/1‑143 (MQ=255)                                                                                          cTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCCTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCCGGTGACCGCAGCTGTCAGCCACATGGCAGCAGATCAGCACCtccccccc                                              >  2:423649/1‑136 (MQ=255)                                                                                                     ttACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGa                                   <  1:54951/143‑1 (MQ=255)                                                                                                                   ggtgaCGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTaaaa                     <  1:31879/141‑1 (MQ=255)                                                                                                                   gatgaCCCTTGCTGACTGGCCGCTTGCCACCATTACTGGCGCAGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCTACGCCAAAGGTaaaa                     <  2:34122/143‑1 (MQ=255)                                                                                                                    atgaCGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGATCAGGTGACAGCAGATGTACGCCAGATGGCAGAAGATCGGCACCTGCTCGCCGCACATCGCGACGCCAAAGGCaaacg                    >  2:472910/1‑141 (MQ=255)                                                                                                                     tgaCGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAAtg                   <  2:490248/143‑1 (MQ=255)                                                                                                                         gCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAACTGTGGa               >  1:123274/1‑143 (MQ=255)                                                                                                                                   tGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACg     >  1:392792/1‑143 (MQ=255)                                                                                                                                   tGGGCGCGTGCCACCATTACTGGCGCGGACAAAGCACGGTCCGTGCAGGGTAAGCAGACAGCAGATGACAGCCAGATGGCACAAAATCAGCACCTGCTCGCCGCCCAGCGCAACGCCAAAGGTAAAAAGTGGAGCAAATTCCg     >  2:447439/1‑143 (MQ=255)                                                                                                                                     ggCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCTGATGTCAGCCGGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATGTACGTc   >  1:89501/1‑143 (MQ=255)                                                                                                                                      gcgcGTGCCACCAGTAGTGGCGCGGACAGCGAAAAATATGTGCTGTGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAATATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCt  <  2:247318/143‑1 (MQ=255)                                                                                                                                          |                                                                                                                                           GCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCT  >  NC_000913/3041261‑3041536

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