Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,898,149	C→T	S43L (TCG→TTG)	cbrC →	UPF0167 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,898,149	0	C	T	100.0%	201.2 / NA	63	S43L (TCG→TTG)	cbrC	UPF0167 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/29);  total (34/29)

TGTTTCTATGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTGAATTTGA  >  NC_000913/3898015‑3898275                                                                                                                                       |                                                                                                                               tGTTTCTATGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGt                                                                                                                            <  1:416511/139‑1 (MQ=255) tGTTTCTATGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGt                                                                                                                            <  2:303660/139‑1 (MQ=255)         tGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTgg                                                                                                                              <  1:552338/129‑1 (MQ=255)         tGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTgg                                                                                                                              >  2:552338/1‑129 (MQ=255)               aaaaTATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGt                                                                                                              >  2:250854/1‑139 (MQ=255)                  atatCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTtga                                                                                                           <  1:223746/139‑1 (MQ=255)                    atCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTtgatg                                                                                                         <  2:753626/139‑1 (MQ=255)                         gCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGtt                                                                                                    >  2:712761/1‑139 (MQ=255)                           cGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGa                                                                                                  >  1:735626/1‑139 (MQ=255)                           cGTTACACCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGa                                                                                                  >  1:679062/1‑139 (MQ=255)                            gTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGaa                                                                                                 <  2:546419/139‑1 (MQ=255)                                 cccAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACAtct                                                                                            <  2:28271/139‑1 (MQ=255)                                   cAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACAtct                                                                                            <  2:179343/137‑1 (MQ=255)                                   cAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACAtct                                                                                            >  1:179343/1‑137 (MQ=255)                                   cAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATctct                                                                                          <  2:12372/139‑1 (MQ=255)                                     aTTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGt                                                                                        <  1:712761/139‑1 (MQ=255)                                     aTTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGt                                                                                        >  2:139483/1‑139 (MQ=255)                                            tatCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTCAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCgtggt                                                                                 >  1:764018/1‑139 (MQ=255)                                               catcCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGt                                                                                                                            >  1:24302/1‑93 (MQ=255)                                               catcCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGt                                                                                                                            <  2:24302/93‑1 (MQ=255)                                                atcCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTAt                                                                             >  1:22487/1‑139 (MQ=255)                                                atcCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTAt                                                                             >  1:686711/1‑139 (MQ=255)                                                 tcCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTcc                                                                                      >  1:352895/1‑129 (MQ=255)                                                 tcCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTcc                                                                                      <  2:352895/129‑1 (MQ=255)                                                 tcCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCgtggt                                                                                 >  2:527798/1‑134 (MQ=255)                                                 tcCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCgtggt                                                                                 <  1:527798/134‑1 (MQ=255)                                                 tcCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTAtt                                                                            >  2:571783/1‑139 (MQ=255)                                                    cAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCg                                                                         <  2:114622/139‑1 (MQ=255)                                                          aCTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGt                                                                   >  1:693632/1‑139 (MQ=255)                                                            tGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGt                                                                                                                            >  2:791085/1‑80 (MQ=255)                                                            tGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGt                                                                                                                            <  1:791085/80‑1 (MQ=255)                                                            tGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTc                                                                 >  1:152444/1‑139 (MQ=255)                                                                  cAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTctgctgc                                                           >  1:203096/1‑139 (MQ=255)                                                                             gAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTg                                                >  2:262670/1‑139 (MQ=255)                                                                                 aGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAgg                                            <  1:250854/139‑1 (MQ=255)                                                                                   gATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTa                                          >  1:117404/1‑139 (MQ=255)                                                                                             tAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTtgatga                                                                                                        <  2:506933/67‑1 (MQ=255)                                                                                             tAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTtgatga                                                                                                        >  1:506933/1‑67 (MQ=255)                                                                                             tAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAg                                                                                                      <  2:412124/69‑1 (MQ=255)                                                                                             tAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAg                                                                                                      >  1:412124/1‑69 (MQ=255)                                                                                                 gTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGgatgat                            <  2:679062/139‑1 (MQ=255)                                                                                                   gCGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGc                          <  1:49374/139‑1 (MQ=255)                                                                                                    cGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGcc                         <  2:693632/139‑1 (MQ=255)                                                                                                    cGATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGcc                         <  2:686711/139‑1 (MQ=255)                                                                                                     gATTGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCa                        >  2:286839/1‑139 (MQ=255)                                                                                                       ttGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTtgatga                                                                                                        >  2:275657/1‑57 (MQ=255)                                                                                                       ttGCTGTGAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTtgatga                                                                                                        <  1:275657/57‑1 (MQ=255)                                                                                                            gtgAACAACAGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGa                 >  2:240082/1‑139 (MQ=255)                                                                                                                    caGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTg         >  1:763537/1‑139 (MQ=255)                                                                                                                    caGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTg         >  2:723135/1‑139 (MQ=255)                                                                                                                    caGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTg         <  2:203096/139‑1 (MQ=255)                                                                                                                    caGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTg         >  2:63659/1‑139 (MQ=255)                                                                                                                    caGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTg         >  2:479970/1‑139 (MQ=255)                                                                                                                    caGACGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTg         >  2:397144/1‑139 (MQ=255)                                                                                                                       aCGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGAc                                                                      >  2:442326/1‑75 (MQ=255)                                                                                                                       aCGTCAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGAc                                                                      <  1:442326/75‑1 (MQ=255)                                                                                                                           cAGTTTATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTGAATTTGa  <  1:262670/139‑1 (MQ=255)                                                                                                                              tttATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGa                                                                       <  2:39000/67‑1 (MQ=255)                                                                                                                              tttATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGa                                                                       >  1:39000/1‑67 (MQ=255)                                                                                                                              tttATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGgat                               <  2:215445/107‑1 (MQ=255)                                                                                                                              tttATTACTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGgat                               >  1:215445/1‑107 (MQ=255)                                                                                                                                   taCTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTcc                                                                                      >  2:64279/1‑47 (MQ=255)                                                                                                                                   taCTTGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTcc                                                                                      <  1:64279/47‑1 (MQ=255)                                                                                                                                       |                                                                                                                               TGTTTCTATGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTGAATTTGA  >  NC_000913/3898015‑3898275

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