Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,051,785	T→C	C34C (TGT→TGC)	ymcF →	protein YmcF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,051,785	0	T	C	100.0%	216.2 / NA	60	C34C (TGT→TGC)	ymcF	protein YmcF
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/24);  total (36/24)

CGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTT  >  NC_000913/1051648‑1051917                                                                                                                                          |                                                                                                                                     cgGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCAt                                                                                                                                    >  1:559443/1‑140 (MQ=255)        gAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGt                                                                                                                             <  1:811305/140‑1 (MQ=255)              tgttACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATcg                                                                                                                       >  2:863918/1‑140 (MQ=255)              tgttACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATcg                                                                                                                       >  2:330560/1‑140 (MQ=255)               gttACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATcga                                                                                                                      <  2:869085/140‑1 (MQ=255)                 tACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATcgacg                                                                                                                    <  2:559443/140‑1 (MQ=255)                     ctctAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACgatga                                                                                                                >  1:282731/1‑140 (MQ=255)                      tctAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACgatgat                                                                                                               >  1:74408/1‑140 (MQ=255)                        tAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCa                                                                                                             >  2:52183/1‑140 (MQ=255)                               gCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACgatgat                                                                                                               >  1:728003/1‑131 (MQ=255)                               gCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACgatgat                                                                                                               <  2:728003/131‑1 (MQ=255)                                ccATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGAt                                                                                                     >  2:639534/1‑140 (MQ=255)                                          tCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGt                                                                                                >  2:444511/1‑135 (MQ=255)                                          tCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGt                                                                                                <  1:444511/135‑1 (MQ=255)                                          tCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCgcgc                                                                                           >  2:756624/1‑140 (MQ=255)                                            aaCATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTg                                                                                         <  1:566158/140‑1 (MQ=255)                                              cATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCa                                                                                       >  2:718267/1‑140 (MQ=255)                                              cATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCa                                                                                       >  1:253414/1‑140 (MQ=255)                                              cATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGAAGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCa                                                                                       >  2:666864/1‑140 (MQ=255)                                                tctcCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCaaa                                                                                     >  1:178633/1‑140 (MQ=255)                                                   ccATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATAc                                                                                  >  2:781291/1‑140 (MQ=255)                                                        tCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCaaa                                                                                     <  2:578877/132‑1 (MQ=255)                                                        tCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCaaa                                                                                     >  1:578877/1‑132 (MQ=255)                                                        tCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACt                                                                             >  1:170477/1‑140 (MQ=255)                                                                ccATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATa                                                                                   >  2:141421/1‑126 (MQ=255)                                                                ccATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATa                                                                                   <  1:141421/126‑1 (MQ=255)                                                                 cATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATgcg                                                                    >  1:144149/1‑140 (MQ=255)                                                                         cATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGa                                                            <  1:666864/140‑1 (MQ=255)                                                                              ttCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTca                                                       <  1:929727/140‑1 (MQ=255)                                                                                  caGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACaaa                                                   <  2:144149/140‑1 (MQ=255)                                                                                      aCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTcac                                                      <  1:187523/133‑1 (MQ=255)                                                                                      aCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTcac                                                      >  2:187523/1‑133 (MQ=255)                                                                                        cGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAataata                                             <  2:178633/140‑1 (MQ=255)                                                                                        cGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAataata                                             >  1:366592/1‑140 (MQ=255)                                                                                          cacaTCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATc                                           >  2:532350/1‑140 (MQ=255)                                                                                          cacaTCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATc                                           >  1:917644/1‑140 (MQ=255)                                                                                          cacaTCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATc                                           >  2:479376/1‑140 (MQ=255)                                                                                          cacaTCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATc                                           >  2:776795/1‑140 (MQ=255)                                                                                             aTCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAgg                                        <  2:253414/140‑1 (MQ=255)                                                                                             aTCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAgg                                        <  2:282731/140‑1 (MQ=255)                                                                                               cGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCATAAAATAATATCAGGCt                                      >  2:550749/1‑140 (MQ=255)                                                                                               cGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAgg                                        >  1:144150/1‑138 (MQ=255)                                                                                               cGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAgg                                        <  2:144150/138‑1 (MQ=255)                                                                                               cGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCt                                      >  1:358546/1‑140 (MQ=255)                                                                                                 gCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAgg                                        >  2:332097/1‑136 (MQ=255)                                                                                                 gCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAgg                                        <  1:332097/136‑1 (MQ=255)                                                                                                  cATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTccc                                   >  1:100444/1‑140 (MQ=255)                                                                                                  cATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTccc                                   <  1:718267/140‑1 (MQ=255)                                                                                                       gATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCa                                          <  2:267948/128‑1 (MQ=255)                                                                                                       gATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCa                                          <  2:263453/128‑1 (MQ=255)                                                                                                       gATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCa                                          >  1:267948/1‑128 (MQ=255)                                                                                                       gATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCa                                          >  1:263453/1‑128 (MQ=255)                                                                                                        aTGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTgg                             <  1:479376/140‑1 (MQ=255)                                                                                                         tgtgACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGa                            >  1:564933/1‑140 (MQ=255)                                                                                                           tgACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAgg                                        >  1:790275/1‑126 (MQ=255)                                                                                                           tgACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAgg                                        <  2:790275/126‑1 (MQ=255)                                                                                                            gACGGAGCGAAACCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGcc                         <  1:532350/140‑1 (MQ=255)                                                                                                                     aaaCCCTTTGGGCGCTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTtata                <  1:756624/140‑1 (MQ=255)                                                                                                                                   cTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTAttttt  >  2:207839/1‑140 (MQ=255)                                                                                                                                   cTAAGTGCATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTAttttt  >  2:112690/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                     CGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTT  >  NC_000913/1051648‑1051917

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