Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,077,174	G→A	G21G (GGG→GGA)	cbeA →	CP4‑44 prophage; cytoskeleton bundling‑enhancing antitoxin CbeA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,077,174	0	G	A	100.0%	94.1 / NA	32	G21G (GGG→GGA)	cbeA	CP4‑44 prophage; cytoskeleton bundling‑enhancing antitoxin CbeA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/10);  total (22/10)

AGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGG  >  NC_000913/2077038‑2077304                                                                                                                                         |                                                                                                                                   cgTGGGAATTCTGATGAGCAGAATTATGGCCACTACCTCACCCTTCATAGCCGATTTTCTGTAAAGATTTTAGCGTGCCAGACACACACCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTg                                                                                                                                 >  2:332879/2‑140 (MQ=255) aGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTg                                                                                                                                 >  2:411374/1‑140 (MQ=255) aGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTg                                                                                                                                 >  1:603749/1‑140 (MQ=255)  gTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGc                                                                                                                                >  2:571921/1‑140 (MQ=255)  gTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGc                                                                                                                                >  1:366903/1‑140 (MQ=255)      gAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGCCCTg                                                                                                                            >  1:365859/1‑140 (MQ=255)             gatgaGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAGTCACGACGGCCTCTAGCGGGGACTGCCCTGCACCGTg                                                                                                                     >  1:198926/1‑140 (MQ=255)                  gCAGAATTATCGCCACTCCCGGCCCATTCTTTGTCAATTGTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGcc                                                                                                                <  2:306665/140‑1 (MQ=255)                      aaTTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGt                                                                                                            >  1:336807/1‑140 (MQ=255)                           tCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCgg                                                                                                       >  2:392559/1‑140 (MQ=255)                                aCTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCAc                                                                                                  >  1:33603/1‑140 (MQ=255)                                  tACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGAATACGCTGCACCGAGACGCCCTGTTTCGGGGCACGt                                                                                                >  2:296536/1‑140 (MQ=255)                                         ccATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGc                                                                                         <  1:475486/140‑1 (MQ=255)                                            ttCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGGCCTGTTTCGGGGCACGCCTGGTGCagg                                                                                      >  2:369847/1‑140 (MQ=255)                                            ttCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGCGACTGCCCGGCACCGTGACGCCCTGCTTCGCGGCACGTCTGGTGCagg                                                                                      >  2:236869/1‑140 (MQ=255)                                                        tttCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGt                                                                          >  2:466050/1‑140 (MQ=255)                                                                 ggATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTAcc                                                                 <  1:392559/140‑1 (MQ=255)                                                                         tCGTGTCAGACACACCCCCCGGGACAACACTTCCCGAAGACAATCACGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGAc                                                         <  1:296536/140‑1 (MQ=255)                                                                             gtCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACcgcg                                                     >  1:432048/1‑140 (MQ=255)                                                                             gtCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTACAGGAGGATAAGCGCTTGCATCACCTTGCCGACcgcg                                                     >  1:222950/1‑140 (MQ=255)                                                                                        tCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCaga                                          >  2:391135/1‑140 (MQ=255)                                                                                               gACAACACTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGt                                   >  2:398723/1‑140 (MQ=255)                                                                                                acaacaCTTCCCGACTACAATCCCGCCCGTCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGtt                                  <  1:223784/140‑1 (MQ=255)                                                                                                 caacaCTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTc                                 <  1:391135/140‑1 (MQ=255)                                                                                                   acacTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAg                               <  2:603749/140‑1 (MQ=255)                                                                                                    cacTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCTACCACGGCGGTATCAGAGGCCTGTTCAGc                              >  1:245053/1‑140 (MQ=255)                                                                                                    cacTTCCCGACGACAATCACGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGc                              >  2:457798/1‑140 (MQ=255)                                                                                                         cccGACGACAATCACGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGc                         >  2:152656/1‑140 (MQ=255)                                                                                                                 cAATCGCGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGt                 <  2:65824/140‑1 (MQ=255)                                                                                                                       cGACCGCCCCTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTaccacc           <  2:336807/140‑1 (MQ=255)                                                                                                                         aCCGCGCCTGGCGGGGACTGCCCTGCACCGTGACGCCATGTTTCGGGGCACGTCTGGGGCAGGAGGGTAACCGGTTGCATTACCTTGCCGAGCGCGCCGGCATCAGAGGCCTGTTCACCGATGCAGATGCGCACCACCTg         >  2:273060/1‑140 (MQ=255)                                                                                                                             ccccTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGAcc     >  2:17156/1‑140 (MQ=255)                                                                                                                                cTGGCGGGGACTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAgg  <  1:152656/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                   AGTGGGAATTCTGATGAGCAGAATTATCGCCACTACCGGACCATTCTTAGCCGATTTTCTGTAAGGATTTTATCGTGTCAGACACACTCCCCGGGACAACACTTCCCGACGACAATCACGACCGCCCCTGGTGGGGGCTGCCCTGCACCGTGACGCCCTGTTTCGGGGCACGTCTGGTGCAGGAGGGTAACCGGTTGCATTACCTTGCCGACCGCGCCGGTATCAGAGGCCTGTTCAGCGATGCAGATGCGTACCACCTGGACCAGG  >  NC_000913/2077038‑2077304

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