Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,828,594	G→A	9.7%	R120H (CGC→CAC)	recX →	regulatory protein RecX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,828,594	0	G	A	9.7%	108.4 / 2.5	41	R120H (CGC→CAC)	recX	regulatory protein RecX
Reads supporting (aligned to +/- strand):  ref base G (18/19);  new base A (2/2);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

GCTCTGGCTATGGCCCTGCGCGTATTCGCGAAGAGCTGTGTCAGCGTGGTTTGGCGCGTGCAGATATCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCTGCTCAGCGGGCGTGATCTCGACGATTATTGAAGCTCGAAACCACAGTGGGTGGCACGTGCCCGTAG  >  NC_002947/1828447‑1828741                                                                                                                                                       |                                                                                                                                                      gCTCTGGCTATGGCCCTGCGCGTATTCGCGAAGAGCTGTGTCAGCGTGGTTTGGCGCGTGCAGATATCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAgcgc                                                                                                                                                      <  4:153090/149‑1 (MQ=255)    ctGGCTATGGCCCTGCGCGTATTCGCGATGAGCTGTGGCAGCGTTGTTTGGCGCGTGCAGATATCGACCAGGCCCTGCGTTAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCg                                                                                                                                                   <  4:43924/149‑1 (MQ=255)       gCTATGGCCCTGCGCGTATTCGCGAAGAGCTGTGTCAGCGTGGTTTGGCTCGTGCAGATATCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCa                                                                                                                                                 <  8:134355/148‑1 (MQ=255)           tGGCCCTGCGCGTATTCGCGAAGAGCTGTGTCAGCGTGGTTTGGCGCGTGCAGATATCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGa                                                                                                                                              <  8:260771/147‑1 (MQ=255)              cccTGCGCGTATTCGCGAAGAGCTGTGTCAGCGTGGTTTGGCGCGTGCAGATATCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCcgc                                                                                                                                         <  6:281693/149‑1 (MQ=255)                       tATTCGCGAAGAGCTGTGTCAGCGTGGTTTGGCGCGTGCAGATATCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCcgcg                                                                                                                               >  7:40404/1‑149 (MQ=255)                             cgAAGAGCTGTGTCAGCGTGGTTTGGCGCGTGCAGATATCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCcag                                                                                                                          >  6:199785/1‑148 (MQ=255)                               aaGAGCTGTGTCAGCGTGGTTTGGCGCGTGCAGATATCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCcag                                                                                                                          >  6:148561/1‑146 (MQ=255)                                agagCTGTGTCAGCGTGGTTTGGCGCGTGCAGATATCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGa                                                                                                                      >  6:48/1‑149 (MQ=255)                                   gCTGTGTCAGCGTGGTTTGGCGCGTGCAGATATCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGAccc                                                                                                                   <  5:48/149‑1 (MQ=255)                                       tgtCAGCGTGGTTTGGCGCGTGCAGATATCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGtt                                                                                                               <  5:199785/149‑1 (MQ=255)                                        gtCAGCGTGGTTTGGCGCGTGCAGATATCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCCAAGATGGTCAGCGCAGGGAGTTTCCGGGCAGCGCCCGCAAGATACGCGC‑CGCCGCGCCCAGCAGCCCCGGt                                                                                                                >  1:105252/1‑149 (MQ=255)                                           aGCGTGGTTTGGCGCGTGCAGATATCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTg                                                                                                           >  6:123767/1‑149 (MQ=255)                                                  tttGGCGCGTGCAGATATCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGAGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTAcc                                                                                                    >  7:247053/1‑149 (MQ=255)                                                     ggCGCGTGCAGATATCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑t                                                                                                 <  2:185291/148‑1 (MQ=255)                                                             cAGATATCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCt                                                                                        <  2:326541/149‑1 (MQ=255)                                                               gATATCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTAt                                                                                      <  5:282319/149‑1 (MQ=255)                                                                atatCGACCAGGCCCTGCGTGAAAGCGGAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTAt                                                                                      >  1:168955/1‑149 (MQ=255)                                                                 tatCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATg                                                                                     <  8:40404/148‑1 (MQ=255)                                                                  atCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCACCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGa                                                                                   >  2:8620/1‑149 (MQ=255)                                                                  atCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCACCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGa                                                                                   >  2:3374/1‑149 (MQ=255)                                                                     gACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACAt                                                                                >  7:37234/1‑149 (MQ=255)                                                                             ccTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGgtc                                                                        >  6:71880/1‑149 (MQ=255)                                                                              cTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCACCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCGCTGGGTTCCCTATGGACATGATCGgtc                                                                        <  1:3374/149‑1 (MQ=255)                                                                              cTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCACCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCGCTGGGTTCCCTATGGACATGATCGgtc                                                                        <  1:8620/149‑1 (MQ=255)                                                                                      aaGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCTGCt                                                                 >  3:236519/1‑147 (MQ=255)                                                                                           gAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCTGCTCAGCgg                                                           <  2:105252/149‑1 (MQ=255)                                                                                                gAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCTGCTCAGCGGGCGt                                                       >  5:260818/1‑148 (MQ=255)                                                                                                                   cgcgACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGCGAGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCGTATGGACAGAAACGGTCGTCTGCTCAGCGGGCGTCATCTCGACGATTATTGat                                    >  1:80661/1‑148 (MQ=255)                                                                                                                   cgcgACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCTGCTCAGCGGGCGTGATCTCGACGATTATTGaa                                    <  6:260818/148‑1 (MQ=255)                                                                                                                    gcgACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCAC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCTGCTCAGCGGGCGTGATCTCGACGATTATTGAAGc                                  >  6:75434/1‑149 (MQ=255)                                                                                                                      gACGTGGGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCTGCTCAGCGGGCGTGATCTCGACGATTATTGAAGCTc                                <  8:247053/149‑1 (MQ=255)                                                                                                                         gtgtGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCTGCTCAGCGGGCGTGATCTCGACGATTATTGAAGCTCGaa                             <  5:71880/149‑1 (MQ=255)                                                                                                                          tgtgGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCTGCTCAGCGGGCGTGATCTcgacga                                            <  7:231438/133‑1 (MQ=255)                                                                                                                          tgtgGCAGCGCAAG‑‑TTTGCCGGGCAGCGACCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCGGCTCAGCGGGCGTGATCTcgacga                                            >  8:231438/1‑133 (MQ=255)                                                                                                                            tgGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCTGCTCAGCGGGCGTGATCTCGACGATTATTGAAGCTCGAAAcc                          <  2:26178/149‑1 (MQ=255)                                                                                                                             ggCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCTGCTCAGCGGGCGTGATCTCGACGATTATTGAAGCTCGAAACca                         >  3:314331/1‑149 (MQ=255)                                                                                                                                  cgcAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCTGCTCAGCGGGCGTGATCTCGACGATTATTGAAGCTCGAAACCACAgtg                    <  5:123767/149‑1 (MQ=255)                                                                                                                                   gcAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCTGCTCAGCGGGCGTGATCTCGACGATTATTGAAGCTCGAAACCACAgtgg                   <  2:168955/149‑1 (MQ=255)                                                                                                                                      aG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCTGCTCAGCGGGCGTGATCTCGACGATTATTGAAGCTCGAAACCACAgtgggtg                >  5:287049/1‑149 (MQ=255)                                                                                                                                      aG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCTGCTCAGCGGGCGTGATCTCGACGATTATTGAAGCTCGAAACCACAgtgggtg                <  4:314331/149‑1 (MQ=255)                                                                                                                                                 ggCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCTGCTCAGCGGGCGTGATCTCGACGATTATTGAAGCTCGAAACCACAGTGGGTGGCACGTGcc       >  5:257310/1‑149 (MQ=255)                                                                                                                                                      cgcCCGCAAGATCCGCAC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCTGCTCAGCGGGCGTGATCTCGACGATTATTGAAGCTCGAAACCACAGTGGGTGACACGTGCCCGTAg  <  5:75434/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                      GCTCTGGCTATGGCCCTGCGCGTATTCGCGAAGAGCTGTGTCAGCGTGGTTTGGCGCGTGCAGATATCGACCAGGCCCTGCGTGAAAGC‑GAGGTGAATTGGAGCGAGCGTATGCGCGACGTGTGGCAGCGCAAG‑‑TTTGCCGGGCAGCGCCCGCAAGATCCGCGC‑AGCCGCGCCCAGCAGACCCGGTTCCTGGCTTACCG‑TGGGTTCCCTATGGACATGATCGGTCGTCTGCTCAGCGGGCGTGATCTCGACGATTATTGAAGCTCGAAACCACAGTGGGTGGCACGTGCCCGTAG  >  NC_002947/1828447‑1828741

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