Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,966,751	C→T	A112V (GCG→GTG)	rho →	transcription termination factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,966,751	0	C	T	100.0%	167.9 / NA	44	A112V (GCG→GTG)	rho	transcription termination factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/21);  total (23/21)

CGTTCCGCAGACAGCTCCTACCTCGCCGGT‑CCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGC  >  NC_000913/3966612‑3966896                                                                                                                                             |                                                                                                                                                  cGTCCCGCAGCCCGCTCCTCCCTC‑CCGGTCCCTGATGACCTCTACGTTCCCCCTAGCAAAACCCGCCGTTTCACCCGCCGAACGGTTAATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTgctgct                                                                                                                                             <  2:663589/146‑1 (MQ=255)    tCCGCAGACAGCTCCTACCTCGCCGGT‑CCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGtt                                                                                                                                      >  2:123018/1‑150 (MQ=255)             aGCTCCTACCTCGCCGGT‑CCTGAGGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGt                                                                                                                              >  2:588512/1‑149 (MQ=255)                     ccccGCCGGT‑CCTGATGACATCTACGGTGCCCCTACCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTc                                                                                                                       <  1:703072/145‑1 (MQ=255)                         gCCGGT‑CCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAc                                                                                                                >  2:213236/1‑151 (MQ=255)                            ggT‑CCTGATCACATCTGCGTTCCCCCTGCCCAATCCCGCCGTTTAACCCTCCGCACTGGTAATACCATCTCTGGTAAGATTCCCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCt                                                                                                              <  1:384982/150‑1 (MQ=255)                              t‑CCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGa                                                                                                            >  1:372392/1‑150 (MQ=255)                                ccTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaa                                                                                                           >  1:40001/1‑150 (MQ=255)                                    atgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGc                                                                                                      >  2:111984/1‑151 (MQ=255)                                             aCGTTTCCCCTAGCCAAACCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaaca                                                                                              <  1:287114/150‑1 (MQ=255)                                                    cccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctc                                                                                      >  1:300449/1‑151 (MQ=255)                                                    cccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCACCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctc                                                                                      >  1:50839/1‑151 (MQ=255)                                                            aaaTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgaga                                                                               >  1:523519/1‑150 (MQ=255)                                                             aacccgccgTTTCACCCTCCGCACTGGTGATCCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAc                                                                             <  1:219621/148‑1 (MQ=255)                                                             aacccgccgTTTCAACCCCCGCACTGGTGATCCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACaaa                                                                                                                 <  1:221021/112‑1 (MQ=255)                                                             aaTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACaaa                                                                                                                 >  2:220959/1‑115 (MQ=255)                                                              aTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACt                                                                            >  2:13437/1‑151 (MQ=255)                                                               tCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACtt                                                                           >  1:378719/1‑151 (MQ=255)                                                                      tttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCg                                                                    >  1:11295/1‑151 (MQ=255)                                                                          accccccGCACTGGTGATCCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCg                                                                    <  2:69101/142‑1 (MQ=255)                                                                           aCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGc                                                                >  1:417613/1‑150 (MQ=255)                                                                            ccTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCAc                                                              >  2:835435/1‑151 (MQ=255)                                                                                        tGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCg                                                   <  1:822683/150‑1 (MQ=255)                                                                                         gATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGt                                                  <  1:706726/150‑1 (MQ=255)                                                                                         gATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTc                                                 >  2:768493/1‑151 (MQ=255)                                                                                          aTACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTc                                                 <  1:485618/150‑1 (MQ=255)                                                                                                 cgcgggTAAGATTGGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGc                                              <  2:726649/142‑1 (MQ=255)                                                                                                          gATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGt                                <  2:586861/151‑1 (MQ=255)                                                                                                             tCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAAc                             >  1:494170/1‑151 (MQ=255)                                                                                                                 ccgccgAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGtt                         >  2:404249/1‑151 (MQ=255)                                                                                                                 ccgccgAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGtt                         >  2:692230/1‑151 (MQ=255)                                                                                                                  cgccgAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTc                        <  1:332851/151‑1 (MQ=255)                                                                                                                   gccgAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCt                       <  2:815134/151‑1 (MQ=255)                                                                                                                            agGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGa                  >  1:92651/1‑147 (MQ=255)                                                                                                                               tGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGAttt            >  1:598282/1‑150 (MQ=255)                                                                                                                                 aaCGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAAc         >  1:158818/1‑151 (MQ=255)                                                                                                                                  aCGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACt        <  2:417480/151‑1 (MQ=255)                                                                                                                                  aCGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACt        <  1:552917/151‑1 (MQ=255)                                                                                                                                   cGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTg       <  1:250136/151‑1 (MQ=255)                                                                                                                                    gCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGc      <  1:213298/151‑1 (MQ=255)                                                                                                                                      tATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTc    >  2:692771/1‑151 (MQ=255)                                                                                                                                       aTTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCt     <  1:112016/149‑1 (MQ=255)                                                                                                                                       aTTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTc    <  2:378595/150‑1 (MQ=255)                                                                                                                                         tttGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGCGGTAACGGTTCTACTGAAGATTTAACTGCTcgc  <  2:19876/150‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                  CGTTCCGCAGACAGCTCCTACCTCGCCGGT‑CCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGC  >  NC_000913/3966612‑3966896

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