Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,018,760	G→A	*476* (TAG→TAA)	rmuC →	DNA recombination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,018,760	0	G	A	100.0%	219.8 / NA	67	*476* (TAG→TAA)	rmuC	DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/34);  total (33/34)

GATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACTTTG  >  NC_000913/4018634‑4018888                                                                                                                               |                                                                                                                                 gaTTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGagt                                                                                                                      >  1:718759/1‑139 (MQ=255)           gATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCg                                                                                                            >  2:757684/1‑139 (MQ=255)            aTTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCgg                                                                                                           >  1:1397533/1‑139 (MQ=255)             tttGTCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCggg                                                                                                          <  2:614555/139‑1 (MQ=255)                ggCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAg                                                                                                       <  2:1104483/139‑1 (MQ=255)                     aaCAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATc                                                                                                  >  1:1044721/1‑139 (MQ=255)                              tGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATc                                                                                         <  1:757684/139‑1 (MQ=255)                                              gTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTa                                                                         <  2:718759/139‑1 (MQ=255)                                                aTCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCCGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTaca                                                                       >  2:767723/1‑139 (MQ=255)                                                aTCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTaca                                                                       >  1:1309496/1‑139 (MQ=255)                                                aTCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTaca                                                                       >  2:1464767/1‑139 (MQ=255)                                                    aCTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTc                                                                   >  1:131829/1‑139 (MQ=255)                                                       tcggtcggTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAt                                                                                                                             <  2:727993/78‑1 (MQ=255)                                                       tcggtcggTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAt                                                                                                                             >  1:727993/1‑78 (MQ=255)                                                                    gAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTtga                                                   <  2:131829/139‑1 (MQ=255)                                                                        agcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCCGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                               <  1:767723/139‑1 (MQ=255)                                                                         gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTCAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGTACAATTTTTTGATGAGc                                              <  2:918160/139‑1 (MQ=255)                                                                         gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCg                                                                                                            <  2:831643/77‑1 (MQ=255)                                                                         gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCg                                                                                                            >  1:831643/1‑77 (MQ=255)                                                                         gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAAt                                                                               <  1:318226/106‑1 (MQ=255)                                                                         gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAAt                                                                               >  2:318226/1‑106 (MQ=255)                                                                         gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                              >  1:280889/1‑139 (MQ=255)                                                                         gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                              <  2:1464827/139‑1 (MQ=255)                                                                         gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                              <  1:907661/139‑1 (MQ=255)                                                                                    aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATc                                                                                         >  2:992262/1‑85 (MQ=255)                                                                                    aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATc                                                                                         <  1:992262/85‑1 (MQ=255)                                                                                    aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTgaga                                   <  2:19893/139‑1 (MQ=255)                                                                                       cTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTggg                                                                                                                         <  2:660181/50‑1 (MQ=255)                                                                                       cTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTggg                                                                                                                         >  1:660181/1‑50 (MQ=255)                                                                                       cTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAggg                                                                                             >  1:763046/1‑78 (MQ=255)                                                                                       cTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAggg                                                                                             <  2:763046/78‑1 (MQ=255)                                                                                          aTCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAt                                                                                                                             <  2:484447/43‑1 (MQ=255)                                                                                          aTCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAt                                                                                                                             >  1:484447/1‑43 (MQ=255)                                                                                           tCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATc                                                                              >  1:121306/1‑89 (MQ=255)                                                                                           tCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATc                                                                              <  2:121306/89‑1 (MQ=255)                                                                                           tCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTg                                      >  1:77346/1‑129 (MQ=255)                                                                                           tCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTg                                      <  2:77346/129‑1 (MQ=255)                                                                                            cAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACTCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGAt                           <  1:236908/139‑1 (MQ=255)                                                                                               cgcgATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAggg                                                                                             >  1:891654/1‑70 (MQ=255)                                                                                               cgcgATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAggg                                                                                             <  2:891654/70‑1 (MQ=255)                                                                                                    tgatgaATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                               <  2:928327/111‑1 (MQ=255)                                                                                                    tgatgaATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                               >  1:928327/1‑111 (MQ=255)                                                                                                    tgatgaATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTcaca                   >  2:1029050/1‑139 (MQ=255)                                                                                                            ataATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAacga           <  2:1397533/139‑1 (MQ=255)                                                                                                            ataATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAacga           <  2:1309496/139‑1 (MQ=255)                                                                                                            ataATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAacga           <  1:449628/139‑1 (MQ=255)                                                                                                              aaTCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                   >  2:1441531/1‑65 (MQ=255)                                                                                                              aaTCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                   <  1:1441531/65‑1 (MQ=255)                                                                                                                tCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCa       >  2:1259616/1‑139 (MQ=255)                                                                                                                  agcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTg                                      >  2:369796/1‑106 (MQ=255)                                                                                                                  agcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTg                                      <  1:369796/106‑1 (MQ=255)                                                                                                                   gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGt                                                                                                         <  1:135208/38‑1 (MQ=255)                                                                                                                   gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGt                                                                                                         >  2:135208/1‑38 (MQ=255)                                                                                                                   gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                               >  2:303644/1‑96 (MQ=255)                                                                                                                   gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                               <  1:303644/96‑1 (MQ=255)                                                                                                                   gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                              <  2:1289184/97‑1 (MQ=255)                                                                                                                   gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                              >  1:1289184/1‑97 (MQ=255)                                                                                                                   gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACtt    >  1:468862/1‑139 (MQ=255)                                                                                                                    cagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCa                                                                                 <  2:874743/61‑1 (MQ=255)                                                                                                                    cagTCGCGCTAACCCATTGGGAGTAGTAAAGCAGGGTAGAAATCTAGGGCATCGACGCCCa                                                                                 >  1:874743/1‑61 (MQ=255)                                                                                                                     agTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCa                                                                                 <  1:1183898/60‑1 (MQ=255)                                                                                                                     agTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCa                                                                                 >  2:1183898/1‑60 (MQ=255)                                                                                                                     agTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCACCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACTTTg  >  2:99315/1‑139 (MQ=255)                                                                                                                       tCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                   <  2:205879/56‑1 (MQ=255)                                                                                                                       tCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                   >  1:205879/1‑56 (MQ=255)                                                                                                                         gcgcTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCa                                         <  2:568095/96‑1 (MQ=255)                                                                                                                         gcgcTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCa                                         >  1:568095/1‑96 (MQ=255)                                                                                                                               |                                                                                                                                 GATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACTTTG  >  NC_000913/4018634‑4018888

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