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.9%	90.0 / NA	32	L29R (CTT→CGT)	ygfZ	folate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (16/15);  total (16/16)
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

TACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTC  >  NC_000913/3041268‑3041535                                                                                                                                   |                                                                                                                                          tACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGTCTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTCCCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTAc                                                                                                                               <  2:141551/143‑1 (MQ=255)       tAGTCATTAATACGCCGGTAATCAACTCAGGGCGCTATCGTGGCTTTGACACCTTTTCCTCACCATCGGCGGAGGTCTGCTGCCCGTTTACCTCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGgcgc                                                                                                                         <  2:683511/143‑1 (MQ=255)          tCATTAATAAGCCGTTAAGCAACTCAGGACGCTACTATGGCGTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGa                                                                                                                      <  1:354130/143‑1 (MQ=255)            aTTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCGCGTTTACCGCTGACGCTGGTGACGCTTGGTGACTGGGCGCGTGACAACATTCCTGGCGCGGACa                                                                                                                    >  1:45029/1‑143 (MQ=255)            aTTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATCACTGGGCACGTGCCACCACTCCTGGCGCGGAcc                                                                                                                    >  1:385427/1‑142 (MQ=255)                          aaGCAACTCAGGACTCTATTATGGCTGTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTTACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATg                                                                                                      <  1:364964/143‑1 (MQ=255)                            gCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCACCCTCCAGCTTCTGCCCGTTTACCTCCGACGCTGATGACGCTTGATGACTGGGCCCGTCCCACCATGCCTGACCCGAACAGAGACAAttatcgcca                                                                                                    >  2:707276/1‑138 (MQ=255)                                     ggctctATTATGGCTTTTACACCTTTCCCCCCCCCGCCGCCTCCGGCTTCTGCCCTTTTACCGCTGACGCCGATGCGGCTTGATGACTGGGCGCGTGCCACTATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGt                                                                                           <  2:709039/141‑1 (MQ=255)                                     gcctctATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCCCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGt                                                                                           <  1:536787/141‑1 (MQ=255)                                          tattatGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGAcag                                                                                      >  2:112732/1‑143 (MQ=255)                                                 gCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTACCCGTTTACTGCTGACGCTGATGACGCTTGTGGACCGGGCCGGTGCCACCATTACTGGCGCGGAGAGCGAAACAAAAGTGCCCGGTCGGGTGACGGCCGaggc                                                                               >  2:131376/1‑140 (MQ=255)                                                     ttACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGc                                                                           <  2:712203/143‑1 (MQ=255)                                                         acCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAAGATATGCCGGGTCAGTTGACAGCCGCTGTCAGCCAgg                                                                       >  1:120866/1‑142 (MQ=255)                                                            tttcCCTCCCCGTCAGCCTACGGCTTCTGCCCGTGTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCAGCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATgg                                                                    <  2:357647/143‑1 (MQ=255)                                                               gcctccCCGTCAGCCTACGGCTTCGGCCCGTCGACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCag                                                                 <  1:422818/142‑1 (MQ=255)                                                                  tccCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGGCAGGTGACAGCCGATGTCAGCCAGATGGCagaag                                                              >  1:325406/1‑143 (MQ=255)                                                                   ccccGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATAGCagaaga                                                             >  1:424614/1‑143 (MQ=255)                                                                     ccGTCAGCCTACGGCCTCGGCCCGTTTACCGCTGACGCGGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGACAGAAGATc                                                           <  2:137887/143‑1 (MQ=255)                                                                        tCAGCCTACGGTTTCTGCCCGTTTACCGCTGACGCGGAAGACTCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCTGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGc                                                        <  2:70489/143‑1 (MQ=255)                                                                                      cttctcGTTTCCCGCTGGCTCTTATTACTTTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCgccgcc                                          <  1:302720/138‑1 (MQ=255)                                                                                         cccGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCAt                                       >  1:465336/1‑143 (MQ=255)                                                                                              ctgccGCGGACGCTGATGCCTCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGGCAGCTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGa                                  <  2:248189/140‑1 (MQ=255)                                                                                               tACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATAAGCAACTGCTCGCCGCCCATTGCGAc                                 >  1:182895/1‑143 (MQ=255)                                                                                                  ggctggccctgaGGACGCTGGAGAACTGGGCGCTTGCCACCATTACGGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGcc                              <  1:384023/142‑1 (MQ=255)                                                                                                             atgaCGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTaaaag                   >  1:267286/1‑142 (MQ=255)                                                                                                             atgaCGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGAAAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGGCCATTGCGACGCCAAAGGTAAAAt                   >  1:511097/1‑143 (MQ=255)                                                                                                                aCGGTTGATGACTGGGCGCGTGCAACCATTACTGGCGCGAACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAAtgtg                <  2:328115/143‑1 (MQ=255)                                                                                                                   cTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAAGATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAg             >  1:374238/1‑143 (MQ=255)                                                                                                                   cTTGATGACTGGGCGCGTGCCACCAGTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAg             <  2:656037/143‑1 (MQ=255)                                                                                                                            tGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACg    >  1:73299/1‑143 (MQ=255)                                                                                                                            tGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCACCGCCCACGGCGACGCCAAAGGACAAATGAGGCGCAATGTACg    >  2:327499/1‑143 (MQ=255)                                                                                                                              ggCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTACAATGTGTCGCCAtttccatc  >  1:389909/1‑138 (MQ=255)                                                                                                                                   |                                                                                                                                          TACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTC  >  NC_000913/3041268‑3041535

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: ‑