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%	201.3 / NA	62	*476* (TAG→TAA)	rmuC	DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/32);  total (30/32)

AAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACTTTG  >  NC_000913/4018626‑4018888                                                                                                                                       |                                                                                                                                 aaaCGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCa                                                                                                                              <  2:1687486/139‑1 (MQ=255)  aaCGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAt                                                                                                                             <  2:1022114/139‑1 (MQ=255)   aCGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAtt                                                                                                                            <  1:60642/139‑1 (MQ=255)   aCGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAtt                                                                                                                            <  1:815157/139‑1 (MQ=255)   aCGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAtt                                                                                                                            >  2:815157/1‑139 (MQ=255)                   gATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCg                                                                                                            >  2:899690/1‑139 (MQ=255)                   gATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCg                                                                                                            >  1:19815/1‑139 (MQ=255)                        ggCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAg                                                                                                       <  1:899690/139‑1 (MQ=255)                        ggCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAg                                                                                                       <  2:149420/139‑1 (MQ=255)                         gCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGtt                                                                                                                  >  2:1640929/1‑127 (MQ=255)                         gCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGtt                                                                                                                  <  1:1640929/127‑1 (MQ=255)                             aaCAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATc                                                                                                  >  2:1566491/1‑139 (MQ=255)                                  gCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAg                                                                                                               >  1:605844/1‑121 (MQ=255)                                  gCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAg                                                                                                               <  2:605844/121‑1 (MQ=255)                                       gAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGa                                                                                                                        >  1:1670742/1‑107 (MQ=255)                                       gAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGa                                                                                                                        <  2:1670742/107‑1 (MQ=255)                                            aGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                   <  2:19815/139‑1 (MQ=255)                                                            aCTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTc                                                                   >  2:616694/1‑139 (MQ=255)                                                               tcggtcggTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTa                                                                                                        <  2:949162/99‑1 (MQ=255)                                                               tcggtcggTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTa                                                                                                        >  1:949162/1‑99 (MQ=255)                                                                                agcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                               <  1:1566491/139‑1 (MQ=255)                                                                                 gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAAccc                                                                                                                               >  1:674384/1‑58 (MQ=255)                                                                                 gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAAccc                                                                                                                               <  2:674384/58‑1 (MQ=255)                                                                                 gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGcgtag                                                                                                                    >  1:53495/1‑69 (MQ=255)                                                                                 gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGagtag                                                                                                                    <  2:53495/69‑1 (MQ=255)                                                                                 gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGaa                                                                                                     >  1:1256784/1‑84 (MQ=255)                                                                                 gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGaa                                                                                                     <  2:1256784/84‑1 (MQ=255)                                                                                 gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                              <  1:448389/139‑1 (MQ=255)                                                                                    gAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAgg                                           <  1:616694/139‑1 (MQ=255)                                                                                          tgaAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTa                                                                         <  2:684760/103‑1 (MQ=255)                                                                                          tgaAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTa                                                                         >  1:684760/1‑103 (MQ=255)                                                                                            aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATc                                                                                                  <  2:129721/76‑1 (MQ=255)                                                                                            aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATc                                                                                                  >  1:129721/1‑76 (MQ=255)                                                                                               cTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGtt                                                                                                                  <  1:1637177/57‑1 (MQ=255)                                                                                               cTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGtt                                                                                                                  >  2:1637177/1‑57 (MQ=255)                                                                                                 tATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGAtggtg                              >  2:854130/1‑139 (MQ=255)                                                                                                 tATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGAtggtg                              >  2:472790/1‑139 (MQ=255)                                                                                                 tATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGAtggtg                              >  2:328946/1‑139 (MQ=255)                                                                                                  aTCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAt                                                                                                                             >  1:250003/1‑43 (MQ=255)                                                                                                  aTCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAt                                                                                                                             <  2:250003/43‑1 (MQ=255)                                                                                                   tCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                   <  2:1047939/84‑1 (MQ=255)                                                                                                   tCAACGCGATGATGAATATAATCAGCAGACGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                   >  1:1047939/1‑84 (MQ=255)                                                                                                    cAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGAt                           <  1:472790/139‑1 (MQ=255)                                                                                                         cgATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCa                                         <  2:225166/120‑1 (MQ=255)                                                                                                         cgATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCa                                         >  1:225166/1‑120 (MQ=255)                                                                                                            tgatgaATATAATCAGCAGTCGCGCCAACCCATTggg                                                                                                                         >  2:764219/1‑37 (MQ=255)                                                                                                            tgatgaATATAATCAGCAGTCGCGCCAACCCATTggg                                                                                                                         <  1:764219/37‑1 (MQ=255)                                                                                                                   tataATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAg                                                                                                       <  1:1269493/48‑1 (MQ=255)                                                                                                                   tataATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAg                                                                                                       >  2:1269493/1‑48 (MQ=255)                                                                                                                    ataATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAacga           <  1:328946/139‑1 (MQ=255)                                                                                                                         cagcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCAc      >  2:1164804/1‑139 (MQ=255)                                                                                                                           gcagTCGCGCTAACCTATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTgaga                                   >  2:1131038/1‑108 (MQ=255)                                                                                                                           gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCaa                                                                                >  1:1048841/1‑63 (MQ=255)                                                                                                                           gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCaa                                                                                <  2:1048841/63‑1 (MQ=255)                                                                                                                           gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAAc                                                             >  2:314083/1‑82 (MQ=255)                                                                                                                           gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAAc                                                             <  1:314083/82‑1 (MQ=255)                                                                                                                           gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTgaga                                   <  1:1131038/108‑1 (MQ=255)                                                                                                                             agTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACTTTg  >  1:1406462/1‑139 (MQ=255)                                                                                                                               tCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTa                                                                         >  1:960548/1‑66 (MQ=255)                                                                                                                               tCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTa                                                                         <  2:960548/66‑1 (MQ=255)                                                                                                                               tCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                               <  1:941285/92‑1 (MQ=255)                                                                                                                               tCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                               >  2:941285/1‑92 (MQ=255)                                                                                                                                       |                                                                                                                                 AAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACTTTG  >  NC_000913/4018626‑4018888

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