Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,870,024	T→C	G216G (GGT→GGC)	cbrA →	colicin M resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,870,024	0	T	C	100.0%	195.9 / NA	55	G216G (GGT→GGC)	cbrA	colicin M resistance protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/19);  total (36/19)

TAAATATGTCGCTATCCAGCAGTGGTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGTGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGAAAAATGCACGGTGCTGTTTCCCTCGCGCTGGCAGGATTT  >  NC_000913/3869892‑3870158                                                                                                                                     |                                                                                                                                       tAAATATGTCGCTATCCAGCAGTGGTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATc                                                                                                                                 <  2:679517/140‑1 (MQ=255)         tCGCTATCCAGCAGTGGTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGaaa                                                                                                                         >  2:519394/1‑140 (MQ=255)           gCTATCCAGCAGTGGTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGa                                                                                                                       >  2:150343/1‑140 (MQ=255)           gCTATCCAGCAGTGGTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGa                                                                                                                       >  1:633685/1‑140 (MQ=255)               tCCAGCAGTGGTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGt                                                                                                                   >  2:21849/1‑140 (MQ=255)               tCCAGCAGTGGTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGt                                                                                                                   >  1:746521/1‑140 (MQ=255)                       tGGTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGAcgcg                                                                                                           >  1:505035/1‑140 (MQ=255)                         gTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGAcgcg                                                                                                           >  2:803685/1‑138 (MQ=255)                         gTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGAcgcg                                                                                                           <  1:803685/138‑1 (MQ=255)                         gTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGtt                                                                                                         >  2:80449/1‑140 (MQ=255)                         gTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGtt                                                                                                         >  1:82745/1‑140 (MQ=255)                         gTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGtt                                                                                                         >  1:594432/1‑140 (MQ=255)                         gTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGCCAGACGCGtt                                                                                                         >  1:5028/1‑140 (MQ=255)                            cgcgGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGAGGCGCCTGTCCGATGAAAGACGGTCAGACGCGTTTCa                                                                                                      >  1:235040/1‑140 (MQ=255)                               ggAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCacga                                                                                                   <  2:594432/140‑1 (MQ=255)                                     aCATCCGGTGCCGTTCTACACCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGa                                                                                             <  2:165899/140‑1 (MQ=255)                                          cGGTGCCGTTCTACTCCGGCATCTTTGATAAATCGATAAGTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAagag                                                                                        <  2:309023/140‑1 (MQ=255)                                             tGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTTAGATGCTTTTCACGACGCTGAAAGAGaaa                                                                                     >  1:259647/1‑140 (MQ=255)                                                 gTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGa                                                                                 >  2:442573/1‑140 (MQ=255)                                                  ttCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAg                                                                                >  1:552343/1‑140 (MQ=255)                                                  ttCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAg                                                                                >  2:741651/1‑140 (MQ=255)                                                    cTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAgcg                                                                              <  1:80449/140‑1 (MQ=255)                                                       cTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCt                                                                           <  2:505035/140‑1 (MQ=255)                                                           tGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCa                                                                       >  2:5766/1‑140 (MQ=255)                                                            gCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGaaa                                                                                     <  2:12611/125‑1 (MQ=255)                                                            gCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGaaa                                                                                     >  1:12611/1‑125 (MQ=255)                                                            gCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAg                                                                                <  2:860622/130‑1 (MQ=255)                                                            gCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAg                                                                                >  1:860622/1‑130 (MQ=255)                                                             cATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGt                                                                     >  1:205385/1‑140 (MQ=255)                                                                  ttGATAATTCGATAACTGACTGGTATTCATGTAGTATCACCAAAGACGCTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAg                                                                <  2:259647/140‑1 (MQ=255)                                                                  ttGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGTGACAATGAGCGCCTTTCAGTTCCAg                                                                >  1:406493/1‑140 (MQ=255)                                                                  ttGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAg                                                                <  1:798479/140‑1 (MQ=255)                                                                    gATAATTCGATAACTAACTGTTTTTCATGGAGTATAAGCGAAGACGGTGATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGtt                                                              <  2:346301/140‑1 (MQ=255)                                                                       aaTTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTgg                                                           >  2:463786/1‑140 (MQ=255)                                                                        aTTCGATATCTAACTGAAATTCATGGAGTAAAAGCCAAGGCGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGt                                                          <  2:406493/140‑1 (MQ=255)                                                                            gATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGa                                                      >  2:17730/1‑140 (MQ=255)                                                                                        ttATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAATAGAAAATGGGCGCCTTACAGTTCCAGTTTGCTAAGACGGTGAAAAGCg                                          >  2:316392/1‑140 (MQ=255)                                                                                        ttATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCg                                          >  1:129075/1‑140 (MQ=255)                                                                                        ttATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAAAGAAAATGAGCGACTTCCAGTTCCACTTTGATAAGACGGTGAAAAGCg                                          >  2:296627/1‑140 (MQ=255)                                                                                          aTTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGaa                                        >  2:115728/1‑140 (MQ=255)                                                                                            tCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGaaaa                                      <  2:5028/140‑1 (MQ=255)                                                                                               tGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGAAAAATg                                   <  2:633685/140‑1 (MQ=255)                                                                                               tGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGAAAAATg                                   <  2:205385/140‑1 (MQ=255)                                                                                                ggAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGAAAAATGc                                  <  1:115728/140‑1 (MQ=255)                                                                                                 gAGTATCAGCAAAGACGATTATTTTATCTTTTGCGACGCCTATCCAATAAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGAAAAATGCa                                 <  2:235040/140‑1 (MQ=255)                                                                                                            aaGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGAAAAATGCACGGTGCTGttt                      >  1:531915/1‑140 (MQ=255)                                                                                                            aaGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGAAAAATGCACGGTGCTGttt                      >  2:517894/1‑140 (MQ=255)                                                                                                            aaGACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGAAAAATGCACGGTGCTGttt                      >  2:703003/1‑140 (MQ=255)                                                                                                              gACGGTTATTTTATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGAAAAATGCACGGTGCTGTTTcc                    >  1:538545/1‑140 (MQ=255)                                                                                                                         tATCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTGAGACGGCGAAAAGCGAGCAATGCACAGGGCGATTTCCCTCGCGCTGGc         >  2:404402/1‑140 (MQ=255)                                                                                                                           tCTTTGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGAAAAATGCACGGTGCTGTTTCCCTCGCGCTGGCAg       >  2:424564/1‑140 (MQ=255)                                                                                                                             tttGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGAAAAATGCACGGTGCTGTTTCCCTcgcg              <  1:772029/131‑1 (MQ=255)                                                                                                                             tttGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGAAAAATGCACGGTGCTGTTTCCCTcgcg              >  2:772029/1‑131 (MQ=255)                                                                                                                             tttGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGAAAAATGCACGGTGCTGTTTCCCTCGCGCTGGCAGGa     >  2:655436/1‑140 (MQ=255)                                                                                                                               tGGCGGCGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGAAAAATGCACGGTGCTGTTTCCCTCGCGCTGGCAGGAtt   >  2:694503/1‑140 (MQ=255)                                                                                                                                ggcggcGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGAATAATGCACGGTGGTGATTCCCTCGCGCTGGCAGGAttt  >  1:321739/1‑140 (MQ=255)                                                                                                                                ggcggcGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGAAAAATGCACGGTGCTGTTTCCCTCGCGCTGGCAGGAttt  >  1:69382/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                       TAAATATGTCGCTATCCAGCAGTGGTTCGCGGAGAAACATCCGGTGCCGTTCTACTCCTGCATCTTTGATAATTCGATAACTAACTGTTATTCATGGAGTATCAGCAAAGACGGTTATTTTATCTTTGGCGGTGCCTATCCAATGAAAGACGGTCAGACGCGTTTCACGACGCTGAAAGAGAAAATGAGCGCCTTTCAGTTCCAGTTTGGTAAGACGGTGAAAAGCGAAAAATGCACGGTGCTGTTTCCCTCGCGCTGGCAGGATTT  >  NC_000913/3869892‑3870158

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: ‑