Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,967,714	C→T	intergenic (+38/‑202)	rho → / → wecA	transcription termination factor/UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,967,714	0	C	T	100.0%	206.5 / NA	56	intergenic (+38/‑202)	rho/wecA	transcription termination factor/UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/37);  total (19/37)

CGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTG  >  NC_000913/3967566‑3967859                                                                                                                                                     |                                                                                                                                                  cgcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGgtgt                                                                                                                                                 <  1:716239/151‑1 (MQ=255)  gcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGtt                                                                                                                                                <  1:607400/151‑1 (MQ=255)  gcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGtt                                                                                                                                                <  1:464251/151‑1 (MQ=255)       aTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCt                                                                                                                                           <  2:321722/151‑1 (MQ=255)            tCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTta                                                                                                                                       <  2:441588/150‑1 (MQ=255)            tCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTta                                                                                                                                       <  1:20073/150‑1 (MQ=255)             cACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTtata                                                                                                                                     <  2:676921/151‑1 (MQ=255)              aCCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTatat                                                                                                                                    <  1:31068/151‑1 (MQ=255)                  gATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGt                                                                                                                                >  2:454056/1‑151 (MQ=255)                   aTGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGt                                                                                                                                >  1:176053/1‑150 (MQ=255)                   aTGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTa                                                                                                                               <  1:498161/151‑1 (MQ=255)                        cGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCtt                                                                                                                          <  2:779542/151‑1 (MQ=255)                          aaaTCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTa                                                                                                                         >  1:223699/1‑150 (MQ=255)                            aTCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATg                                                                                                                      <  2:12883/151‑1 (MQ=255)                            aTCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATg                                                                                                                      >  2:281313/1‑151 (MQ=255)                                       ggAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGc                                                                                                            <  1:593054/150‑1 (MQ=255)                                       ggAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCg                                                                                                           <  1:318996/151‑1 (MQ=255)                                            tCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGtt                                                                                                      >  2:741759/1‑151 (MQ=255)                                             ccTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTa                                                                                                     <  2:724044/151‑1 (MQ=255)                                              cTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGtt                                                                                                      >  1:594431/1‑149 (MQ=255)                                                 aTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGaa                                                                                                 <  1:824895/151‑1 (MQ=255)                                                    aataaACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAAt                                                                                              <  2:636366/151‑1 (MQ=255)                                                       aaaCTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTcc                                                                                           >  1:261139/1‑151 (MQ=255)                                                        aaCTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCt                                                                                          <  1:855564/151‑1 (MQ=255)                                                         aCTTGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTg                                                                                         <  1:55901/151‑1 (MQ=255)                                                            ggCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGa                                                                                       <  1:125845/150‑1 (MQ=255)                                                               aaTGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAAttt                                                                                   <  2:757084/151‑1 (MQ=255)                                                                aTGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTg                                                                                  <  1:179995/151‑1 (MQ=255)                                                                  gaccaagaccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCt                                                                                >  2:162895/1‑151 (MQ=255)                                                                  gaccaagaccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCt                                                                                <  1:59254/151‑1 (MQ=255)                                                                    ccaagaccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTgg                                                                              >  2:796008/1‑151 (MQ=255)                                                                            aaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTAt                                                                      <  2:265628/151‑1 (MQ=255)                                                                                  gATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAGTCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAAtt                                                                <  1:498300/151‑1 (MQ=255)                                                                                      tcttcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGc                                                             >  1:12745/1‑150 (MQ=255)                                                                                      tcttcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGc                                                             >  2:676709/1‑150 (MQ=255)                                                                                         tcGAAATGATGAAACGCTCATAGATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCa                                                            <  2:561145/148‑1 (MQ=255)                                                                                         tcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCatat                                                         >  2:560608/1‑151 (MQ=255)                                                                                         tcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCata                                                          <  1:553291/150‑1 (MQ=255)                                                                                         tcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCa                                                            <  1:551899/148‑1 (MQ=255)                                                                                           gAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATc                                                        >  1:615008/1‑150 (MQ=255)                                                                                             aaTGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCaaa                                                     <  1:577323/151‑1 (MQ=255)                                                                                                        gCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAAtttt                                          <  2:434340/151‑1 (MQ=255)                                                                                                            aTAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGc                                       >  2:205550/1‑150 (MQ=255)                                                                                                              aaaTTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGcaca                                    <  1:361062/151‑1 (MQ=255)                                                                                                               aaTTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCaa                                                                  >  1:50357/1‑120 (MQ=255)                                                                                                               aaTTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCaa                                                                  <  2:50379/120‑1 (MQ=255)                                                                                                                  ttGTCTTATGCCAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGa                                 <  1:429042/150‑1 (MQ=255)                                                                                                                    gTCTTATGCCAACAACGCCACGTGTGTACGGGGTGCTTTGCTTTTATATCTGTAATCTTAATGCCTCGCTGGCGATGTTACGAAAATTCCTGGAATTTGCTGGCATGTTATGCCAGTTTCATGTCAAATGGTTAATTTTTGCACAGGagt                               <  1:46072/150‑3 (MQ=255)                                                                                                                             cAAAAACGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGCGGGttt                      >  1:664175/1‑150 (MQ=255)                                                                                                                              aaaaaCGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTgg                    <  1:290938/151‑1 (MQ=255)                                                                                                                               aaaaCGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTgg                    >  1:15852/1‑150 (MQ=255)                                                                                                                                aaaCGCCACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGaa                  <  2:662059/151‑1 (MQ=255)                                                                                                                                      cACGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACt            >  1:363492/1‑151 (MQ=255)                                                                                                                                        cGTGTTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCGAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACttt          >  2:375859/1‑151 (MQ=255)                                                                                                                                         gtgtTTACGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTc         <  1:647700/151‑1 (MQ=255)                                                                                                                                               aCGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGGGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTc    <  1:480424/150‑1 (MQ=255)                                                                                                                                                cGTGGTGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTg  <  1:426827/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTG  >  NC_000913/3967566‑3967859

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