Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	REL606	728,198	T→C	12.9%	P119P (CCA→CCG)	abrB ←	predicted regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	728,198	0	T	C	12.9%	72.9 / 4.6	32	P119P (CCA→CCG)	abrB	predicted regulator
Reads supporting (aligned to +/- strand):  ref base T (14/13);  new base C (2/2);  total (17/15)
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.42e-01

CTTCAGCGTTATCACCCAGCATCATGCGGGTGACCAGAACAGCGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAGCGCACCAATAACCAACCCACGATGGCGCTGGAAAGCAGCGTCACCAGTAAAATAGCAAGCACGATCGGCCAGTTGGCGGCCAGCGTGGTGAGAATC  >  REL606/728053‑728339                                                                                                                                                  |                                                                                                                                              cTTCAGCGTTATCACCCAGCATCATGCGGGTGACCAGAACAGCGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCCGgag                                                                                                                                           >  2:383408/1‑150 (MQ=255) cTTCAGCGTTATCACCCAGCATCATGCGGGTGACCAGAACAGCGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCCGgag                                                                                                                                           >  2:274653/1‑150 (MQ=255)  tgcagcctgtccACCCGGCAGCATGCGTGTGACCTGGAGAGCGGCCCCCTCGACGCACAGCACTCGCAGCTATTGCACAAGCGCCCCCGGGCGAATACCTGCGCCGTAATCATGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGgaga                                                                                                                                          <  1:140782/139‑1 (MQ=255)              aCCCAGCATCATGCGGGTGACCAGAACAGCGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGc                                                                                                                              >  2:619723/1‑150 (MQ=255)                   gCATCATGCGGGTGACCAGAACAGCGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCgg                                                                                                                         >  2:201690/1‑150 (MQ=255)                      tcatGCGGGTGACCAGAACAGCGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTAt                                                                                                                      >  2:369529/1‑150 (MQ=255)                      tcatGCGGGTGACCAGAACAGCGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCAACGGTAt                                                                                                                      >  1:280273/1‑150 (MQ=255)                             ggTGACCAGAACAGCGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAgg                                                                                                               >  2:249826/1‑150 (MQ=255)                                gACCAGAACAGCGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAg                                                                                                            >  1:507601/1‑150 (MQ=255)                                gACCAGAACAGCGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCCGCGCCGCCTGGGGAAGAACCCCAGGCACCGGTATTTCCAAGCAg                                                                                                            >  1:257197/1‑150 (MQ=255)                                     gAACAGCGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGACTATCCGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCCGGCACCGGCACTTCCAGGCAGTGAGc                                                                                                       >  1:558401/1‑150 (MQ=255)                                          gCGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAg                                                                                                  <  2:560020/150‑1 (MQ=255)                                          gCGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAg                                                                                                  <  1:201111/150‑1 (MQ=255)                                          gCGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAg                                                                                                  <  1:249115/150‑1 (MQ=255)                                          gCGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCCGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAg                                                                                                  <  1:382296/150‑1 (MQ=255)                                          gCGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCCGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAg                                                                                                  <  1:273871/150‑1 (MQ=255)                                          gCGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCAGTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAg                                                                                                  <  2:118816/150‑1 (MQ=255)                                           cGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAgc                                                                                                 <  2:397253/150‑1 (MQ=255)                                               ccccGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAGCGCAc                                                                                             <  2:509049/150‑1 (MQ=255)                                                      acgaACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAGCGCACCAATaac                                                                                      >  2:105658/1‑150 (MQ=255)                                                       cgaACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAGCGCACCAATaacc                                                                                     <  2:281073/150‑1 (MQ=255)                                                       cgaACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAGCGCACCAATaacc                                                                                     <  1:356899/150‑1 (MQ=255)                                                       cgaACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAGCGCACCAATaacc                                                                                     <  2:276389/150‑1 (MQ=255)                                                                  cGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGGAATCTTGCGCCATAGCGACCATGGCCGCCGCGACGCCTGGAGAAGAGCCCCGGGGCACGCTAGTGCCAGGCAGTGAGCTATAGCGCACCAATacacacaccacgcgg                                                                          >  1:254997/1‑136 (MQ=255)                                                                      gATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAGCGCACCAATAACCAAACCACGATGgcgc                                                                      >  1:542302/1‑150 (MQ=255)                                                                              aTAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAGCGCACCAATAACCAACCCACGATGGCGCTGGAAagc                                                              <  2:257931/150‑1 (MQ=255)                                                                                                 tatCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAAACACAGGCACCGGTATTTCCAGGCAGTGAGCTATAGGGCACCAATAACCAACCCACGATGGCGCTGGAAAGCAGCGTCACCAGGAAAAtga                                           >  1:476483/1‑148 (MQ=255)                                                                                                              aTCTTGCGCCATAGCGACCATTGCCGATGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAGCGCACCAATAACCAACCCACGATGGCGCTGGAAAGCAGCGTCACCAATAAAATAGCAAGCACGATCgg                              >  1:359041/1‑150 (MQ=255)                                                                                                              aTCTTGCGCCATAGCGACCATTGCCGATGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAGCGCACCAATAACAAACCCACGATGGCGCTGGAAAGCAGCGTCACCAGGAAAATAGCAAGCACGATCgg                              >  1:369461/1‑150 (MQ=255)                                                                                                                cTTGCGCCATAGCGACAATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAGCGAACCAATAACCAACCCACGATGGCGCTGGAAAGCAGCGTCACCAGTAAAATAGCAAGCACGATCGGcc                            <  1:330001/150‑1 (MQ=255)                                                                                                                      ccATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAGCGCACCAATAACCAACCCACGATGGCGCTGGAAAGCAGCGTCACCAGTAAAATAGCAAGCACGATCGGCCAGTTgg                      <  2:364300/150‑1 (MQ=255)                                                                                                                        aTAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCCGGCCGGGAGATAGAGCGCACCAATAACCAACCCACGATGGCGCTGGAAAGGAGAGTAGCCAGTAAAATAGCAAGCACGAACCGCCAGTTggcg                    >  1:196193/1‑150 (MQ=255)                                                                                                                                          gcggcgCCAGGAGAAGAACCACAGGCACCGGTCTTTCCAGGCGGTGAGCTAAAGCGCACCAATAACCAACCCACGAAGGCGCCGGAAAGCGACGTCACCAGCAAGATAGCAAACAGGGGCGGCCAGTTGCCGGACAGcaacgtggacata  >  1:582332/1‑138 (MQ=255)                                                                                                                                                  |                                                                                                                                              CTTCAGCGTTATCACCCAGCATCATGCGGGTGACCAGAACAGCGGCCCCCGCGACGAACAGCACTCGCAGATATTGCATAAACGCCACCAGGCGAATATCTGCGCCGTAATCTTGCGCCATAGCGACCATTGCCGCTGCGCCGCCTGGAGAAGAACCCCAGGCACCGGTATTTCCAGGCAGTGAGCTATAGCGCACCAATAACCAACCCACGATGGCGCTGGAAAGCAGCGTCACCAGTAAAATAGCAAGCACGATCGGCCAGTTGGCGGCCAGCGTGGTGAGAATC  >  REL606/728053‑728339

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