Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,531,179	A→C	V97G (GTG→GGG)	rpoH ←	RNA polymerase sigma factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,531,179	0	A	C	100.0%	225.7 / NA	53	V97G (GTG→GGG)	rpoH	RNA polymerase sigma factor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (34/19);  total (34/19)

AGTTACGCAGGACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCA  >  REL606/3531108‑3531251                                                                        |                                                                         aGTTACGCAGGACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTc                                                                      >  1:107529/1‑76 (MQ=255)   ttACGCAGGACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCg                                                                    <  1:159452/76‑1 (MQ=255)    tACGCAGGACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCgg                                                                   >  2:230571/1‑76 (MQ=255)     aCGCAGGACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCggg                                                                  >  2:73554/1‑76 (MQ=255)     aCGCAGGACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCggg                                                                  >  2:199135/1‑76 (MQ=255)      cGCAGGACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGt                                                                 <  2:997274/76‑1 (MQ=255)      cGCAGGACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGt                                                                 <  2:581145/76‑1 (MQ=255)      cGCAGGACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGt                                                                 <  1:490848/76‑1 (MQ=255)       gCAGGACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGtt                                                                <  2:147554/76‑1 (MQ=255)            aCGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAAc                                                           <  1:456989/76‑1 (MQ=255)             cGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACg                                                          >  1:545031/1‑76 (MQ=255)              gTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACgg                                                         >  1:922436/1‑76 (MQ=255)                aTTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGCTAAAACGgcg                                                       >  2:762391/1‑76 (MQ=255)                     tGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTg                                                  >  1:1609765/1‑76 (MQ=255)                      ggATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGc                                                 >  1:559331/1‑76 (MQ=255)                       gATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCt                                                >  2:623177/1‑76 (MQ=255)                        aTCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCtt                                               >  1:160722/1‑76 (MQ=255)                          ctctGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTtc                                             >  2:1511011/1‑76 (MQ=255)                           tctGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTtca                                            >  1:1204828/1‑76 (MQ=255)                              gCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTtcatca                                         <  1:216593/76‑1 (MQ=255)                              gCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTtcatca                                         >  1:811053/1‑76 (MQ=255)                               cTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAg                                        <  2:966389/76‑1 (MQ=255)                               cTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAg                                        <  2:231249/76‑1 (MQ=255)                                 ttGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGc                                      >  1:488611/1‑76 (MQ=255)                                 ttGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGc                                      >  1:1286897/1‑76 (MQ=255)                                  tGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGcc                                     >  1:1584249/1‑76 (MQ=255)                                   gATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCg                                    >  1:1049182/1‑76 (MQ=255)                                   gATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCg                                    <  2:792246/76‑1 (MQ=255)                                      ccAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATg                                 >  1:1151391/1‑76 (MQ=255)                                       cAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGt                                <  1:1607111/76‑1 (MQ=255)                                       cAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGt                                >  2:1356425/1‑76 (MQ=255)                                       cAGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGt                                <  2:16616/76‑1 (MQ=255)                                        aGTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGtt                               <  2:559610/76‑1 (MQ=255)                                         gTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCATTGCTTTCATCAGGCCGATGTTa                              <  2:1626223/76‑1 (MQ=255)                                         gTGAACGGCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTa                              <  1:992981/76‑1 (MQ=255)                                                gCGAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCt                       <  1:516051/76‑1 (MQ=255)                                                  gAAGGAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGa                     >  2:1160649/1‑76 (MQ=255)                                                     ggAGACCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATc                  >  1:142633/1‑76 (MQ=255)                                                      gagaCCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCa                 >  1:558289/1‑76 (MQ=255)                                                      gagaCCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCa                 >  1:604080/1‑76 (MQ=255)                                                       agaCCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCaa                >  1:666409/1‑76 (MQ=255)                                                       agaCCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCaa                >  2:1360970/1‑76 (MQ=255)                                                        gaCCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAACCaaa               >  2:556680/1‑76 (MQ=255)                                                         aCCAGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAAt              <  1:1329875/76‑1 (MQ=255)                                                            aGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATccn           <  1:1032991/76‑2 (MQ=255)                                                            aGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCg           >  1:93204/1‑76 (MQ=255)                                                            aGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCg           >  1:1189678/1‑76 (MQ=255)                                                            aGGCGCACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCg           >  1:430298/1‑76 (MQ=255)                                                               cgcACACCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCt        >  2:565495/1‑76 (MQ=255)                                                                 cacaCCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCtgt      >  2:1221891/1‑76 (MQ=255)                                                                   cacCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTgg    >  2:328172/1‑76 (MQ=255)                                                                   cacCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTgg    <  2:560798/76‑1 (MQ=255)                                                                   cacCCCCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTgg    >  2:5792/1‑76 (MQ=255)                                                                     ccccCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCa  <  1:1580176/76‑1 (MQ=255)                                                                     ccccCTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCa  <  1:938393/76‑1 (MQ=255)                                                                        |                                                                         AGTTACGCAGGACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTAAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCA  >  REL606/3531108‑3531251

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