Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,405,326	A→G	T196A (ACA→GCA)	abgR →	putative DNA‑binding transcriptional regulator of abgABT operon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,405,326	0	A	G	94.9%	122.5 / ‑0.8	39	T196A (ACA→GCA)	abgR	putative DNA‑binding transcriptional regulator of abgABT operon
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (19/18);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.54e-01

TTCACCATCAATACCTATTATCAGGGACCGTACGACCACGAATTTACTTTTGAGAAATTACTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGACAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCGAGACGTTCTCAGCCTGTATCAGTCTGGTGGCAAAAAGCGATTTTCTCAGCAAAC  >  NC_000913/1405179‑1405470                                                                                                                                                    |                                                                                                                                                 ttCACCATCAATACCTATTATCAGGGACCGTACGACCACGAATTTACTTTTGAGAAATTACTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCa                                                                                                                                                <  2:287381/150‑1 (MQ=255)             aCCTATTATCAGGGACCGTACGACCACGAATTTACTTTTGAGAAATTACTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCcac                                                                                                                                   >  1:305811/1‑151 (MQ=255)              ccTATTATCAGGGACCGTACGACCACGAATTTACTTTTGAGAAATTACTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCcaca                                                                                                                                  <  1:107197/151‑1 (MQ=255)              ccTATTATCAGGGACCGTACGACCACGAATTTACTTTTGAGAAATTACTAGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCcac                                                                                                                                   >  2:153855/1‑150 (MQ=255)                    atCAGGGACCGTACGACCACGAATTTACTTTTGAGAAATTACTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGc                                                                                                                              >  1:168149/1‑149 (MQ=255)                      cAGGGACCGTACGACCACGAATTTACTTTTGAGAAATTACTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGct                                                                                                                          >  2:348688/1‑151 (MQ=255)                         ggACCGTACGACCACGAATTTACTTTTGAGAAATTACTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGctact                                                                                                                       <  2:45128/151‑1 (MQ=255)                         ggACCGTACGACCACGAATTTACTTTTGAGAAATTACTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGACAATGCCGACGCCACACGGCAGctact                                                                                                                       <  1:388867/151‑1 (MQ=255)                                aCGACCACGAATTTACTTTTGAGAAATTACTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAAc                                                                                                                 <  2:515759/150‑1 (MQ=255)                                 cGACCACGAATTTACTTTTGAGAAATTACTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAAc                                                                                                                 >  1:192035/1‑149 (MQ=255)                                    ccACGAATTTACTTTTGAGAAATTACTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGc                                                                                                                              <  1:83089/133‑1 (MQ=255)                                    ccACGAATTTACTTTTGAGAAATTACTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGc                                                                                                                              >  2:83089/1‑133 (MQ=255)                                    ccACGAATTTACTTTTGAGAAATTACTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTg                                                                                                            <  1:403836/151‑1 (MQ=255)                                          aTTTACTTTTGAGAAATTACTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGaa                                                                                                      <  2:125004/151‑1 (MQ=255)                                           tttACTTTTGAGAAATTACTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGaa                                                                                                      <  2:197213/150‑1 (MQ=255)                                              aCTTTTGAGAAATTACTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGc                                                                                                  >  1:496043/1‑151 (MQ=255)                                                         aTTACTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACCCCCCCCCATTGGTGCCGGTTCGATCAAACAGTAACTGAATAACAGCGGGGCAATGCGGAGGCCACAGGGCAGCTACTAAAAACGGTTGAGGGAATTGCTTgccaaccgg                                                                                       >  2:557155/1‑143 (MQ=255)                                                            aCTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGgcg                                                                                    >  1:539132/1‑151 (MQ=255)                                                                 aaaaGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAAc                                                                               <  1:156797/151‑1 (MQ=255)                                                                        aTTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCcac                                                                          >  2:120460/1‑149 (MQ=255)                                                                               aTCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGgtg                                                                 >  1:369606/1‑151 (MQ=255)                                                                                   tttGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGACAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGgt                                                                  <  1:519233/146‑1 (MQ=255)                                                                                         gccCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCg                                                        >  1:400661/1‑150 (MQ=255)                                                                                           ccGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCg                                                        >  1:463569/1‑148 (MQ=255)                                                                                             gggACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCGAGACg                                                   <  1:89468/151‑1 (MQ=255)                                                                                                 cacCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCGAGACGTtct                                               >  1:503134/1‑151 (MQ=255)                                                                                                    cccGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCGAGACGTtctc                                              >  2:534303/1‑149 (MQ=255)                                                                                                          aTTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCGAGACGTTCTCAGCCTGTa                                      <  1:595430/151‑1 (MQ=255)                                                                                                            tGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCGAGACGTTCTCAGCCTGTAt                                     <  2:496033/150‑1 (MQ=255)                                                                                                                       cGATCAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCGAGACGTTCTCAGCCTGTATCAGTCtggtgg                          <  1:555630/150‑1 (MQ=255)                                                                                                                           cAAACAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCGAGACGTTCTCAGCCTGTATCAGTCTGGTGGCaaa                      >  2:51837/1‑150 (MQ=255)                                                                                                                             aaCAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCGAGACGTTCTCAGCCTGTATCAGTCTGGTGGCAAAAAg                   <  1:591962/151‑1 (MQ=255)                                                                                                                               cAGTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCGAGACGTTCTCAGCCTGTATCAGTCTGGTGGCAAAAAGCg                 <  2:112111/151‑1 (MQ=255)                                                                                                                                 gTTACTGGATTTCAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACTATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCGAGACGTTCTCAGCCTGTATCAGTCTGGTGGCAAAAAGCGAt               <  2:486097/151‑1 (MQ=255)                                                                                                                                 gTTACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCGAGACGTTCTCAGCCTGTATCAGTCTGGTGGCAAAAAGCGAt               <  2:332158/151‑1 (MQ=255)                                                                                                                                  ttACTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCGAGACGTTCTCAGCCTGTATCAGTCTGGTGGCAAAAAGCGAtt              >  1:84955/1‑151 (MQ=255)                                                                                                                                     cTGGATTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCGAGACGTTCTCAGCCTGTATCAGTCTGGTGGCAAAAAGCGATTTtc           >  1:501749/1‑151 (MQ=255)                                                                                                                                         aTTACAGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCGAGACGTTCTCAGCCTGTATCAGTCTGGTGGCAAAAAGCGATTTTCTCAg       <  1:582083/151‑1 (MQ=255)                                                                                                                                              aGCTGGGCAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCGAGACGTTCTCAGCCTGTATCAGTCTGGTGGCAAAAAGCGATTTTCTCAGCAAAc  >  2:174895/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                 TTCACCATCAATACCTATTATCAGGGACCGTACGACCACGAATTTACTTTTGAGAAATTACTGGAAAAGCAATTCGCGATCTTTTGCCGCCCGGGACACCCCGCCATTGGTGCCCGTTCGATCAAACAGTTACTGGATTACAGCTGGACAATGCCGACGCCACACGGCAGCTACTACAAACAGTTGAGTGAATTGCTTGACGATCAGGCGCAAACGCCACAGGTCGGTGTAGTCTGCGAGACGTTCTCAGCCTGTATCAGTCTGGTGGCAAAAAGCGATTTTCTCAGCAAAC  >  NC_000913/1405179‑1405470

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