Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	4,122,984	C→T	R119H (CGT→CAT)	ftsN ←	essential cell division protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	4,122,984	0	C	T	100.0%	182.5 / NA	55	R119H (CGT→CAT)	ftsN	essential cell division protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/27);  total (28/27)

GGCGTTGCTGTTCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGT  >  NC_000913_3_bme_pgi/4122838‑4123129                                                                                                                                                   |                                                                                                                                                  ggCGTTGCTGTTCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGCCGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAACAAGCTGATgt                                                                                                                                                 >  2:474983/1‑149 (MQ=255)     gtGCTGTTCCGCTAATTGCTGCGCGTGACGTTGGGGCTGTAGGGTTTGCTGACGCTGTTCCGGCGTCTGTTCACTCCACGGAACTTCAACCAGCTGCGTTTGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATgttgtt                                                                                                                                             <  2:235477/148‑1 (MQ=255)       gCTGTTCCGCTAATTGCTGCGCCTGATGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCt                                                                                                                                           <  2:786098/149‑1 (MQ=255)       gCTGTTCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCt                                                                                                                                           >  1:33037/1‑149 (MQ=255)        cTGTTCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTg                                                                                                                                          <  1:567999/149‑1 (MQ=255)        cTGTTCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTg                                                                                                                                          >  1:748158/1‑149 (MQ=255)          ggtCCGCTAATTGCTGTCCGGGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGgt                                                                                                                                        <  2:221233/147‑1 (MQ=255)          gTTCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGgt                                                                                                                                        >  2:429632/1‑149 (MQ=255)               gcttcttGCTGCGCCTCACGTTGGCGTTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCAGATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAg                                                                                                                                   <  1:780467/144‑1 (MQ=255)               gCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAGGCTGATGTTGTTCTAgtgtgct                                                                                                                                   >  2:607521/1‑146 (MQ=255)               gCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTCGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAg                                                                                                                                   >  2:50489/1‑149 (MQ=255)                 tAAGTGCTGCGCCTGACGTTGTCGCTGTAGCGTTTGCTGGCGCTGTTCCGGGGTCTGTTCATTCCACGGCACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCTGCCTGCATTTGTTTAAGAAGCTGATGTTGTTCTGGTGTCAGtt                                                                                                                                 <  2:52298/149‑1 (MQ=255)                   aTTGCTGCGCCTGACGTTGGCGCTATAGCGATTGATGACGCTGTTCCGGCGTCTCTCCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCCGCCTGCATTTGTTCAAGCAGCTGATTTTGTTCCGTTGACAGtttc                                                                                                                               >  2:519300/1‑147 (MQ=255)                      gctgcGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTcc                                                                                                                            >  2:352668/1‑149 (MQ=255)                        tgcGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATGAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTGGCTCCgg                                                                                                                          >  1:4510/1‑149 (MQ=255)                        tgcGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCgg                                                                                                                          >  1:485179/1‑149 (MQ=255)                         gcgcCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGc                                                                                                                         <  1:656961/149‑1 (MQ=255)                            ccTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGtt                                                                                                                      <  1:338365/149‑1 (MQ=255)                                 cGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCTAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTttcac                                                                                                                 <  2:741616/149‑1 (MQ=255)                                 cGTTGGCGCTGTAGCGTGTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGGTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTttcac                                                                                                                 <  1:309392/149‑1 (MQ=255)                                 cGTTGGCGCTGTAGCGTCTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTttcac                                                                                                                 <  1:473973/149‑1 (MQ=255)                                                 tttGCTGACGCTGTTCCGGTGTCTGTTCATGCCACGGAACTTCAACCAGCTGTGTTGGCTGCTGGCGCATATCGGCGTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGaa                                                                                                 <  2:124312/149‑1 (MQ=255)                                                 tttGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGaa                                                                                                 <  1:352668/149‑1 (MQ=255)                                                           cTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGATTCTGTgg                                                                                       <  2:706637/149‑1 (MQ=255)                                                             gTTCCGGCGTCTGTGCATTCTACGGATCTTTAACCAGGTGCGTTGGCTGTTGGCGCATATCAGCCTGAATTTGTTCAAGAAGTTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGgc                                                                                     <  2:105878/149‑1 (MQ=255)                                                                     gTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCgcacgcac                                                                             >  2:418835/1‑149 (MQ=255)                                                                     gTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCgcacgcac                                                                             >  1:276567/1‑149 (MQ=255)                                                                      tCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCaccacc                                                                                                        <  2:233451/121‑1 (MQ=255)                                                                      tCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCaccacc                                                                                                        >  1:233451/1‑121 (MQ=255)                                                                      tCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACt                                                                            >  1:94091/1‑149 (MQ=255)                                                                           tCATTCCACGGAGCTGCAACGCGGTCCGTTGGGTGTTGGCGCATATCGGCCTGCATTTGTTCAGGAAGCTGATGTTGTTCGGGTGTCAGTTGCCCCGGGGTTTTCATTTCACCACCGGCAGAAGGTTCTGTGGAcgcaggcattgttgg                                                                       <  2:97975/149‑8 (MQ=255)                                                                              ttCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTg                                                                    <  2:748158/149‑1 (MQ=255)                                                                              ttCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTCCTCCGGCGTTTTGACTTCACCACCTGCAGAAGGTTCTGTGGGCGCTCGTCCTCCCGGCtt                                                                    >  2:657756/1‑148 (MQ=255)                                                                               tCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTgg                                                                   <  1:307572/149‑1 (MQ=255)                                                                               tCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTgg                                                                   <  1:58517/149‑1 (MQ=255)                                                                                 cACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCg                                                                 >  2:745914/1‑149 (MQ=255)                                                                                  aCGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGttt                                                                                                                     >  2:554333/1‑96 (MQ=255)                                                                                  aCGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGttt                                                                                                                     <  1:554333/96‑1 (MQ=255)                                                                                  aCGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACCCCCGGCTGGCGa                                                                >  2:468198/1‑149 (MQ=255)                                                                                      gaCTTCAACGAGCTGCGTTAGCTTCCGGCGAATCTCAGCCTGCATTTGTTCAAGTATCTGATGTTGTTCTGGTGTCAGTTGCGCCGTCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACttt                                                            <  2:563944/148‑1 (MQ=255)                                                                                         ttCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCa                                                         >  1:547215/1‑149 (MQ=255)                                                                                               cAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCCATCCTGGCTGGAGAGTTTCCGGCTCtt                                                   >  1:236237/1‑149 (MQ=255)                                                                                                              ctgGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGc                                    <  1:429632/149‑1 (MQ=255)                                                                                                              ctgGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCCCCAGc                                    >  2:47343/1‑149 (MQ=255)                                                                                                               tgGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTccc                                                                         >  1:164510/1‑111 (MQ=255)                                                                                                               tgGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTccc                                                                         <  2:164510/111‑1 (MQ=255)                                                                                                                 gcgcATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGtt                                 >  1:152836/1‑149 (MQ=255)                                                                                                                     atatCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGttcttc                             <  1:745914/149‑1 (MQ=255)                                                                                                                     atatCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGttcttc                             <  1:339798/149‑1 (MQ=255)                                                                                                                          aGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGtt                        >  1:200693/1‑149 (MQ=255)                                                                                                                                 aTTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtggtgg                 >  2:394899/1‑149 (MQ=255)                                                                                                                                         aagaagCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCCGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGt         >  1:52731/1‑149 (MQ=255)                                                                                                                                            aagCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTc      <  1:409514/149‑1 (MQ=255)                                                                                                                                            aagCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTACCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTc      >  1:698275/1‑149 (MQ=255)                                                                                                                                                tGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGt  >  1:561616/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GGCGTTGCTGTTCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGT  >  NC_000913_3_bme_pgi/4122838‑4123129

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑