Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,966,600	T→C	F62L (TTT→CTT)	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,600	0	T	C	83.0%	109.3 / 13.2	47	F62L (TTT→CTT)	rho	transcription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (4/4);  new base C (19/20);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

GATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAA  >  NC_000913/3966385‑3966840                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  taTTAACCCCAAGTTTACGAACTCACACCACTATGAATCTTACCGAATTAACGACTACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCATGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgat                                                                                                                                                                                                    <  2:89810/262‑1 (MQ=255)                                tATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATAATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGAGGGACTGGGTTTCCTCCGTTCCGCAGACATGCCCTACCCTGCCGGTCACGATGACATCTACGTGTTCCATGTGCAAAACCCCCCCTTTCAACcccc                                                                                                                                                        >  2:445617/1‑276 (MQ=255)                                   gAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCTCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATccg                                                                                                                                                                      <  1:316748/258‑1 (MQ=255)                                       cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGc                                                                                                                                                                                                                                >  2:302208/1‑197 (MQ=255)                                          aCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGTCTTGGTTtcc                                                                                                                                                                                                                                          >  1:14106/1‑184 (MQ=255)                                          aCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGa                                                                                                                                                                                                                             <  1:302211/197‑1 (MQ=255)                                             gAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTtcctcc                                                                                                                                                                                                                                       <  2:14106/184‑1 (MQ=255)                                                                          gATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGc                                                                                                                                                                                                                                >  2:112365/1‑162 (MQ=255)                                                                             cACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGa                                                                                                                                                                                                                             <  1:112365/162‑1 (MQ=255)                                                                                             atGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGGTATCGTGGGTGTTGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAGTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgcg                                                                                          <  1:445621/276‑1 (MQ=255)                                                                                                                cTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTg                                                                                                        >  1:224212/1‑243 (MQ=255)                                                                                                                cTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCt                                                                                                                                                                                                       >  2:456192/1‑149 (MQ=255)                                                                                                                   gTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGACATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgat                                                                                                                                                                                                    <  1:456196/149‑1 (MQ=255)                                                                                                                             gCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTTGTGATGGCGTACTTGTGATATTGCAGTATGGACTTGGTTTCCTCCGTTACTCAGACAGCTCCTCCCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATccg                                                                                                                                                                      <  1:331072/168‑1 (MQ=255)                                                                                                                                   cATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAAt                                                                                                                                                                         >  2:85798/1‑159 (MQ=255)                                                                                                                                      tattTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGt                                                                                                                                                                  >  2:372943/1‑163 (MQ=255)                                                                                                                                      tattTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtg                                                                                                                                                                                                      >  1:403545/1‑128 (MQ=255)                                                                                                                                      tattTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATccg                                                                                                                                                                      <  1:85798/159‑1 (MQ=255)                                                                                                                                         ttttGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGTTTc                                                                                                                                                               <  1:372947/163‑1 (MQ=255)                                                                                                                                         ttttGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgatg                                                                                                                                                                                                   <  2:403541/128‑1 (MQ=255)                                                                                                                                                            caCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTcc                                                                                                                                                                                                                      >  2:414303/1‑90 (MQ=255)                                                                                                                                                            caCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCt                                                                                                                                                                                                       >  2:181005/1‑105 (MQ=255)                                                                                                                                                             aCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACGTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTCGCCAATTCCGCCGTTTCAACCTCCGCACTGGTGTTACCATCTCTGGTAAGATTcgcccgcc                                                                                                                   >  1:439632/1‑187 (MQ=255)                                                                                                                                                               gCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTAc                                                                                                                                                                                                                   <  1:414307/90‑1 (MQ=255)                                                                                                                                                               gCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgat                                                                                                                                                                                                    <  1:181006/105‑1 (MQ=255)                                                                                                                                                                cAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgatga                                                                                                                                                                                                  >  1:75637/1‑106 (MQ=255)                                                                                                                                                                   agagTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACAt                                                                                                                                                                                               <  2:75637/106‑1 (MQ=255)                                                                                                                                                                         gCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAAt                                                                                                                                                                         >  2:260640/1‑121 (MQ=255)                                                                                                                                                                            aaGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATccg                                                                                                                                                                      <  1:260642/121‑1 (MQ=255)                                                                                                                                                                               atatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATcc                                                                                                                                                                       >  1:315461/1‑117 (MQ=255)                                                                                                                                                                               atatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGTTTCAACCTcc                                                                                                                                                        >  2:123085/1‑132 (MQ=255)                                                                                                                                                                               atatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTACGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGTTTCAACCTcc                                                                                                                                                        >  1:279641/1‑132 (MQ=255)                                                                                                                                                                                  tCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATccgcc                                                                                                                                                                    <  2:315458/117‑1 (MQ=255)                                                                                                                                                                                  tCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCa                                                                                                                                                     <  1:123085/132‑1 (MQ=255)                                                                                                                                                                                  tCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCa                                                                                                                                                     <  2:279638/132‑1 (MQ=255)                                                                                                                                                                                  tCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTGTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATccg                                                                                                                                                                      <  2:286250/115‑1 (MQ=255)                                                                                                                                                                                         tGATGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCg                                                              >  2:286602/1‑212 (MQ=255)                                                                                                                                                                                           aTGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCa                                                                                                                                         >  2:198050/1‑135 (MQ=255)                                                                                                                                                                                            tGGCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCGACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACa                                                           <  1:286605/212‑1 (MQ=255)                                                                                                                                                                                             ggCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTACCGAAGTTAACTTCGACAAACCTGAAACCGCCCGCAACAAAATCCCCTTTTAGCACTTAAacccccgtcccgctt  >  2:21318/1‑257 (MQ=255)                                                                                                                                                                                              gCGTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAtct                                                                                                                                      <  1:198051/135‑1 (MQ=255)                                                                                                                                                                                                gTACTGGAGATATTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTga                         <  1:21318/242‑1 (MQ=255)                                                                                                                                                                                                           atTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCt                                                                                                                                                                                                       >  1:449391/1‑58 (MQ=255)                                                                                                                                                                                                           atTGCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAATAAAATCCTCtt                            >  2:253573/1‑228 (MQ=255)                                                                                                                                                                                                              gCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgat                                                                                                                                                                                                    <  2:449387/58‑1 (MQ=255)                                                                                                                                                                                                              gCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGt                                                                                                                                                                  >  1:308728/1‑91 (MQ=255)                                                                                                                                                                                                              gCAGGATGGACTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAATAAAATCCTCTTTga                         <  1:253575/228‑1 (MQ=255)                                                                                                                                                                                                                 ggatggacTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGTTTc                                                                                                                                                               <  2:308725/91‑1 (MQ=255)                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  GATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGT‑TTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAA  >  NC_000913/3966385‑3966840

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑