Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,145,846	A→C	intergenic (+177/‑18)	dnaG → / → rpoD	DNA primase/RNA polymerase sigma factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,145,846	0	A	C	100.0%	172.3 / NA	41	intergenic (+177/‑18)	dnaG/rpoD	DNA primase/RNA polymerase sigma factor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (17/24);  total (17/24)

AAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGC  >  REL606/3145773‑3145918                                                                          |                                                                         aaaCAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAg                                                                        <  1:981071/76‑1 (MQ=255) aaaCAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAg                                                                        >  2:546935/1‑76 (MQ=255)  atcaacCGTATAATATCAGCGTTCACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAgt                                                                       <  2:1557531/74‑1 (MQ=255)        cGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGAt                                                                 >  1:382245/1‑76 (MQ=255)          ttGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCACTCGCCCAACACCAACCTCATGAAATCAGTGTGGATAc                                                               <  1:1473612/76‑1 (MQ=255)           tgtaTTTCAGCGTTAGCGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATAcc                                                              >  1:159973/4‑76 (MQ=255)               tttCAGCGTTAACGGCTNAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTATGCATACCGTCt                                                          >  1:227243/1‑76 (MQ=255)                ttCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCtt                                                         <  2:1641651/76‑1 (MQ=255)                ttCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCtt                                                         <  2:565278/76‑1 (MQ=255)                ttCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCtt                                                         <  2:558312/76‑1 (MQ=255)                 tCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTa                                                        <  2:849321/76‑1 (MQ=255)                 tCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTa                                                        <  2:1150244/76‑1 (MQ=255)                 tCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTa                                                        <  2:1150236/76‑1 (MQ=255)                  cAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACNGTCTTAt                                                       <  2:1424414/76‑1 (MQ=255)                    gCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCCGTGTGGATACCGTCTTATgg                                                     <  2:148947/76‑1 (MQ=255)                      gTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAg                                                   >  1:573408/1‑76 (MQ=255)                       ttAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGc                                                  >  2:576263/1‑76 (MQ=255)                          aCGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCaaa                                               <  1:988248/76‑1 (MQ=255)                          aCGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCaaa                                               <  2:799152/76‑1 (MQ=255)                            ggCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAAc                                             >  1:301983/1‑76 (MQ=255)                               tGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCg                                          >  1:1272666/1‑76 (MQ=255)                                gAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCGc                                         <  2:331269/76‑1 (MQ=255)                                gAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCGc                                         >  2:1520844/1‑76 (MQ=255)                                 aaGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCGCa                                        >  2:185271/1‑76 (MQ=255)                                 aaGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCGCa                                        >  2:642511/1‑76 (MQ=255)                                   ggACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGt                                      <  2:1378525/76‑1 (MQ=255)                                    gACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGTGCAAAACCCGCAGTc                                     <  2:152531/76‑1 (MQ=255)                                    gACATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTc                                     <  2:1657439/76‑1 (MQ=255)                                     aCATCGGGTCAATCGCCCAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTca                                    <  2:1524432/76‑1 (MQ=255)                                                    cccAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCTCAGTCACAGCTGAAActtctt                     <  2:158935/76‑1 (MQ=255)                                                    cccAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAActtctt                     <  1:334543/76‑1 (MQ=255)                                                     ccAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCAGAAACTTCTTg                    <  1:50900/76‑1 (MQ=255)                                                      cAACACCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGt                   >  1:366383/1‑76 (MQ=255)                                                       nacacCAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTc                  >  1:1369544/2‑76 (MQ=255)                                                            cAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCg             <  2:1600609/76‑1 (MQ=255)                                                            cAACCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCg             <  2:1146584/76‑1 (MQ=255)                                                              aCCTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTg           >  1:1360232/1‑76 (MQ=255)                                                               ccTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTgg          >  2:1483831/1‑76 (MQ=255)                                                               ccTCATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTgg          <  2:1043548/76‑1 (MQ=255)                                                                  cATGAAATCAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTaa       >  1:1394514/1‑76 (MQ=255)                                                                       aaTCAGTGTGGATATCGTCTTATGGAGNAAAACCCGAAGGCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGc  >  1:1553203/1‑76 (MQ=255)                                                                          |                                                                         AAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGTGGTAAGGAGC  >  REL606/3145773‑3145918

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