Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,869,607	C→A	S77R (AGC→AGA)	cbrA →	colicin M resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,869,607	0	C	A	100.0%	126.9 / NA	37	S77R (AGC→AGA)	cbrA	colicin M resistance protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (18/19);  total (18/19)

TGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGT  >  NC_000913/3869485‑3869743                                                                                                                           |                                                                                                                                         tgtgGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGa                                                                                                                            >  1:211352/1‑137 (MQ=255) tgtgGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGa                                                                                                                            <  2:211352/137‑1 (MQ=255) tgtgGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGt                                                                                                                         <  2:474971/140‑1 (MQ=255) tgtgGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGt                                                                                                                         >  2:259941/1‑140 (MQ=255)  gtgGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTc                                                                                                                        >  1:660620/1‑140 (MQ=255)        cTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCa                                                                                                                  <  1:202291/140‑1 (MQ=255)         tGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCAt                                                                                                                 <  1:97597/140‑1 (MQ=255)           aaGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCg                                                                                                               >  2:646529/1‑140 (MQ=255)           aaGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCg                                                                                                               >  2:439502/1‑140 (MQ=255)            aGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGc                                                                                                              <  1:520921/140‑1 (MQ=255)             ggTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCt                                                                                                             <  2:179205/140‑1 (MQ=255)               tttCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCa                                                                                                           <  2:660620/140‑1 (MQ=255)                ttCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCAc                                                                                                          >  2:786620/1‑140 (MQ=255)                ttCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCAc                                                                                                          >  1:313579/1‑140 (MQ=255)                ttCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCAc                                                                                                          <  1:665318/140‑1 (MQ=255)                ttCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCAc                                                                                                          <  1:259941/140‑1 (MQ=255)                ttCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCAc                                                                                                          >  2:665318/1‑140 (MQ=255)                  cAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCAcgc                                                                                                        >  2:558043/1‑140 (MQ=255)                  cAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCAcgc                                                                                                        <  1:646529/140‑1 (MQ=255)                             gcggcgGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAg                                                                                             <  1:234223/140‑1 (MQ=255)                             gcggcgGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAg                                                                                             >  2:754502/1‑140 (MQ=255)                              cggcggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGc                                                                                            <  2:251875/140‑1 (MQ=255)                                 cggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGaa                                                                                         <  1:733946/140‑1 (MQ=255)                                      tgctgGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCtat                                                                                    <  1:558043/140‑1 (MQ=255)                                                  aCGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAAt                                                                        >  1:277951/1‑140 (MQ=255)                                                         gCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATg                                                                 <  1:786620/140‑1 (MQ=255)                                                                    ttCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTtg                                                      >  1:791894/1‑140 (MQ=255)                                                                        cgATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGa                                                  <  2:277951/140‑1 (MQ=255)                                                                        cgATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGa                                                  >  2:626159/1‑140 (MQ=255)                                                                                       tCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATcc                                   >  2:261099/1‑140 (MQ=255)                                                                                        ccTGTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATccc                                  >  2:262862/1‑140 (MQ=255)                                                                                           gTCGATGTGATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGc                               <  1:262862/140‑1 (MQ=255)                                                                                                  tgATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGtt                        <  2:791894/140‑1 (MQ=255)                                                                                                  tgATCGCTAATCCGCAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGtt                        >  2:41548/1‑140 (MQ=255)                                                                                                                cAGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGa          >  2:135856/1‑140 (MQ=255)                                                                                                                 aGATTTTCAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGAc         >  1:131649/1‑140 (MQ=255)                                                                                                                        cAGAGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGTCTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCtgt  <  1:754502/140‑1 (MQ=255)                                                                                                                           |                                                                                                                                         TGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGT  >  NC_000913/3869485‑3869743

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