Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,967,093	A→G	E226G (GAA→GGA)	rho →	transcription termination factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,967,093	0	A	G	100.0%	134.5 / NA	38	E226G (GAA→GGA)	rho	transcription termination factor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/16);  total (22/16)

GGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTT  >  NC_000913/3966963‑3967236                                                                                                                                   |                                                                                                                                                ggTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTg                                                                                                                             <  2:46790/151‑1 (MQ=255)      aaCCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGaa                                                                                                                        <  1:382678/151‑1 (MQ=255)          aTGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAAccc                                                                                                                     >  2:513934/1‑150 (MQ=255)            gctgctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCa                                                                                                                  <  1:71151/151‑1 (MQ=255)               gctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCAtc                                                                                                                <  2:9991/150‑1 (MQ=255)                 tgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCg                                                                                                             >  1:36240/1‑151 (MQ=255)                    aGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGcc                                                                                                           >  2:536409/1‑150 (MQ=255)                           tGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCa                                                                                                    <  2:89617/150‑1 (MQ=255)                            gCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGCACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAgg                                                                                                  >  2:509230/1‑151 (MQ=255)                             cTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGt                                                                                                 >  1:5138/1‑151 (MQ=255)                             cTCAGAGCATTGCTTACAACCACCCGGATTGGGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCAAGATGCAGCGTATGGTAAAAGGTGGAGGTGTTGCGTCTACCTTTGACGAACCCGCATCTAGCCACGTTCCGGt                                                                                                 >  1:318921/1‑151 (MQ=255)                                 gagCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCg                                                                                             >  2:313997/1‑151 (MQ=255)                                     aTTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGaa                                                                                          >  1:380033/1‑150 (MQ=255)                                        gCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATgg                                                                                      <  1:654458/151‑1 (MQ=255)                                        gCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATgg                                                                                      <  2:118518/151‑1 (MQ=255)                                         cTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGt                                                                                     >  2:249381/1‑151 (MQ=255)                                            aCAACCACCCGGATTGTGTGCTGATGGTTCCGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGAt                                                                                  >  2:14733/1‑151 (MQ=255)                                                   cccGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCgaga                                                                             >  1:471363/1‑149 (MQ=255)                                                        aTTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCaa                                                                      <  1:536405/151‑1 (MQ=255)                                                            tgtgCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGc                                                                  >  2:299356/1‑151 (MQ=255)                                                               gCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCt                                                                <  1:39698/150‑1 (MQ=255)                                                               gCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCt                                                                <  2:62101/150‑1 (MQ=255)                                                               gCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCt                                                                <  1:506144/150‑1 (MQ=255)                                                                    tGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTg                                                           >  1:13874/1‑150 (MQ=255)                                                                      gTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGc                                                        >  2:479438/1‑151 (MQ=255)                                                                           gctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaa                                                    <  1:44540/150‑1 (MQ=255)                                                                                         gTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCAt                                     >  1:344841/1‑151 (MQ=255)                                                                                         gTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGGGATCGAGAAGGCCAAACGCCTCGTTGAGCACAAGACAGACGTTATCAt                                     >  1:17918/1‑151 (MQ=255)                                                                                          tCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCAtt                                    >  2:48617/1‑151 (MQ=255)                                                                                             ggAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTg                                 <  1:70392/151‑1 (MQ=255)                                                                                                     tAACCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTc                         >  2:455081/1‑151 (MQ=255)                                                                                                      aaCCGAGATGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTcc                        >  1:53869/1‑151 (MQ=255)                                                                                                           agaTGCAGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCAc                   <  2:418389/151‑1 (MQ=255)                                                                                                                 aGCGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAAAGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCt             >  1:148664/1‑151 (MQ=255)                                                                                                                   cGTCTGGTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCt             >  1:261523/1‑149 (MQ=255)                                                                                                                        ggTAAAAGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGgcgcg       <  2:344844/150‑1 (MQ=255)                                                                                                                            aaaGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCtt  <  1:457411/151‑1 (MQ=255)                                                                                                                             aaGGTGGAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCtt  >  1:449105/1‑150 (MQ=255)                                                                                                                                   |                                                                                                                                                GGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTT  >  NC_000913/3966963‑3967236

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