Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,967,246	T→C	V277A (GTT→GCT)	rho →	transcription termination factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,967,246	0	T	C	100.0%	184.4 / NA	51	V277A (GTT→GCT)	rho	transcription termination factor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/26);  total (25/26)

GCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTA  >  NC_000913/3967101‑3967392                                                                                                                                                  |                                                                                                                                                   gCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTc                                                                                                                                               <  1:875057/151‑1 (MQ=255) gCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTc                                                                                                                                               >  1:575701/1‑151 (MQ=255)   ttCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCg                                                                                                                                             >  2:171533/1‑151 (MQ=255)    tCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCgg                                                                                                                                            >  1:149494/1‑151 (MQ=255)         cTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGGTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTACGGCg                                                                                                                                          >  2:786125/1‑148 (MQ=255)            tGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTACGGCGTCAGGt                                                                                                                                    >  2:718641/1‑151 (MQ=255)              aCGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTa                                                                                                                                   >  1:888224/1‑150 (MQ=255)                 aaCCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAgt                                                                                                                               <  1:693684/151‑1 (MQ=255)                  aCCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAgtg                                                                                                                              <  2:577089/151‑1 (MQ=255)                  aCCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAgtg                                                                                                                              <  2:516224/151‑1 (MQ=255)                           tcGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCggt                                                                                                                     >  1:243457/1‑151 (MQ=255)                            cGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCggtg                                                                                                                    <  2:888153/151‑1 (MQ=255)                            cGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCggtg                                                                                                                    <  2:276106/151‑1 (MQ=255)                                     cAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATg                                                                                                           <  2:183774/151‑1 (MQ=255)                                             ggAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGcc                                                                                                   <  2:863610/151‑1 (MQ=255)                                               aaaTGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCt                                                                                                 <  1:90262/151‑1 (MQ=255)                                                aaTGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTg                                                                                                >  1:194713/1‑151 (MQ=255)                                                      gATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACGCCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCg                                                                                           <  2:787047/150‑1 (MQ=255)                                                       aTCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTc                                                                                         >  2:479193/1‑151 (MQ=255)                                                         cGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCg                                                                                       >  2:297272/1‑151 (MQ=255)                                                         cGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCg                                                                                       <  1:871334/151‑1 (MQ=255)                                                           agaAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGa                                                                                      >  1:415149/1‑150 (MQ=255)                                                            gaAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGaaa                                                                                    <  2:76574/151‑1 (MQ=255)                                                               ggCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGc                                                                                 >  1:594676/1‑151 (MQ=255)                                                                gCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGct                                                                                >  1:164670/1‑151 (MQ=255)                                                                       gCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTgg                                                                         >  2:907765/1‑151 (MQ=255)                                                                         cTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTg                                                                       >  1:846631/1‑151 (MQ=255)                                                                           ggTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTg                                                                       <  1:627263/149‑1 (MQ=255)                                                                                    caAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAAc                                                            >  1:526818/1‑151 (MQ=255)                                                                                             cGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAag                                                     >  2:426784/1‑149 (MQ=255)                                                                                                    aTTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCa                                            <  1:56289/151‑1 (MQ=255)                                                                                                       cTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGcc                                         >  1:854851/1‑151 (MQ=255)                                                                                                            cGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGAcc                                    >  2:44023/1‑151 (MQ=255)                                                                                                                  cATCACTCGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATc                              >  1:288937/1‑151 (MQ=255)                                                                                                                         cGTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACgg                       <  2:409471/151‑1 (MQ=255)                                                                                                                          gTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGgc                      <  1:656550/151‑1 (MQ=255)                                                                                                                          gTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGgc                      <  1:329849/151‑1 (MQ=255)                                                                                                                          gTCTGGCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGgc                      <  1:789621/151‑1 (MQ=255)                                                                                                                              ggCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGgcg                     <  2:594627/148‑1 (MQ=255)                                                                                                                              ggCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCt                   <  2:367537/150‑1 (MQ=255)                                                                                                                              ggCGCGCGCTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCt                   <  2:215822/150‑1 (MQ=255)                                                                                                                                cgcgcgcTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTAt                <  2:323959/151‑1 (MQ=255)                                                                                                                                 gcgcgcTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATc               >  2:856240/1‑151 (MQ=255)                                                                                                                                     gcTTACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATa           <  2:72151/151‑1 (MQ=255)                                                                                                                                        tACAACACCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATAc          <  1:556036/149‑1 (MQ=255)                                                                                                                                          caacaCCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGt      >  2:685582/1‑151 (MQ=255)                                                                                                                                          caacaCCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGt      >  1:752844/1‑151 (MQ=255)                                                                                                                                          caacaCCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGt      <  2:142985/151‑1 (MQ=255)                                                                                                                                           aacaCCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGtt     <  1:713921/151‑1 (MQ=255)                                                                                                                                           aacaCCGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGtt     >  1:748962/1‑151 (MQ=255)                                                                                                                                               ccGCTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGTGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTa  >  1:207573/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTGGTGTGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTA  >  NC_000913/3967101‑3967392

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