Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,966,727	C→T	P104L (CCG→CTG)	rho →	transcription termination factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,966,727	0	C	T	100.0%	172.5 / NA	45	P104L (CCG→CTG)	rho	transcription termination factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/26);  total (19/26)

TGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACG  >  NC_000913/3966580‑3966870                                                                                                                                                    |                                                                                                                                                tGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGa                                                                                                                                               >  2:532775/1‑150 (MQ=255)     gaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAaaga                                                                                                                                           >  2:437304/1‑150 (MQ=255)        atTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGt                                                                                                                                       >  2:24678/1‑151 (MQ=255)         ttGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGt                                                                                                                                       <  2:383509/150‑1 (MQ=255)          tGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTg                                                                                                                                      <  1:408871/150‑1 (MQ=255)          tGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTg                                                                                                                                      >  2:86278/1‑150 (MQ=255)           gCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGaa                                                                                                                                    <  1:359357/151‑1 (MQ=255)                   gatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTAttt                                                                                                                            >  1:51258/1‑151 (MQ=255)                   gatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCAGATGACATCGACGTGTCCCCTAGCCAAATCCGCCGTTTCAACCACCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTAttt                                                                                                                            <  1:112404/151‑1 (MQ=255)                    atTTGGTTTCCTCCGTTCCGCAGCAAGCTCCTCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTAttt                                                                                                                            <  1:383014/150‑1 (MQ=255)                      ttGGTTTCCTCCGTTCCGCAGCCAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTg                                                                                                                          <  2:50494/150‑1 (MQ=255)                      ttGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTgc                                                                                                                         >  2:70761/1‑151 (MQ=255)                      ttGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTg                                                                                                                          <  2:276747/150‑1 (MQ=255)                      ttGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAAACCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTgc                                                                                                                         >  2:271503/1‑151 (MQ=255)                        ggTTTCCCCCGTTCCGCAGACAGCTCCTCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTgcgc                                                                                                                       <  1:227697/151‑1 (MQ=255)                          tttCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCgctg                                                                                                                     >  1:515386/1‑151 (MQ=255)                             cctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCgctgct                                                                                                                   <  1:322977/150‑1 (MQ=255)                                ccGTTCCGCAGAAAGCTCCTACCTCGCCGGTCCTGATGACATCTGCGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaaa                                                                                                               <  1:407227/151‑1 (MQ=255)                                  gTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGt                                                                                                             >  2:370986/1‑151 (MQ=255)                                     ccGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTa                                                                                                           >  1:191497/1‑150 (MQ=255)                                        cAGACAGCTCCTCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                                                                        <  1:466101/150‑1 (MQ=255)                                              gCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTaa                                                                                                 <  2:364090/151‑1 (MQ=255)                                               cTCCTACCTCGCCGGTCCTGATGACAGCTCCGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTa                                                                                                  <  1:401213/149‑1 (MQ=255)                                                 ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACt                                                                                               <  2:504205/150‑1 (MQ=255)                                                 ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACt                                                                                               <  1:208827/150‑1 (MQ=255)                                                  cTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACt                                                                                               <  1:74801/149‑1 (MQ=255)                                                        cGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaa                                                                                        <  2:416232/150‑1 (MQ=255)                                                        cGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACa                                                                                         <  2:93044/149‑1 (MQ=255)                                                               ccTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaa                                                                                 <  2:365521/150‑1 (MQ=255)                                                                cTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaa                                                                                >  1:326341/1‑150 (MQ=255)                                                                       cATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                                                                        >  1:198063/1‑151 (MQ=255)                                                                            aCGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaac                                                                     <  1:356125/149‑1 (MQ=255)                                                                                         ccAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgag                                                      <  2:31284/151‑1 (MQ=255)                                                                                                     tttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCg                                          >  1:355144/1‑151 (MQ=255)                                                                                                          aCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCa                                     <  1:86252/151‑1 (MQ=255)                                                                                                                   cTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGCCAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAActc                            >  1:286473/1‑151 (MQ=255)                                                                                                                   cTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAActc                            >  2:323914/1‑151 (MQ=255)                                                                                                                    tGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGTTATTTTGCGCTGCTGAAAGTTAACGCAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAActc                            >  1:238647/1‑150 (MQ=255)                                                                                                                        gATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTc                       <  2:253042/151‑1 (MQ=255)                                                                                                                                ctctGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATg               >  2:240765/1‑151 (MQ=255)                                                                                                                                  ctGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTAt                <  1:14493/148‑1 (MQ=255)                                                                                                                                   tGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGa             >  1:527779/1‑150 (MQ=255)                                                                                                                                       aaGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGt         <  1:126089/150‑1 (MQ=255)                                                                                                                                        aGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTgg       >  1:504125/1‑151 (MQ=255)                                                                                                                                            tCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTaa    <  1:370890/150‑1 (MQ=255)                                                                                                                                             cgcccgctGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACg  <  2:391844/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                TGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACG  >  NC_000913/3966580‑3966870

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: ‑