Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,967,132	A→G	H239R (CAC→CGC)	rho →	transcription termination factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,967,132	0	A	G	100.0%	143.2 / NA	41	H239R (CAC→CGC)	rho	transcription termination factor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/24);  total (17/24)

TGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTG  >  NC_000913/3966989‑3967271                                                                                                                                                |                                                                                                                                            tGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTc                                                                                                                                        >  2:145011/1‑149 (MQ=255)   cTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAg                                                                                                                                      <  2:96743/149‑1 (MQ=255)          cATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCgg                                                                                                                               >  1:514262/1‑149 (MQ=255)          cATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGa                                                                                                                              >  2:130435/1‑150 (MQ=255)                 tACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTg                                                                                                                       <  1:672530/150‑1 (MQ=255)                 tACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTg                                                                                                                       <  1:385653/150‑1 (MQ=255)                                tgtgtgCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCaaa                                                                                                        <  2:252608/150‑1 (MQ=255)                                 gtgtgCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACg                                                                                                      <  2:206903/151‑1 (MQ=255)                                     gCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCt                                                                                                   <  2:639896/150‑1 (MQ=255)                                           ggTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGa                                                                                             <  1:59966/150‑1 (MQ=255)                                              tCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGca                                                                                          <  1:99012/150‑1 (MQ=255)                                               ctgctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGcac                                                                                         >  1:456869/1‑150 (MQ=255)                                                    gATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaa                                                                                    >  1:401992/1‑150 (MQ=255)                                                       cgacgaACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGAc                                                                                <  1:37231/151‑1 (MQ=255)                                                        gacgaACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACg                                                                               <  2:627780/151‑1 (MQ=255)                                                           gaACGTCAGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGtt                                                                             >  1:69555/1‑150 (MQ=255)                                                             aCGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGCGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGACATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACACGAAAGACGTTATc                                                                          >  1:633677/1‑151 (MQ=255)                                                                   ggAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTg                                                                    <  2:448773/151‑1 (MQ=255)                                                                   ggAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTg                                                                    <  2:295922/151‑1 (MQ=255)                                                                    gaagaaGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGc                                                                   <  2:249038/151‑1 (MQ=255)                                                                      agaagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTc                                                                 <  1:557429/151‑1 (MQ=255)                                                                                gagaTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCa                                                       <  1:12421/151‑1 (MQ=255)                                                                                 agaTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCAc                                                      >  1:417213/1‑151 (MQ=255)                                                                                   aTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTc                                                    <  1:334599/151‑1 (MQ=255)                                                                                          gTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGgc                                             >  1:474253/1‑151 (MQ=255)                                                                                            cTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGgcg                                            >  1:344475/1‑150 (MQ=255)                                                                                                     ggTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCTCGCGCTTaca                                  <  1:425523/151‑1 (MQ=255)                                                                                                     ggTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTac                                   >  1:661255/1‑150 (MQ=255)                                                                                                        gAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGATCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCTCTCGTCTGGCGCGCGCTTacaaca                               <  1:424477/151‑1 (MQ=255)                                                                                                        gAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTacaac                                <  1:611840/150‑1 (MQ=255)                                                                                                             tgttgCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCgtt                          >  2:605297/1‑151 (MQ=255)                                                                                                              gttgCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCgttg                         <  1:287164/151‑1 (MQ=255)                                                                                                                  cTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTcc                     >  2:105782/1‑151 (MQ=255)                                                                                                                    tCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTcc                     <  1:490559/149‑1 (MQ=255)                                                                                                                    tCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCgg                   >  1:501895/1‑151 (MQ=255)                                                                                                                    tCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCg                    <  1:195941/150‑1 (MQ=255)                                                                                                                        ccTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGCGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACTTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCTCGCTTACAACACCCTTGTTCCGGCg                 <  1:299807/149‑1 (MQ=255)                                                                                                                          tttGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTc               >  1:448889/1‑149 (MQ=255)                                                                                                                           ttGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCa              <  1:188558/149‑1 (MQ=255)                                                                                                                                      gCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTg  >  2:418221/1‑150 (MQ=255)                                                                                                                                      gCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTg  >  2:527949/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                            TGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTG  >  NC_000913/3966989‑3967271

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