Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,041,634	A→C	T108P (ACC→CCC)	ygfZ →	folate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,041,634	0	A	C	93.8%	68.3 / NA	33	T108P (ACC→CCC)	ygfZ	folate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (21/9);  total (23/10)
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.96e-01

TGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTG‑AACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAG  >  NC_000913/3041518‑3041748                                                                                                                      |                                                                                                                   tgtgCAGCAATTTACTTCTGTGCCGCGACGGCGATGGCTTTGCCTGGATTGAACGGCGCCGCGTGCGTG‑AACCGCAGCCGACTGAACTGAAAAAATATTCGGTATTCTCTAAAGTGTCCATCGCGCCAg                                                                                                        <  2:284287/129‑1 (MQ=255)  gtgGAGCAATTTACGTCTGTTCCGCGACGGCGATAGCTTTGCATGGATTGAACGGCGCAGCGTGCGTG‑AACCGCAGCTGACTGAACTGAAAAAATATTCGGTATTCTCTAAAGTGCACATAGCCCCaca                                                                                                       >  2:61423/1‑127 (MQ=255)    ggaggaATTTCCGTCTGTTCCGCGACGGCGATGGCTTTGCTTGGATTGAACGGCGCAGCGTGCGTG‑ATCTGCAGCTGATTGATTTTTTAATATATGCGGTATGCTCTAAAGGGCCCATCGCGCCAgacg                                                                                                     <  1:325604/129‑1 (MQ=255)         aaTTTACGTCTGTTCCGCGACGGCAATGGCATTACATGGAATAAACGACGCAGCGTGCGTG‑AACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAg                                                                                                >  1:481971/1‑129 (MQ=255)           tttACGTCTGTTCCGCGACGGCGATGGCTTTGCAAAAATTGAACGGCGCAGAATACATG‑AACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCg                                                                                              >  1:419849/1‑129 (MQ=255)            ttACGTCTGTTCCGCGTCGGCGATGGCTTTGCATAAATTGAACGGCGCAGCGTGCGTG‑AACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCgt                                                                                             >  1:354040/1‑129 (MQ=255)            ttACGTCTGTTCCGCGACGGCGATGGCTTTGCATGAATTAAACAGAGAAACATGCGTG‑AACAGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCgt                                                                                             >  1:11312/1‑129 (MQ=255)                    gTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTG‑AACCGAAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTgctgct                                                                                     >  1:128509/1‑129 (MQ=255)                    gTTCCGCGACGGCGATGGCTTTGCATGGATTGAACAAAACAGCGAACGTA‑AACAGCAGCTGACAGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTgctgct                                                                                     >  1:20695/1‑129 (MQ=255)                     ttacgcgACGGCGATGGCTTTGCATGGATTGAACGAC‑CACGCTGAGAGTAATCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGagctgcgt                                                                                    >  2:115441/4‑127 (MQ=255)                            aCGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTG‑AACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTACAGTGCCCATCGCGCCCGACGACGACCGTGTGCTGCTTGGCGTTg                                                                             >  2:66439/1‑129 (MQ=255)                                     gCTTTGCATGGATTGAACAGCGCAGCGTGCGTG‑AACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCGGGtatta                                                                    >  2:303446/1‑128 (MQ=255)                                      cTTTGCATGGATTGAACGGCGCAGCGTGCGTG‑AACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCa                                                                   >  2:241447/1‑129 (MQ=255)                                       tttGCATGGATTGAACGGCGCAGCGTGCGTG‑AACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATTGCTTCATTCGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAg                                                                  <  1:227047/129‑1 (MQ=255)                                          gCATGGATTGAACGGCGCAGCGTGCGTG‑AACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGGGTTGCCGGTTTTCAGgcg                                                               >  2:436442/1‑129 (MQ=255)                                               gATTGAACGGCGCAGCGTGCGTG‑AACCGCAGCTGACTGAACTGAAAAAATATACGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCACgcgcgcgc                                                          >  2:164564/1‑129 (MQ=255)                                               gATTGAACGGCGCAGCGTGCGTG‑AACCGCAACTAACAGAAAAAAAAAAATAAGCGGAAAACACTAAAGAGCCCAACGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGgcgcgcgg                                                          >  1:23092/1‑128 (MQ=255)                                                aTTGAACGGCGCAGCGTGCGTG‑AACCGCAGCTGACTGAACTAAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGcgcgcc                                                         >  1:255459/1‑129 (MQ=255)                                                      cGGCGCAGCGTGCGTG‑AACCCCAGCTGACTGAACTAAAAAAATAAGCGATATTCTCTAAAATAACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTg                                                   >  1:405121/1‑129 (MQ=255)                                                       ggCGCAGCGTGCGTG‑AACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATAGCGCCAGACAACGAGCGTGTGCTGCTTGGTGTTGCCGGCTTCCAGGCGCGCGCCGCGCCgg                                                  >  2:529306/1‑129 (MQ=255)                                                          gcAGCGTGCGTG‑ACCCGCAGCTGACTGAACTGAAAAAATATGCGGTATACTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGCGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCaa                                               >  1:187344/1‑129 (MQ=255)                                                            aGCGTGCGTG‑AACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCTCATCGCGTCATACTATTTTCTTTTTCTGCTTGTTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAAt                                             <  1:230030/129‑1 (MQ=255)                                                            aGCGTGCGTG‑AACCGCAGCTGACTGAACTGAAAAAATATGCGGAATTCTCTAAAGAAACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCCGGCGCGCGCCGCGCTGGCAAAt                                             >  1:427006/1‑129 (MQ=255)                                                                 gcgtg‑AACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTTTTTGTTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCtt                                        <  1:633916/129‑1 (MQ=255)                                                                 gcgtc‑AACCGCCGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGTCGACGAGTGTGTGCTGCTTTGTTTTGCTGGTTTTCGGGCTAGCGCCGCGCTGGCAAAGCTCtt                                        <  1:424568/129‑1 (MQ=255)                                                                           gCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGTGTGTTCTTCTTGTTTTTGTCGGTTTTCAGGAGCGCGCCGCGCAGGCAAAGCTCTTTAGCGAACt                               <  1:281871/129‑1 (MQ=255)                                                                            cAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCTTTTGTGCGTGACTACAAGCGTGTGCTGCTTGGTTTTGCTGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTg                              <  1:457164/129‑1 (MQ=255)                                                                                  actgaactgaAAAAATATGCGGTATTCTCTAAAGTGCCCATAGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGACGGGCCGCAGGGGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCg                        >  1:124383/1‑129 (MQ=255)                                                                                    tgaactgaaAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGTTGGTTTTTAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaa                      <  1:241447/129‑1 (MQ=255)                                                                                    tgaactgaaAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCATTAGCGAACTGCCTACgaa                      >  2:508679/1‑129 (MQ=255)                                                                                    tgaactgaaAAAATATGCGGTATTCTCTAAAGTGCCCATCGCACCAGACGACGAGCTTGAACTGATTGGGGTTCCCGGCTTTCAGGCGCGCGCCACGCTGGCCAATCGCTCCACCGGACTGCCATccac                      >  2:230851/1‑126 (MQ=255)                                                                                                atagGCGGTCTTCGCTAGAGGGCCCATCCCGCCTGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGt          <  2:552546/129‑1 (MQ=255)                                                                                                        tATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCATAAGCTAAATGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAaag  >  1:547616/1‑129 (MQ=255)                                                                                                                      |                                                                                                                   TGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTG‑AACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAG  >  NC_000913/3041518‑3041748

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