Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,942,993	G→A	G111E (GGG→GAG)	znuC →	zinc ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,942,993	0	G	A	100.0%	223.5 / NA	69	G111E (GGG→GAG)	znuC	zinc ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/37);  total (32/37)

CGCATCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGGGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACTCAAGGCGTGGATGTA  >  NC_000913/1942863‑1943120                                                                                                                                   |                                                                                                                                cgcATCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGa                                                                                                                         >  1:1306846/1‑139 (MQ=255)   cATCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGAtt                                                                                                                       >  2:880699/1‑139 (MQ=255)               aCCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATg                                                                                                           <  2:1306846/139‑1 (MQ=255)                   cAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCaaa                                                                                                       <  2:609727/139‑1 (MQ=255)                        gCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGctc                                                                                                  >  2:1593992/1‑139 (MQ=255)                        gCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGctc                                                                                                  >  2:591581/1‑139 (MQ=255)                                     accacGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGaaa                                                                                     <  2:67532/139‑1 (MQ=255)                                      ccacGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAAc                                                                                    >  1:1686399/1‑139 (MQ=255)                                         cGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCa                                                                                 >  1:475586/1‑139 (MQ=255)                                         cGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCa                                                                                 <  2:864190/139‑1 (MQ=255)                                         cGTTGCCACTGACCGTAAACCG‑TTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAg                                                                                >  1:1199746/1‑139 (MQ=255)                                          gTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAg                                                                                <  1:119869/139‑1 (MQ=255)                                          gTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAg                                                                                >  1:39802/1‑139 (MQ=255)                                           ttGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGc                                                                               <  1:1593992/139‑1 (MQ=255)                                           ttGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGc                                                                               <  1:591581/139‑1 (MQ=255)                                              ccACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCgtg                                                                            <  2:1435687/139‑1 (MQ=255)                                                    aCCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCAc                                                                                                                             <  1:1389131/84‑1 (MQ=255)                                                    aCCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCAc                                                                                                                             >  2:1389131/1‑84 (MQ=255)                                                      cGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCaaa                                                                                                       <  1:1452369/104‑1 (MQ=255)                                                      cGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCaaa                                                                                                       >  2:1452369/1‑104 (MQ=255)                                                             ccgtttttACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCGAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAg                                                             <  2:1199746/136‑1 (MQ=255)                                                                tttttACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTaa                                                                                                                     >  1:1206809/1‑80 (MQ=255)                                                                tttttACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTaa                                                                                                                     <  2:1206809/80‑1 (MQ=255)                                                                            cgcCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGcc                                              <  2:475586/139‑1 (MQ=255)                                                                            cgcCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGcc                                              <  2:39802/139‑1 (MQ=255)                                                                             gcCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGccc                                             <  2:1606179/139‑1 (MQ=255)                                                                                 tGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGccc                                             >  1:635589/1‑135 (MQ=255)                                                                                 tGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGccc                                             <  2:635589/135‑1 (MQ=255)                                                                                 tGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCaa                                         >  1:1396989/1‑139 (MQ=255)                                                                                   gTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCAc                                                                                                                             >  1:542415/1‑53 (MQ=255)                                                                                   gTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCAc                                                                                                                             <  2:542415/53‑1 (MQ=255)                                                                                      cacaTAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCCGGCCGAGCACCTGATTAACGCCCCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAattatt                                    >  2:318558/1‑139 (MQ=255)                                                                                         aTAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTaa                                                                                                                     >  1:999664/1‑55 (MQ=255)                                                                                         aTAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTaa                                                                                                                     <  2:999664/55‑1 (MQ=255)                                                                                           aaaGAAGATATTTTTCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCa                                                                                 <  2:431615/89‑1 (MQ=255)                                                                                           aaaGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAACAGCTCTCTGGTGGCGAAACGCa                                                                                 >  1:431615/1‑89 (MQ=255)                                                                                             agaagaTATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAg                                                                                                     >  1:1084341/1‑67 (MQ=255)                                                                                             agaagaTATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAg                                                                                                     <  2:111174/67‑1 (MQ=255)                                                                                             agaagaTATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAg                                                                                                     <  2:1084341/67‑1 (MQ=255)                                                                                             agaagaTATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAg                                                                                                     >  1:111174/1‑67 (MQ=255)                                                                                             agaagaTATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTg                             >  1:1100428/1‑139 (MQ=255)                                                                                                   tatTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGaa                       >  2:767753/1‑139 (MQ=255)                                                                                                        tGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCAc                  <  1:61372/139‑1 (MQ=255)                                                                                                        tGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCAc                  <  2:1100428/139‑1 (MQ=255)                                                                                                         gCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCAc                                                                                                                <  1:1364243/44‑1 (MQ=255)                                                                                                         gCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCAc                                                                                                                >  2:1364243/1‑44 (MQ=255)                                                                                                            tGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACg                                                                                   >  1:1531057/1‑70 (MQ=255)                                                                                                            tGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACg                                                                                   <  2:1531057/70‑1 (MQ=255)                                                                                                            tGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCa                                                                                 >  2:729741/1‑72 (MQ=255)                                                                                                            tGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCa                                                                                 <  1:729741/72‑1 (MQ=255)                                                                                                              cACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGc                                                                                  >  2:103665/1‑69 (MQ=255)                                                                                                              cACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGc                                                                                  <  1:103665/69‑1 (MQ=255)                                                                                                              cACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAg                                                                                >  1:1210162/1‑71 (MQ=255)                                                                                                              cACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAg                                                                                <  2:1210162/71‑1 (MQ=255)                                                                                                               aCTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCa                                                                                 <  1:1186811/69‑1 (MQ=255)                                                                                                               aCTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCa                                                                                 >  2:1186811/1‑69 (MQ=255)                                                                                                                cTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACTCAAGGCg          <  1:880699/139‑1 (MQ=255)                                                                                                                    aaCGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACTCAAGGCGTGGa      <  2:1686399/139‑1 (MQ=255)                                                                                                                    aaCGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACTCAAGGCGTGGa      <  2:1396989/139‑1 (MQ=255)                                                                                                                        tgtCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCa                                                                                 <  1:1659818/60‑1 (MQ=255)                                                                                                                        tgtCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCa                                                                                 >  2:1659818/1‑60 (MQ=255)                                                                                                                        tgtCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACTCAAGGCGTGGATGTa  <  1:987867/139‑1 (MQ=255)                                                                                                                        tgtCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACTCAAGGCGTGGATGTa  >  2:118895/1‑139 (MQ=255)                                                                                                                           ccAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCgtg                                                                            >  1:654737/1‑62 (MQ=255)                                                                                                                           ccAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCgtg                                                                            <  2:654737/62‑1 (MQ=255)                                                                                                                             aGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGa                           >  1:856285/1‑109 (MQ=255)                                                                                                                             aGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGa                           <  2:856285/109‑1 (MQ=255)                                                                                                                               gCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAg                                                                                >  1:46555/1‑54 (MQ=255)                                                                                                                               gCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAg                                                                                <  2:46555/54‑1 (MQ=255)                                                                                                                                   |                                                                                                                                CGCATCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGGGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACTCAAGGCGTGGATGTA  >  NC_000913/1942863‑1943120

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