Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,161,093	C→T	S9S (AGC→AGT)	fabR →	DNA‑binding transcriptional repressor FabR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,161,093	0	C	T	96.8%	88.5 / ‑5.8	31	S9S (AGC→AGT)	fabR	DNA‑binding transcriptional repressor FabR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/15);  minor base A (0/1);  total (15/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 = 9.66e-01

CAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTT‑ACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTG  >  NC_000913/4160855‑4161303                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                    cAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTT‑ACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGc                                                                                                                                                                                                         <  1:215808/251‑1 (MQ=255)             gTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTT‑ACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑ATTCAGAc                                                                                                                                                                                            <  1:227859/251‑1 (MQ=255)                              atCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGGTTTATTGCGTT‑ACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑AGTGAGACATCGTGATGGGCGTAag                                                                                                                                                                           <  2:407547/251‑1 (MQ=255)                                   gAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTT‑ACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGTCAAGG‑ATTCAGACATCGTGATGGGCGTAAGAgcg                                                                                                                                                                       <  1:123193/250‑1 (MQ=255)                                                      tGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTATCCGTTTACGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTT‑ACCGTTCATTCACAATACTGCAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGAAGCAAAAAGAAAAAACccg                                                                                                                                                   >  2:309125/1‑251 (MQ=255)                                                         ttCTTACTTTTGCCCTTACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTT‑ACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgcc                                                                                                                                                 <  1:220232/250‑1 (MQ=255)                                                                                   tGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTT‑ACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGcc                                                                                                                      >  1:244471/1‑251 (MQ=255)                                                                                   tGGATTTTTATCCACGATACTCTCCGTTTAAGCGGCAGGTTCCCGCTGAACGTAAAATACCCGGCCAAAAAAGTGAAAAAAATATTTTTTTTTGCGTT‑CCCTTTATTTCACAATACGGGAGCAACCCATTA‑GTTCCTTTCTCTGGTATAGTCCCAGTATTATTTTGGGAAAGA‑ATTCAGACATCGTTGTGGGCGTAAGAGCTCATCAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGc                                                                                                                       <  2:188963/250‑1 (MQ=255)                                                                                            aTCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTT‑ACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTg                                                                                                             >  1:152770/1‑251 (MQ=255)                                                                                                  gaTACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTT‑ACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAAc                                                                                                       >  1:188496/1‑251 (MQ=255)                                                                                                  gaTACTCTCCGTTTAAGCGGAAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTT‑ACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAAc                                                                                                       >  1:150006/1‑251 (MQ=255)                                                                                                  aaTACTCACCTTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGACAAAAGAATGGAAGTAATTATGTTTTATT‑CGTTTCCCGTTCATCCCAAATCCGGGAGAAATCCAGTATGTTA‑ATTATCTGGTTTAGTGCAATTATTACTTTTGGCAAGG‑ATTCAGACATCGTGAGGGGAGTAAGACCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAAc                                                                                                       <  2:215970/250‑1 (MQ=255)                                                                                                    tACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTT‑ACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGc                                                                                                     >  2:71668/1‑251 (MQ=255)                                                                                                                aaGCGGCAGGTTCCCGCTGTGCGGAAAAAACCCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTT‑ACCGTTCATTCACAATACGGGAGCAACCCAGTATGTGC‑ATACTCTGGTATAGTGCAGGAAGTCATTTTGGCAAGT‑TTTCAGATATCGTGATGGGCGAAAGATCGCAAGAAAAAAAAAAAAGCCGCCGTTCTCTGGTGGAAGCCTCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAg                                                                                          <  1:386709/250‑1 (MQ=255)                                                                                                                      cAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTT‑ACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGt                                                                                   >  2:127777/1‑251 (MQ=255)                                                                                                                                   aactaaaaaGACCCGGCCAAAGAATTGCTGTAAATTTTTTTTTTTGCTTT‑CCCGTTCATTCACAATACGGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTATTACTTTTGGAAAGG‑ATTCAGACATCGTTATGGGCGAAAGAGCGCAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAGGTgg                                                                      <  2:325530/246‑1 (MQ=255)                                                                                                                                     agaaaaaGAACCGGCCAAAGAATTGAAGTAAATATTTTTTATGGCGTT‑ACCGTTCATTCACAATACTGGAGCAATCAAGTTTGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑ATTCAGACATCGTGAGGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGgcg                                                                    <  2:188496/248‑1 (MQ=255)                                                                                                                                       tAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTT‑ACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGcgcg                                                                  >  1:126990/1‑251 (MQ=255)                                                                                                                                                cGGCAAAAAAATTGCAATAAATATGTTTTATTGCGTT‑CCCTTTCATTCACAATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTATTATTTTTGGCAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCggg                                                         <  1:127777/251‑1 (MQ=255)                                                                                                                                                              cAGTAAATATGTTTTATTGCGTT‑ACCGTCCATTCACAATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑ATTCAGACATCGCGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGACGCGGGCATTGCTCCCACct                                           >  1:318447/1‑251 (MQ=255)                                                                                                                                                                 aaaaTATGTTTTATTGCGTT‑ACCGTTCATTCACAATACTGGGGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGAAAGG‑ATTCAGACATCGTGATGGGCGAAAGAGCGCAAAAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCtt                                        <  2:266587/250‑1 (MQ=255)                                                                                                                                                                    attttttttATTGCGTT‑CCCTTTAATTCACAATACTGGAGCATCCCAGTATTTTC‑ATTCTCTGGTATAGTACCAGTAGTCCTTTTGGAAAGG‑ATTAAGACTTCGGGATGGGCGAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCttttt                                     <  1:87766/246‑1 (MQ=255)                                                                                                                                                                             aTTGCGTT‑ACCTTTCATGCACAATACTGGAGAAATCCAGTATTTTC‑ATTCGCTGTTAGAGTCCGAGTAGTACTTT‑GGAAAGGGATTAAGACATCGTGAGGGGCGTAAGAGCGTAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCAtt                            <  1:263964/251‑1 (MQ=255)                                                                                                                                                                              ttGCGTT‑ACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCAttt                           >  1:109502/1‑251 (MQ=255)                                                                                                                                                                                      ccccttcattcaCAATCCTGGACCAATCCATTTTTTTA‑ATACTCTGGATTAGTGCGAGTATTTCTTTGGGCAAGG‑TTTCAGACGTCGTGAGGGGGGTAAGCGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGAc                    <  2:297031/247‑1 (MQ=255)                                                                                                                                                                                       ccGTTCATTCAAAATACTGGAGCAATCCATTATGTTC‑ATTCTCTGTTATAGTGCCAGTAGTACTTTGGGAAAGG‑ATTCAAACATCGTGATGGGCGTAAGAGCGAAACAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACg                   <  1:372730/251‑1 (MQ=255)                                                                                                                                                                                         gTTCATTCACAATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGt                  >  2:83072/1‑250 (MQ=255)                                                                                                                                                                                                  caATACTGGATCAATCCAGACTGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑ATTCAGACATCGTGATGCGCGTAAGATCGCAACAAAAAGACAAAACCCGCCGTTCGCTGGTGGCAGCCGCATTTAGCCAATTAAGTGCTGAACGCCGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCCCCTCTTTTTATCGTCATTTCCGCTACGTAGACGAACTg        >  2:353842/1‑251 (MQ=255)                                                                                                                                                                                                  caATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTg        >  2:200042/1‑251 (MQ=255)                                                                                                                                                                                                    aTACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTggg      >  2:243378/1‑251 (MQ=255)                                                                                                                                                                                                        tGGAGCAATCCAGTATGTTA‑ATTCTCTGGTATAGTGCCAGTAGTACTTTTGGCAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTg  >  1:150810/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                    CAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTT‑ACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTC‑ATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGG‑ATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTG  >  NC_000913/4160855‑4161303

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