Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,316,792	G→A	I245I (ATC→ATT)	nusA ←	transcription termination/antitermination protein NusA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,316,792	0	G	A	100.0%	229.7 / NA	64	I245I (ATC→ATT)	nusA	transcription termination/antitermination protein NusA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/29);  total (35/29)

TTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGG  >  NC_000913/3316665‑3316923                                                                                                                                |                                                                                                                                    ttGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCg                                                                                                                         <  1:478667/140‑1 (MQ=255) ttGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCg                                                                                                                         <  1:269461/140‑1 (MQ=255)  tGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGt                                                                                                                        >  2:279006/1‑140 (MQ=255)                aCTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACCCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACcgcg                                                                                                          >  1:128917/1‑140 (MQ=255)                    cgcCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGAttt                                                                                                      >  1:590643/1‑140 (MQ=255)                    cgcCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGAttt                                                                                                      >  1:76160/1‑140 (MQ=255)                    cgcCGGGTTATCATCCCACAGGACGATATCGATACACTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGAttt                                                                                                      >  2:719185/1‑140 (MQ=255)                       cGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCg                                                                                                   <  2:392933/140‑1 (MQ=255)                          gTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCAc                                                                                                >  1:192698/1‑140 (MQ=255)                           ttATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACg                                                                                               <  2:430022/140‑1 (MQ=255)                           ttATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACg                                                                                               <  2:192698/140‑1 (MQ=255)                           ttATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACACAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACg                                                                                               <  1:231092/140‑1 (MQ=255)                               catcCCACAGGACGATATCGATACACTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGaa                                                                                           <  1:719185/140‑1 (MQ=255)                                        ggACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATcg                                                                                  >  2:250819/1‑140 (MQ=255)                                        ggACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATcg                                                                                  >  1:897981/1‑140 (MQ=255)                                               atCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGGGCGCCACGCATACCTTCGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCg                                                                           <  2:626499/140‑1 (MQ=255)                                                        cTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTaa                                                                  >  1:273807/1‑140 (MQ=255)                                                            ccACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGGTCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTa                                                                   >  1:903773/1‑135 (MQ=255)                                                            ccACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGGTCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTa                                                                   <  2:903773/135‑1 (MQ=255)                                                                ccAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATc                                                          >  2:982053/1‑140 (MQ=255)                                                                     tCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCActtc                                                     >  2:290232/1‑140 (MQ=255)                                                                       aGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCActtctt                                                   >  2:1003906/1‑140 (MQ=255)                                                                       aGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCActtctt                                                   >  1:527549/1‑140 (MQ=255)                                                                            acacCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCg                                              <  2:741031/140‑1 (MQ=255)                                                                            acacCGCCTGAACACGCGCGCCACGCATACCTACCCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCg                                              <  2:128917/140‑1 (MQ=255)                                                                             cacCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGa                                             >  1:682369/1‑140 (MQ=255)                                                                              acCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGAt                                            >  1:582813/1‑140 (MQ=255)                                                                              acCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGAt                                            >  2:13205/1‑140 (MQ=255)                                                                              acCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGAt                                            >  1:1063750/1‑140 (MQ=255)                                                                                cGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGAttt                                          >  1:703432/1‑140 (MQ=255)                                                                                cGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGAttt                                          >  2:376295/1‑140 (MQ=255)                                                                                cGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGAttt                                          >  2:373361/1‑140 (MQ=255)                                                                                cGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGAttt                                          >  1:12910/1‑140 (MQ=255)                                                                                 gCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTc                                         >  2:1060601/1‑140 (MQ=255)                                                                                  ccTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCt                                        <  2:273807/140‑1 (MQ=255)                                                                                  ccTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCt                                        >  1:838309/1‑140 (MQ=255)                                                                                   cTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTg                                       <  1:13205/140‑1 (MQ=255)                                                                                   cTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTg                                       >  2:582592/1‑140 (MQ=255)                                                                                    tGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTgg                                      <  1:373361/140‑1 (MQ=255)                                                                                     gAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGc                                     >  2:852111/1‑140 (MQ=255)                                                                                      aaCACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCa                                    >  1:533455/1‑140 (MQ=255)                                                                                       acacGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAc                                   <  1:852111/140‑1 (MQ=255)                                                                                         acGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACtt                                 <  2:533455/140‑1 (MQ=255)                                                                                         acGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACtt                                 <  1:279006/140‑1 (MQ=255)                                                                                         acGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACtt                                 <  2:897981/140‑1 (MQ=255)                                                                                         acGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACtt                                 <  1:290232/140‑1 (MQ=255)                                                                                         acGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACtt                                 <  1:376295/140‑1 (MQ=255)                                                                                          cgcgcgCCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTc                                >  2:973688/1‑140 (MQ=255)                                                                                           gcgcgcCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCa                               >  1:256239/1‑140 (MQ=255)                                                                                             gcgcCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAAt                             <  2:838309/140‑1 (MQ=255)                                                                                              cgcCACGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATa                            <  2:1063750/140‑1 (MQ=255)                                                                                                  aCGCATACCTACGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGa                        <  1:582592/140‑1 (MQ=255)                                                                                                             cGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATc             <  1:250819/140‑1 (MQ=255)                                                                                                             cGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATc             >  2:274801/1‑140 (MQ=255)                                                                                                             cGCAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACCTCAATACGGAACAGTTCGATc             >  2:915541/1‑140 (MQ=255)                                                                                                               cAAGCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAg           >  2:508644/1‑140 (MQ=255)                                                                                                                  gCACCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCAt        >  2:345033/1‑140 (MQ=255)                                                                                                                    aCCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGCGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCAttt      <  2:666274/140‑1 (MQ=255)                                                                                                                    aCCTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCAttt      <  1:982053/140‑1 (MQ=255)                                                                                                                     ccTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTc     <  1:973688/140‑1 (MQ=255)                                                                                                                      cTACCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTcc    <  2:582813/140‑1 (MQ=255)                                                                                                                        aCCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCAttt      >  2:373494/1‑136 (MQ=255)                                                                                                                        aCCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCAttt      <  1:373494/136‑1 (MQ=255)                                                                                                                        aCCGGATCAATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCgg  >  1:110086/1‑140 (MQ=255)                                                                                                                                |                                                                                                                                    TTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGG  >  NC_000913/3316665‑3316923

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