Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,966,437	G→T	K7N (AAG→AAT)	rho →	transcription termination factor Rho

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,966,437	0	G	T	100.0%	175.5 / NA	50	K7N (AAG→AAT)	rho	transcription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (29/21);  total (29/21)

AGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATG  >  NC_000913/3966306‑3966573                                                                                                                                    |                                                                                                                                         aGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGcc                                                                                                                                  <  2:121812/140‑1 (MQ=255) aGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGcc                                                                                                                                  <  1:451632/140‑1 (MQ=255) aGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGcc                                                                                                                                  >  1:150960/1‑140 (MQ=255)   gCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCgg                                                                                                                                >  2:267408/1‑140 (MQ=255)       cGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTc                                                                                                                            >  2:707074/1‑140 (MQ=255)        gAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCt                                                                                                                           >  2:672017/1‑140 (MQ=255)        gAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCt                                                                                                                           <  2:282333/140‑1 (MQ=255)          aaaGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGa                                                                                                                         <  1:707074/140‑1 (MQ=255)          aaaGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGa                                                                                                                         <  1:376592/140‑1 (MQ=255)           aaGACATGATTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAg                                                                                                                        <  2:760861/140‑1 (MQ=255)             gACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCt                                                                                                                      >  1:513232/1‑140 (MQ=255)                      cGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCActct                                                                                                             >  1:570692/1‑140 (MQ=255)                                      gCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATgggg                                                                                             >  1:178904/1‑140 (MQ=255)                                      gCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATgggg                                                                                             >  1:451185/1‑140 (MQ=255)                                          ggATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTgg                                                                                         >  2:56983/1‑140 (MQ=255)                                                tCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAAcc                                                                                   <  1:99429/140‑1 (MQ=255)                                                tCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAAcc                                                                                   >  1:312782/1‑140 (MQ=255)                                                     ccATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCt                                                                              <  1:218864/140‑1 (MQ=255)                                                        tACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt                                                                           <  1:672017/140‑1 (MQ=255)                                                          ccATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt                                                                         >  1:36326/1‑140 (MQ=255)                                                          ccATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt                                                                         >  2:414806/1‑140 (MQ=255)                                                             ttCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCg                                                                      >  1:185328/1‑140 (MQ=255)                                                                  aacaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc                                                                 <  2:185328/140‑1 (MQ=255)                                                                    caTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCG‑AAGCAgg                                                              >  1:307291/1‑140 (MQ=255)                                                                              gagaTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAttt                                                     >  1:239723/1‑140 (MQ=255)                                                                              gagaTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGACCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTATGCAGGACATTAttt                                                     >  2:34829/1‑140 (MQ=255)                                                                               agaTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt                                                    >  1:525892/1‑140 (MQ=255)                                                                               agaTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt                                                    >  1:453749/1‑140 (MQ=255)                                                                                 aTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTg                                                  <  1:267408/140‑1 (MQ=255)                                                                                  tttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGc                                                 <  2:525892/140‑1 (MQ=255)                                                                                    tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCa                                               >  1:154062/1‑140 (MQ=255)                                                                                       cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATcc                                            >  1:615333/1‑140 (MQ=255)                                                                                         cAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTg                                          <  2:239723/140‑1 (MQ=255)                                                                                            gTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAag                                       >  2:356989/1‑140 (MQ=255)                                                                                              ttAAGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAagca                                     <  2:150960/140‑1 (MQ=255)                                                                                                 aGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTg                                                  >  2:207200/1‑124 (MQ=255)                                                                                                 aGAACTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTg                                                  <  1:207200/124‑1 (MQ=255)                                                                                                     cTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCa                              >  1:36032/1‑140 (MQ=255)                                                                                                     cTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCa                              <  2:570692/140‑1 (MQ=255)                                                                                                     cTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCa                              >  2:592729/1‑140 (MQ=255)                                                                                                     cTCACACCACTATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCa                              <  2:615333/140‑1 (MQ=255)                                                                                                               tATGAATCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCg                    <  2:513232/140‑1 (MQ=255)                                                                                                                    aTCTTACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGat               >  1:217047/1‑140 (MQ=255)                                                                                                                        tACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGa                >  1:650433/1‑135 (MQ=255)                                                                                                                        tACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGa                <  2:650433/135‑1 (MQ=255)                                                                                                                        tACCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCt           >  1:196782/1‑140 (MQ=255)                                                                                                                         aCCGAATTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCtt          >  1:331299/1‑140 (MQ=255)                                                                                                                             aaTTAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGt      >  2:202254/1‑140 (MQ=255)                                                                                                                               ttAAATAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGa    <  1:414806/140‑1 (MQ=255)                                                                                                                                 aaaTAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATg  <  2:178904/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                         AGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATG  >  NC_000913/3966306‑3966573

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