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

CGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACTTTGGT  >  NC_000913/4018631‑4018890                                                                                                                                  |                                                                                                                                   cgAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTggg                                                                                                                           >  1:147950/1‑139 (MQ=255)             ttaTTTGCCTGAACAGGCGGTGAGTCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCTAATGATGAAGCTTATCAACGCGATGATGAATATTATCACTAGTCGCTCTAACCCATTGGGAGTAGTTAAGcc                                                                                                               <  2:707232/137‑1 (MQ=255)             ggATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGcc                                                                                                               <  2:1037298/139‑1 (MQ=255)                                                    tCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTggg                                                                                                                           <  2:349078/88‑1 (MQ=255)                                                    tCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTggg                                                                                                                           >  1:349078/1‑88 (MQ=255)                                                       aCTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTc                                                                     >  1:126934/1‑139 (MQ=255)                                                       aCTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTc                                                                     >  2:29293/1‑139 (MQ=255)                                                            ggtcggtTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGaa                                                                >  2:1101826/1‑139 (MQ=255)                                                                            gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAt                                                                                                                               >  2:963998/1‑60 (MQ=255)                                                                            gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAt                                                                                                                               <  1:963998/60‑1 (MQ=255)                                                                            gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCggg                                                                                                            <  1:953266/79‑1 (MQ=255)                                                                            gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCggg                                                                                                            >  2:953266/1‑79 (MQ=255)                                                                            gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                                 <  2:1114781/138‑1 (MQ=255)                                                                            gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                                 >  1:1114781/1‑138 (MQ=255)                                                                            gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                                <  1:29293/139‑1 (MQ=255)                                                                                aaTGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCg                                                                                                              >  1:309818/1‑73 (MQ=255)                                                                                aaTGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCg                                                                                                              <  2:309818/73‑1 (MQ=255)                                                                                       aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCggg                                                                                                            >  2:1031406/1‑68 (MQ=255)                                                                                       aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCggg                                                                                                            <  1:1031406/68‑1 (MQ=255)                                                                                       aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAAt                                                                                 <  2:443466/95‑1 (MQ=255)                                                                                       aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAAt                                                                                 >  1:443466/1‑95 (MQ=255)                                                                                       aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTgaga                                     <  2:147950/139‑1 (MQ=255)                                                                                       aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTgaga                                     >  1:977292/1‑139 (MQ=255)                                                                                              tCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCTAATCTGTTACACTTCTGGaa                                                                <  2:667631/105‑1 (MQ=255)                                                                                              tCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGaa                                                                >  1:667631/1‑105 (MQ=255)                                                                                              tCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTtgatgatc                                                <  2:73607/121‑3 (MQ=255)                                                                                              tCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                                >  1:73607/1‑121 (MQ=255)                                                                                               cAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGAt                             <  2:62400/139‑1 (MQ=255)                                                                                               cAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGAt                             <  1:1101826/139‑1 (MQ=255)                                                                                                   gcgATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAg                                                                                                         >  2:155849/1‑59 (MQ=255)                                                                                                   gcgATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAg                                                                                                         <  1:155849/59‑1 (MQ=255)                                                                                                      atgatgAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTgaga                                     <  1:304027/124‑1 (MQ=255)                                                                                                      atgatgAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTgaga                                     >  2:304027/1‑124 (MQ=255)                                                                                                             atatAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                                 >  2:288414/1‑105 (MQ=255)                                                                                                             atatAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                                 <  1:288414/105‑1 (MQ=255)                                                                                                                     agcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACt       >  1:521308/1‑139 (MQ=255)                                                                                                                      gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAg                                                                                                         >  1:527434/1‑40 (MQ=255)                                                                                                                      gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAg                                                                                                         <  2:527434/40‑1 (MQ=255)                                                                                                                      gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGa                              >  1:316480/1‑115 (MQ=255)                                                                                                                      gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGa                              <  2:316480/115‑1 (MQ=255)                                                                                                                      gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGAt                             <  1:1003865/116‑1 (MQ=255)                                                                                                                      gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGAt                             >  2:1003865/1‑116 (MQ=255)                                                                                                                      gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGCGCATGGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATCGAGATGGTCGATAACTCACAAGAAACGACGCACtt      >  2:395002/1‑139 (MQ=255)                                                                                                                       cagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCaa                                                                                  <  2:6072/62‑1 (MQ=255)                                                                                                                       cagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCaa                                                                                  >  1:6072/1‑62 (MQ=255)                                                                                                                         gTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                     >  2:131897/1‑57 (MQ=255)                                                                                                                         gTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                     <  1:131897/57‑1 (MQ=255)                                                                                                                          tCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACTTTGGt  >  1:968236/1‑139 (MQ=255)                                                                                                                            gcgcTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                                <  2:299071/91‑1 (MQ=255)                                                                                                                            gcgcTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                                >  1:299071/1‑91 (MQ=255)                                                                                                                                  |                                                                                                                                   CGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACTTTGGT  >  NC_000913/4018631‑4018890

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