Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,422,588	A→C	K51N (AAA→AAC)	yhdZ →	putative amino acid ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,422,588	0	A	C	100.0%	202.8 / NA	54	K51N (AAA→AAC)	yhdZ	putative amino acid ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (38/16);  total (38/16)

GCCAAATTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGAAGTGGGAATGGTCTTTCAGCATTTCAATCT  >  NC_000913/3422440‑3422731                                                                                                                                                     |                                                                                                                                                gCCAAATTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTc                                                                                                                                               >  2:321989/1‑151 (MQ=255)     aaTTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGa                                                                                                                                             >  2:772937/1‑149 (MQ=255)     aaTTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACa                                                                                                                                           <  1:410598/151‑1 (MQ=255)        tttACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAAcc                                                                                                                                        >  2:913617/1‑151 (MQ=255)        tttACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAAcc                                                                                                                                        >  2:154681/1‑151 (MQ=255)          tACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCAt                                                                                                                                      >  1:227597/1‑151 (MQ=255)               cAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTc                                                                                                                                    <  2:710725/148‑1 (MQ=255)               cAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGtt                                                                                                                                 <  1:572796/151‑1 (MQ=255)                   cTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTAt                                                                                                                             >  2:758888/1‑151 (MQ=255)                        aaCGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAatc                                                                                                                        <  1:869844/151‑1 (MQ=255)                        aaCGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAatc                                                                                                                        >  2:233225/1‑151 (MQ=255)                         aCGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAatc                                                                                                                        >  1:579186/1‑150 (MQ=255)                            cgATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCt                                                                                                                    >  1:518124/1‑151 (MQ=255)                             gatgatTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTg                                                                                                                   >  1:518538/1‑151 (MQ=255)                                    aCGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAAc                                                                                                            <  1:206528/151‑1 (MQ=255)                                      gCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACAt                                                                                                          >  2:128459/1‑151 (MQ=255)                                      gCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACAt                                                                                                          >  1:829191/1‑151 (MQ=255)                                       cTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATc                                                                                                         >  2:726510/1‑151 (MQ=255)                                        tGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATc                                                                                                         >  1:800502/1‑150 (MQ=255)                                            aaaCGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAAc                                                                                                      >  1:127232/1‑149 (MQ=255)                                               cGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAgg                                                                                                   >  2:835869/1‑149 (MQ=255)                                                gTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGAc                                                                                                <  2:584143/151‑1 (MQ=255)                                                       aaTGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGt                                                                                         >  2:481455/1‑151 (MQ=255)                                                        aTGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTg                                                                                        >  2:364968/1‑151 (MQ=255)                                                               ggACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATg                                                                                 >  2:346718/1‑151 (MQ=255)                                                                gACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGAt                                                                                  <  2:178974/149‑1 (MQ=255)                                                                 aCAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGc                                                                               >  2:459018/1‑151 (MQ=255)                                                                      tCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCg                                                                           <  2:620077/150‑1 (MQ=255)                                                                         aTGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACt                                                                       <  1:484282/151‑1 (MQ=255)                                                                               tGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATg                                                                  >  2:634495/1‑150 (MQ=255)                                                                                        taAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCg                                                        >  1:899657/1‑151 (MQ=255)                                                                                         aaaTTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGc                                                       >  2:595093/1‑151 (MQ=255)                                                                                           aTTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCaa                                                     <  1:529360/151‑1 (MQ=255)                                                                                               aaCCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATAtt                                                 >  2:525443/1‑151 (MQ=255)                                                                                                aCCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTg                                                >  2:524180/1‑151 (MQ=255)                                                                                                 ccGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTg                                                >  1:925439/1‑150 (MQ=255)                                                                                                 ccGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTg                                                >  2:256410/1‑150 (MQ=255)                                                                                                  cGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAg                                              <  2:708699/151‑1 (MQ=255)                                                                                                   gTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAgc                                             >  1:152338/1‑151 (MQ=255)                                                                                                      cAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGcgcg                                          >  2:686950/1‑151 (MQ=255)                                                                                                       aaCCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGt                                         >  2:646646/1‑151 (MQ=255)                                                                                                        aCCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTc                                        >  1:108175/1‑151 (MQ=255)                                                                                                            ggAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTcaggc                                    >  2:702377/1‑151 (MQ=255)                                                                                                             gagaACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTcaggca                                   <  1:322927/151‑1 (MQ=255)                                                                                                              agaACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTcaggcag                                  <  1:495055/151‑1 (MQ=255)                                                                                                              agaACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTcaggcag                                  <  1:76365/151‑1 (MQ=255)                                                                                                                aaCGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTcaggcagg                                 >  2:321200/1‑150 (MQ=255)                                                                                                                aaCGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTcaggcagg                                 >  2:600858/1‑150 (MQ=255)                                                                                                                  cGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGa                                >  1:43657/1‑149 (MQ=255)                                                                                                                           cTGTGTGGCCCTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGAAGTGGGAATgg                     >  2:435835/1‑151 (MQ=255)                                                                                                                                    ccTTCAGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGAAGTGGGAATGGTCTTTCa              >  2:905633/1‑149 (MQ=255)                                                                                                                                         aGGTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGAAGTGGGAATGGTCTTTCAGCAttt        >  2:89126/1‑150 (MQ=255)                                                                                                                                          ggTTCCGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGAAGTGGGAATGGTCTTTCAGCATTTCa      <  2:857077/151‑1 (MQ=255)                                                                                                                                              ccGGTAACTCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGAAGTGGGAATGGTCTTTCAGCATTTCAAtct  <  1:726528/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                GCCAAATTTTACTGCAACCTGCTAACGCGATGATTACGCTGGAAAACGTCAATAAATGGTATGGACAATTCCATGTTTTGAAAAATATAAATTTAACCGTGCAACCGGGAGAACGGATCGTTCTGTGTGGCCCTTCAGGTTCCGGTAAATCGACAACCATTCGTTGTATTAATCATCTGGAAGAACATCAACAGGGACGGATCGTGGTAGATGGCATCGAACTTAATGAAGATATCCGCAATATTGAGCGCGTCAGGCAGGAAGTGGGAATGGTCTTTCAGCATTTCAATCT  >  NC_000913/3422440‑3422731

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