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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,966,598	0	G	A	100.0%	189.2 / NA	47	G61E (GGA→GAA)	rho	transcription termination factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/24);  total (23/24)

CTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGT  >  NC_000913/3966456‑3966737                                                                                                                                               |                                                                                                                                            cTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGt                                                                                                                                      <  2:38518/150‑1 (MQ=255) cTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGt                                                                                                                                      <  2:205115/150‑1 (MQ=255) cTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGt                                                                                                                                      <  1:98534/150‑1 (MQ=255)      cACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTtcctc                                                                                                                                <  2:529281/151‑1 (MQ=255)      cACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTtcctc                                                                                                                                <  1:183608/151‑1 (MQ=255)                  aaaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGa                                                                                                                     >  1:477667/1‑150 (MQ=255)                            cTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTAcc                                                                                                          <  2:619644/151‑1 (MQ=255)                            cTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTAcc                                                                                                          <  2:289744/151‑1 (MQ=255)                            cTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTAcc                                                                                                          <  2:140733/151‑1 (MQ=255)                            cTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTAcc                                                                                                          <  1:73076/151‑1 (MQ=255)                            cTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTAcc                                                                                                          <  1:224066/151‑1 (MQ=255)                                aaaaCCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTAGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCg                                                                                                       <  2:427331/150‑1 (MQ=255)                                           cGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgatg                                                                                           <  2:545016/151‑1 (MQ=255)                                               tGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACAt                                                                                       >  1:382893/1‑151 (MQ=255)                                                    aaGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACg                                                                                  >  2:519231/1‑151 (MQ=255)                                                    aaGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACg                                                                                  <  1:217626/151‑1 (MQ=255)                                                    aaGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACg                                                                                  <  2:202864/151‑1 (MQ=255)                                                         ggACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTc                                                                              >  1:299969/1‑150 (MQ=255)                                                         ggACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGGCCTGATGACATCTACGTTTc                                                                              >  2:255885/1‑150 (MQ=255)                                                                attTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGc                                                                       >  1:609955/1‑150 (MQ=255)                                                                     tGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATc                                                                 <  1:540875/151‑1 (MQ=255)                                                                         aTCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCCCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccgc                                                              >  2:17092/1‑150 (MQ=255)                                                                          tCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccgccg                                                            >  1:440055/1‑151 (MQ=255)                                                                           ccTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGt                                                           >  1:597598/1‑151 (MQ=255)                                                                                     caCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAAcc                                                    <  1:123971/148‑1 (MQ=255)                                                                                       cGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCg                                                >  1:416840/1‑150 (MQ=255)                                                                                         cAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCCCCTCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCAc                                             <  1:148496/151‑1 (MQ=255)                                                                                             gagTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGt                                         >  2:152894/1‑151 (MQ=255)                                                                                              agTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTCCGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTgg                                          <  2:477824/149‑1 (MQ=255)                                                                                              agTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTgg                                          <  2:331305/149‑1 (MQ=255)                                                                                               gTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCCCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTg                                        <  2:39891/150‑1 (MQ=255)                                                                                                   cGAAGATATCTTTGGTGATGGCGTACTGGAGATTTTGCAGGATGAATTTGTTTTCCTCCGTTCCGCAGACAGCTCCTCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATAcc                                   <  2:329801/151‑1 (MQ=255)                                                                                                   cGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATAcc                                   <  1:396035/151‑1 (MQ=255)                                                                                                   cGAAGATATCTTGGGTGATGGCGTACGGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATAc                                    <  1:519069/150‑1 (MQ=255)                                                                                                      aGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAt                                 >  2:15006/1‑150 (MQ=255)                                                                                                       gATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAt                                 >  1:504483/1‑149 (MQ=255)                                                                                                           tCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTgg                           >  1:527401/1‑151 (MQ=255)                                                                                                           tCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTgg                           >  1:143602/1‑151 (MQ=255)                                                                                                               tGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAg                       >  2:256713/1‑151 (MQ=255)                                                                                                                ggTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGa                      >  2:353586/1‑151 (MQ=255)                                                                                                                  tGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTCCGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGAtt                    <  2:478494/151‑1 (MQ=255)                                                                                                                       gCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcc                >  1:320403/1‑150 (MQ=255)                                                                                                                       gCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcc                >  1:297742/1‑150 (MQ=255)                                                                                                                         gTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcccg              >  1:80311/1‑150 (MQ=255)                                                                                                                           aCTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCccgccg           >  2:474526/1‑151 (MQ=255)                                                                                                                                  atatTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaagaag    >  1:449064/1‑151 (MQ=255)                                                                                                                                   tatTGCAGGATGAATTTGG‑TTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGt  >  2:528098/1‑151 (MQ=255)                                                                                                                                               |                                                                                                                                            CTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGT  >  NC_000913/3966456‑3966737

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