Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,230,727	C→A	intergenic (‑333/‑40)	hlyE ← / → umuD	hemolysin E/translesion error‑prone DNA polymerase V subunit; RecA‑activated auto‑protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,230,727	0	C	A	100.0%	187.3 / NA	48	intergenic (‑333/‑40)	hlyE/umuD	hemolysin E/translesion error‑prone DNA polymerase V subunit; RecA‑activated auto‑protease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (15/33);  total (15/33)

AGCCGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGT  >  NC_000913/1230580‑1230867                                                                                                                                                    |                                                                                                                                             aGCCGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGaata                                                                                                                                            >  1:694395/1‑150 (MQ=255) aGCCGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGaata                                                                                                                                            <  1:568334/150‑1 (MQ=255)    cGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGt                                                                                                                                        <  1:221659/151‑1 (MQ=255)      ggTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGtat                                                                                                                                      <  2:453620/151‑1 (MQ=255)       gTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGtat                                                                                                                                      >  1:271203/1‑150 (MQ=255)              ttttCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATaaaaa                                                                                                                               <  1:353806/150‑1 (MQ=255)                  cATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGt                                                                                                                           <  1:12548/150‑1 (MQ=255)                     ttCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGtata                                                                                                                        <  1:220977/150‑1 (MQ=255)                         gctgGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACtt                                                                                                                    <  1:599220/150‑1 (MQ=255)                         gctgGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTc                                                                                                                   >  2:107364/1‑151 (MQ=255)                            ggATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAgg                                                                                                                >  1:352401/1‑151 (MQ=255)                            ggATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAgg                                                                                                                <  1:382072/151‑1 (MQ=255)                             gATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGc                                                                                                               >  1:566660/1‑151 (MQ=255)                                    gtcgtcGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGattat                                                                                                        <  1:130782/151‑1 (MQ=255)                                       gtcgCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATtattat                                                                                                     >  1:71325/1‑151 (MQ=255)                                          gCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTAtgtt                                                                                                  <  2:304460/151‑1 (MQ=255)                                          gCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTAtgtt                                                                                                  <  2:142344/151‑1 (MQ=255)                                           ccAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTAtgttg                                                                                                 <  1:502644/151‑1 (MQ=255)                                              gAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGttt                                                                                              <  2:116182/151‑1 (MQ=255)                                              gAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGttt                                                                                              <  1:726559/151‑1 (MQ=255)                                               aaGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGttt                                                                                              >  2:2191/1‑150 (MQ=255)                                                aGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTAt                                                                                            <  2:446137/151‑1 (MQ=255)                                                    cACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAg                                                                                        <  1:311854/151‑1 (MQ=255)                                                       gTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCt                                                                                     >  1:480580/1‑151 (MQ=255)                                                        tGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTg                                                                                    >  2:100218/1‑151 (MQ=255)                                                         gAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTg                                                                                    <  1:664042/150‑1 (MQ=255)                                                         gAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGc                                                                                   <  2:684235/151‑1 (MQ=255)                                                            cacaAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCgg                                                                                 >  2:311930/1‑150 (MQ=255)                                                                  aTAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCcg                                                                          >  2:539798/1‑151 (MQ=255)                                                                           cGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATtgt                                                                 <  2:582891/151‑1 (MQ=255)                                                                           cGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATtgt                                                                 <  2:705663/151‑1 (MQ=255)                                                                             aaaaaTCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGa                                                               <  2:386266/151‑1 (MQ=255)                                                                                           tATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTAtt                                                  <  1:582703/150‑1 (MQ=255)                                                                                            aTGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTAtt                                                  <  1:110509/149‑1 (MQ=255)                                                                                                aGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGa                                            <  1:608982/151‑1 (MQ=255)                                                                                                 gCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGa                                            <  2:538134/150‑1 (MQ=255)                                                                                                 gCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGAt                                           <  1:231782/151‑1 (MQ=255)                                                                                                  cGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATCGTGACTTTTCCGCTATTTAGCGATc                                          <  2:142231/151‑1 (MQ=255)                                                                                                       aaTATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGtt                                     <  2:329326/151‑1 (MQ=255)                                                                                                          atTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCa                                   <  2:725133/150‑1 (MQ=255)                                                                                                              aTAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGtgtg                              >  2:380206/1‑151 (MQ=255)                                                                                                                aGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGc                            >  1:145078/1‑151 (MQ=255)                                                                                                                      aaTCAGTATTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCt                      <  2:284539/151‑1 (MQ=255)                                                                                                                             aTTGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCg                <  1:383000/150‑1 (MQ=255)                                                                                                                               tGATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGgc              >  2:194912/1‑150 (MQ=255)                                                                                                                                gATCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGgc              >  2:655959/1‑149 (MQ=255)                                                                                                                                 aTCTGCTGGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGgca             <  1:398018/149‑1 (MQ=255)                                                                                                                                     gctgGCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGAt       <  2:592662/151‑1 (MQ=255)                                                                                                                                        ggCAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTa     <  2:165926/150‑1 (MQ=255)                                                                                                                                          cAAGAACAGAATACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGt  >  1:691259/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                             AGCCGGGTGATATTTTTCATTTCTGCTGGATGAGCGTCGTCGCCAGAAGGCCACGTGAGCACAAGATAAGAGAACGAAAAATCAGCAGCCTATGCAGCGACAAATATTGATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGT  >  NC_000913/1230580‑1230867

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