Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,966,727	C→T	P104L (CCG→CTG)	rho →	transcription termination factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,966,727	0	C	T	100.0%	209.5 / NA	52	P104L (CCG→CTG)	rho	transcription termination factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/26);  total (26/26)

GGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTT  >  NC_000913/3966581‑3966873                                                                                                                                                   |                                                                                                                                                   ggAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGaaa                                                                                                                                                <  2:379523/151‑1 (MQ=255)          gCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGaa                                                                                                                                       <  1:417307/151‑1 (MQ=255)             ggttggatTTGCTTTCCCCCGTCCCGCAGACAGCTCCTTCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACg                                                                                                                                     <  2:282202/147‑1 (MQ=255)             ggatggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGc                                                                                                                                    >  2:113357/1‑151 (MQ=255)              gatggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCt                                                                                                                                   >  2:924358/1‑151 (MQ=255)                  gatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGCCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTAttt                                                                                                                               <  2:274097/151‑1 (MQ=255)                  gatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTAttt                                                                                                                               <  1:435194/151‑1 (MQ=255)                       ggTTTCCTCCGTTCCGCAGACAGCTCCTCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTgcg                                                                                                                           <  2:327218/150‑1 (MQ=255)                        gTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCgct                                                                                                                         >  1:618430/1‑151 (MQ=255)                          ttCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCgctgc                                                                                                                       >  1:95819/1‑151 (MQ=255)                          ttCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCgctgc                                                                                                                       >  2:839937/1‑151 (MQ=255)                            cctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGctgctg                                                                                                                     <  1:64377/151‑1 (MQ=255)                               ccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaaa                                                                                                                  <  1:943447/151‑1 (MQ=255)                                cGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaa                                                                                                                   >  1:238194/1‑149 (MQ=255)                                    ccGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTa                                                                                                              <  2:157955/150‑1 (MQ=255)                                      gCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                                                                           >  1:305677/1‑151 (MQ=255)                                      gCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                                                                           >  1:2788/1‑151 (MQ=255)                                       cAGACAGCTCCTACCTCGCCGGCCCTGATGACATCTACGTTTCCCCTAGCCAATTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                                                                           <  2:34769/150‑1 (MQ=255)                                        aGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGaa                                                                                                         <  1:422332/151‑1 (MQ=255)                                        aGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGa                                                                                                          >  1:554219/1‑150 (MQ=255)                                        aGACAGCTCCAACCTCGCCGGTCCTGATGGCATCTACGTTCCCCCTAGCAAAAGCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGaa                                                                                                         <  2:368319/151‑1 (MQ=255)                                         gACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGa                                                                                                          >  1:774843/1‑149 (MQ=255)                                           cAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGtt                                                                                                      <  1:140319/151‑1 (MQ=255)                                                ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAAc                                                                                                   <  1:849390/149‑1 (MQ=255)                                                ccTACCTCCCCGGTCCTGATGACATTTTCGTTTCCCCTAGCCAAATCCGGCGTTTCAACCTCCGCACTTGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                                                                           <  1:619058/141‑1 (MQ=255)                                                    ccTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTc                                                                                                <  1:787376/148‑1 (MQ=255)                                                    ccTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                                                                             <  1:198057/151‑1 (MQ=255)                                                    ccTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                                                                             <  2:155586/151‑1 (MQ=255)                                                           ggTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                                                                                      >  2:214240/1‑151 (MQ=255)                                                              ccTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaa                                                                                    >  2:605254/1‑150 (MQ=255)                                                                tgatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAAc                                                                                 >  2:801226/1‑151 (MQ=255)                                                                            cGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGGTGAAAGTTAACGAATTTAACTTCGACAAGCCTGAAAACGCCCGCAACaaa                                                                     <  2:520605/151‑1 (MQ=255)                                                                             gTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaa                                                                     >  1:852742/1‑150 (MQ=255)                                                                                    cTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCtt                                                             >  2:881567/1‑151 (MQ=255)                                                                                      aGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt                                                            <  1:522042/150‑1 (MQ=255)                                                                                      aGCCAAATCCGCCGCTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt                                                            <  1:197983/150‑1 (MQ=255)                                                                                      aGCCAAATCCGCCGATTCAACCTCCGCACTGGTGATACCACCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg                                                           >  2:20324/1‑151 (MQ=255)                                                                                       gCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg                                                           >  1:374690/1‑150 (MQ=255)                                                                                       gCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACACCCGCAACAAAATCCTCTTTg                                                           >  1:844284/1‑150 (MQ=255)                                                                                        ccAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt                                                            >  2:512348/1‑148 (MQ=255)                                                                                            aTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACt                                                     <  1:91353/151‑1 (MQ=255)                                                                                               cgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTa                                                   >  2:739834/1‑150 (MQ=255)                                                                                                    tttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCg                                             >  2:947637/1‑151 (MQ=255)                                                                                                     ttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCg                                             <  2:265924/150‑1 (MQ=255)                                                                                                          ccTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCAc                                       >  1:479291/1‑151 (MQ=255)                                                                                                            tCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGc                                     >  1:284432/1‑151 (MQ=255)                                                                                                               gCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCaa                                   >  2:130062/1‑150 (MQ=255)                                                                                                                    ggTGATACCATCTCTGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACtctc                             <  1:706455/151‑1 (MQ=255)                                                                                                                                tctGGTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGACCGCAACAAAATCCTCTTTGAGAACTTAAACCCGCTGCACGCAAACTCTCGTCTGCGTATgg                 >  2:329562/1‑151 (MQ=255)                                                                                                                                   ggTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATgg                 <  1:554671/148‑1 (MQ=255)                                                                                                                                   ggTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAAc              >  1:502354/1‑151 (MQ=255)                                                                                                                                    gTAAGATTCGCCCGCTGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACg             <  1:365037/151‑1 (MQ=255)                                                                                                                                            cgcccgctGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACg     <  1:737737/151‑1 (MQ=255)                                                                                                                                             gcccgctGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACgg    >  1:419178/1‑151 (MQ=255)                                                                                                                                               ccgctgAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGtt  >  1:712626/1‑151 (MQ=255)                                                                                                                                               ccgctgAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGtt  >  1:387458/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTT  >  NC_000913/3966581‑3966873

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