Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	5,948,434	G→A	L14F (CTT→TTT)	rho ←	transcription termination factor Rho

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,948,434	0	G	A	100.0%	170.1 / NA	50	L14F (CTT→TTT)	rho	transcription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/22);  total (28/22)

AACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGA  >  NC_002947/5948290‑5948577                                                                                                                                                 |                                                                                                                                                aaccATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATc                                                                                                                                              >  6:21075/1‑148 (MQ=255)  accATCCTGGAGAAGCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCg                                                                                                                                             <  4:122184/148‑1 (MQ=255)      tCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTaa                                                                                                                                         >  1:34400/1‑148 (MQ=255)      tCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTaa                                                                                                                                         >  8:101864/1‑148 (MQ=255)         tGGATAATCTCCAGCACTCCGTAACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATc                                                                                                                                       <  8:151440/147‑1 (MQ=255)         tGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCg                                                                                                                                      >  2:34755/1‑148 (MQ=255)           gagaATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGc                                                                                                                                    >  7:81706/1‑148 (MQ=255)           gagaATCTCCAGCACGCCGTCAACCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGTGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATTCCCATCTGTTCTGCCATTTCCAAAAAATCGGTAATCGGc                                                                                                                                    >  8:80787/1‑148 (MQ=255)                gctcCATCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCTTGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTTCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTg                                                                                                                               <  7:122620/147‑1 (MQ=255)                  tcCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCt                                                                                                                             >  2:121556/1‑148 (MQ=255)                    cAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCt                                                                                                                             >  1:139865/1‑146 (MQ=255)                    cAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCt                                                                                                                             >  8:157903/1‑146 (MQ=255)                     aGCTCGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTGCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGa                                                                                                                          <  4:121164/148‑1 (MQ=255)                     aGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGa                                                                                                                          <  4:155483/148‑1 (MQ=255)                     aGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGa                                                                                                                          <  7:118688/148‑1 (MQ=255)                     aGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGa                                                                                                                          <  6:6679/148‑1 (MQ=255)                        aCGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGt                                                                                                                        >  2:81551/1‑147 (MQ=255)                                  ccGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCgg                                                                                                                                     >  2:160088/1‑124 (MQ=255)                                   cGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGc                                                                                                                                    <  1:160088/124‑1 (MQ=255)                                          tctTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGCGCCATGTTTTCGATGCCCCTCTGTTCGGCCATTTCAAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATaatg                                                                                                      >  2:23875/1‑145 (MQ=255)                                                     tttCGCATGCTTCTTCAGCAGGGCGAAAATCAAGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGCGGAATGACGTa                                                                                           >  8:128029/1‑147 (MQ=255)                                                        cGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATc                                                                                        <  3:99037/147‑1 (MQ=255)                                                          cATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCAt                                                                                      <  3:104443/147‑1 (MQ=255)                                                          cATGCTTCTTCAGAAGGGCGAAAATCACGTCCTGTTTCCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTATTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCAt                                                                                      >  3:134561/1‑147 (MQ=255)                                                           aTGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCgg                                                                                                                                     >  7:84135/1‑99 (MQ=255)                                                           aTGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATg                                                                                     >  3:51471/1‑147 (MQ=255)                                                            tGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGc                                                                                                                                    <  8:84135/99‑1 (MQ=255)                                                               ttcttcAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTaa                                                                                  >  6:168493/1‑146 (MQ=255)                                                               ttcttcAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTAAGTCAGGTTCATAAGGGGAATTCCGTAATCATGTaa                                                                                  >  2:21/1‑146 (MQ=255)                                                                         gggCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAAtt                                                                      <  4:79058/148‑1 (MQ=255)                                                                          ggCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATg                                                                                     >  8:115454/1‑132 (MQ=255)                                                                           gCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGt                                                                                    <  7:115454/132‑1 (MQ=255)                                                                           gCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTaa                                                                    <  3:80558/148‑1 (MQ=255)                                                                                  tCACGTACTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGACAATTAAGCTTCTg                                                             <  6:90036/148‑1 (MQ=255)                                                                                   cACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTgg                                                            <  7:69913/148‑1 (MQ=255)                                                                                        ccTGTTTGCGCGAACGGGCCATGTTTTCTATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTaa                                                       <  1:154354/148‑1 (MQ=255)                                                                                        ccTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTaa                                                       <  2:23545/148‑1 (MQ=255)                                                                                          tGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTaa                                                       >  6:20577/1‑146 (MQ=255)                                                                                                 gcgAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGgccgc                                              >  1:135445/1‑148 (MQ=255)                                                                                                   gAACTGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGgccgcgc                                            >  5:86833/1‑148 (MQ=255)                                                                                                   gAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGCAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGgccgcgc                                            >  4:57136/1‑148 (MQ=255)                                                                                                       gggCAATCTTTTCGATTACCATCTCTACGGCCATTTCCAAAAAATCGGAAATCGGATTTTGCTTGAGTTCAGTCAGGTTCATAAGAGGAATGACGTAATCATGTAAGAATGGAGAATTAACCTTATGCCTTAATGAGGCCGCGCCGCt                                        >  8:22253/1‑148 (MQ=255)                                                                                                          ccATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGa                                     >  4:132518/1‑148 (MQ=255)                                                                                                             tGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATgg                                  <  4:27097/148‑1 (MQ=255)                                                                                                                   cGATGCCCATCTGTTCGGCCATTTCCAAAAAATCTGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACa                             <  8:143153/147‑1 (MQ=255)                                                                                                                    gATGCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCAGCGCCGCTGGATGGCGACAgg                           >  2:48235/1‑148 (MQ=255)                                                                                                                       gCCCATCTGTTCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATc                        <  7:19671/148‑1 (MQ=255)                                                                                                                                ttCGGCCATTTCCAAAAAATCGGTAATTGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGa               >  4:88135/1‑148 (MQ=255)                                                                                                                                 tCGGCCATTTCCAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGa               <  5:21075/147‑1 (MQ=255)                                                                                                                                            cAAAAAATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAggg   >  5:175351/1‑148 (MQ=255)                                                                                                                                             aaaaaaTCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGa  >  6:4447/1‑148 (MQ=255)                                                                                                                                             aaaaaaTCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGa  >  8:79574/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                                AACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCATGTAAGAAGGGAGAATTAAGCTTCTGGCTTAATGAGGCCGCGCCGCTGGATGGCGACAGGATCGCGTACTGATTCGAATTAGGGA  >  NC_002947/5948290‑5948577

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