Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,865,432	G→T	S281* (TCG→TAG)	yidE ←	putative transport protein YidE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,865,432	0	G	T	100.0%	188.9 / NA	54	S281* (TCG→TAG)	yidE	putative transport protein YidE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (30/24);  total (30/24)

CCCGGCACGGTTCAGGCGCGAGATAACAACGTCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCG  >  NC_000913/3865296‑3865559                                                                                                                                         |                                                                                                                                cccGGCACGGTTCAGGCGCGAGATAACAACGTCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAg                                                                                                                              <  1:477275/140‑1 (MQ=255)   cGGCACGGTTCAGGCGCGAGATAACAACGTCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTa                                                                                                                            >  1:582392/1‑140 (MQ=255)     gCACGGTTCAGGCGCGAGATAACAACGTCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATc                                                                                                                          >  2:48151/1‑140 (MQ=255)     gCACGGTTCAGGCGCGAGATAACAACGTCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATc                                                                                                                          >  2:415380/1‑140 (MQ=255)          gTTCAGGCGCGAGATAACAACGTCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGAcct                                                                                                                     <  2:282340/140‑1 (MQ=255)          gTTCAGGCGCGAGATAACAACGTCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGAcct                                                                                                                     >  1:666489/1‑140 (MQ=255)               ggCGCGAGATAACAACGTCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGa                                                                                                                >  1:382847/1‑140 (MQ=255)                  gcgAGATAACAACGTCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCa                                                                                                             <  2:48577/140‑1 (MQ=255)                   cgAGATAACAACGTCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCaa                                                                                                            >  2:88672/1‑140 (MQ=255)                   cgAGATAACAACGTCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCaa                                                                                                            >  2:83132/1‑140 (MQ=255)                     agaTAACAACGTCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATc                                                                                                                          >  2:111014/1‑124 (MQ=255)                     agaTAACAACGTCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATc                                                                                                                          <  1:111014/124‑1 (MQ=255)                     agaTAACAACGTCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATc                                                                                                          >  1:604932/1‑140 (MQ=255)                     agaTAACAACGTCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATc                                                                                                          >  1:270928/1‑140 (MQ=255)                       aTAACAACGTCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCAc                                                                                                        <  1:657587/140‑1 (MQ=255)                              cGTCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTgc                                                                                                 >  2:540274/1‑140 (MQ=255)                                tCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTgcgc                                                                                               <  1:484054/140‑1 (MQ=255)                                tCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTgcgc                                                                                               >  2:660939/1‑140 (MQ=255)                                 cATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCa                                                                                              <  2:270928/140‑1 (MQ=255)                                   gaGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGGACGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCAtt                                                                                            >  2:583288/2‑140 (MQ=255)                                              aaaGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCg                                                                                 <  1:380432/140‑1 (MQ=255)                                                aGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCgct                                                                               >  2:107025/1‑140 (MQ=255)                                                  tGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCgctgg                                                                             >  1:105800/1‑140 (MQ=255)                                                   gCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCgctggc                                                                            <  2:666489/140‑1 (MQ=255)                                                    cAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCgctggct                                                                           >  2:225641/1‑140 (MQ=255)                                                    cAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCgctggct                                                                           >  2:77674/1‑140 (MQ=255)                                                     aGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCgctggctg                                                                          >  2:398205/1‑140 (MQ=255)                                                       gTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGAc                                                                        >  2:354929/1‑140 (MQ=255)                                                          gcgAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACcca                                                                     <  2:223550/140‑1 (MQ=255)                                                            gAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACccacc                                                                   <  1:191702/140‑1 (MQ=255)                                                                   ttttCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATgca                                                            >  2:279881/1‑140 (MQ=255)                                                                   ttttCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATgca                                                            <  2:726934/140‑1 (MQ=255)                                                                        cGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCaaa                                                       >  2:112249/1‑140 (MQ=255)                                                                           gCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCg                                                    <  1:354929/140‑1 (MQ=255)                                                                             aCGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGcc                                                  <  1:438302/140‑1 (MQ=255)                                                                             aCGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGcc                                                  <  1:279881/140‑1 (MQ=255)                                                                             aCGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGcc                                                  <  2:393657/140‑1 (MQ=255)                                                                                ttttCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAg                                               <  1:345740/140‑1 (MQ=255)                                                                                ttttCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAg                                               <  2:618246/140‑1 (MQ=255)                                                                                ttttCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAg                                               <  2:582392/140‑1 (MQ=255)                                                                                     aTTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGa                                          <  1:775106/140‑1 (MQ=255)                                                                                         gTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAAt                                      >  2:316336/1‑140 (MQ=255)                                                                                            accaccACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGt                                   >  2:749191/1‑140 (MQ=255)                                                                                                    acGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCg                           <  1:660939/140‑1 (MQ=255)                                                                                                    acGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATAGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCg                           >  2:59546/1‑140 (MQ=255)                                                                                                        tCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAgg                       <  1:583288/140‑1 (MQ=255)                                                                                                        tCGACGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAgg                       >  1:152629/1‑140 (MQ=255)                                                                                                           aCGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAAc                    >  2:196048/1‑140 (MQ=255)                                                                                                           aCGCGCAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAAc                    >  1:405684/1‑140 (MQ=255)                                                                                                                cAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAgg                       <  1:476852/132‑1 (MQ=255)                                                                                                                cAAATCAGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAgg                       >  2:476852/1‑132 (MQ=255)                                                                                                                      aGTGCCTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGttt         >  1:50700/1‑140 (MQ=255)                                                                                                                           cTTTCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTcg    <  1:225641/140‑1 (MQ=255)                                                                                                                             ttCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTcgcg  >  1:464924/1‑140 (MQ=255)                                                                                                                             ttCGTGGACAGCTAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATACAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTcgcg  >  1:574286/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                CCCGGCACGGTTCAGGCGCGAGATAACAACGTCATAGCGTTCTTTAAAGTGCAGGTCGCGAATACGTTTTCCGAGCACGTTTTCATTGGTGACCACCACACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCG  >  NC_000913/3865296‑3865559

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