Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,347,069	G→A	*164* (TAG→TAA)	ptsN →	sugar‑specific enzyme IIA component of PTS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,347,069	0	G	A	100.0%	202.6 / NA	62	*164* (TAG→TAA)	ptsN	sugar‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/32);  total (30/32)

GCACCTGCATACTCTGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAGTTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGTTTTTACTG  >  NC_000913/3346934‑3347197                                                                                                                                        |                                                                                                                                 gCACCTGCATACTCTGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTa                                                                                                                               <  2:1220895/139‑1 (MQ=255)               tGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGtctc                                                                                                                 >  2:523889/1‑139 (MQ=255)                gTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTctct                                                                                                                <  1:199696/139‑1 (MQ=255)                 tCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCtt                                                                                                               >  1:78648/1‑139 (MQ=255)                  cGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCttt                                                                                                              >  1:687175/1‑139 (MQ=255)                     tggtggCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAg                                                                                                           >  2:702620/1‑139 (MQ=255)                       gtggCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGAc                                                                                                         <  2:78648/139‑1 (MQ=255)                       gtggCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGAc                                                                                                         >  1:97346/1‑139 (MQ=255)                                cGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTgagg                                                                                                <  2:147792/139‑1 (MQ=255)                                cGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATATCTCTTTTAGACGTTGTgagg                                                                                                <  2:1052626/139‑1 (MQ=255)                                    tGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGaa                                                                                            <  2:1134915/139‑1 (MQ=255)                                     ggCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGaaa                                                                                           <  1:523889/139‑1 (MQ=255)                                     ggCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGaaa                                                                                           <  1:702620/139‑1 (MQ=255)                                            aaaaCCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                    <  2:687175/139‑1 (MQ=255)                                                  aTCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTAc                                                                              <  2:97346/139‑1 (MQ=255)                                                          ccgTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGt                                                                                       <  1:413994/122‑1 (MQ=255)                                                          ccgTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGt                                                                                       >  2:413994/1‑122 (MQ=255)                                                                               cGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACgttcaggttca                                                 <  2:1257073/139‑1 (MQ=255)                                                                                  tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAg                                                                                                           <  1:393008/78‑1 (MQ=255)                                                                                  tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAg                                                                                                           >  2:393008/1‑78 (MQ=255)                                                                                  tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                    <  1:457354/101‑1 (MQ=255)                                                                                  tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                    >  2:457354/1‑101 (MQ=255)                                                                                  tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTtcaggttcaggt                                              >  2:1282271/1‑139 (MQ=255)                                                                                  tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTtcaggttcaggt                                              >  2:249713/1‑139 (MQ=255)                                                                                            tATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCg                                    >  2:469625/1‑139 (MQ=255)                                                                                              tCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAattatt                                                                                                                             >  2:885981/1‑48 (MQ=255)                                                                                              tCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAattatt                                                                                                                             <  1:885981/48‑1 (MQ=255)                                                                                                    cATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTa                                                                                      <  2:1096765/81‑1 (MQ=255)                                                                                                    cATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTa                                                                                      >  1:1096765/1‑81 (MQ=255)                                                                                                     aTTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCtgcgtg                           <  1:1282271/139‑1 (MQ=255)                                                                                                      ttACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                    >  1:148057/1‑81 (MQ=255)                                                                                                      ttACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                    <  2:148057/81‑1 (MQ=255)                                                                                                      ttACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACgtt                                                         <  2:274985/108‑1 (MQ=255)                                                                                                      ttACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACgtt                                                         >  1:274985/1‑108 (MQ=255)                                                                                                            aTACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGa                                                             <  1:287233/98‑1 (MQ=255)                                                                                                            aTACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGa                                                             >  2:287233/1‑98 (MQ=255)                                                                                                            aTACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTtcaggttcaggt                                              >  2:208535/1‑113 (MQ=255)                                                                                                            aTACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTtcaggttcaggt                                              <  1:208535/113‑1 (MQ=255)                                                                                                            aTACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGt                                      >  2:674390/1‑121 (MQ=255)                                                                                                            aTACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGt                                      <  1:674390/121‑1 (MQ=255)                                                                                                              aCCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTa                                             <  2:238265/112‑1 (MQ=255)                                                                                                              aCCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTa                                             >  1:238265/1‑112 (MQ=255)                                                                                                              aCCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAg                  >  1:1304319/1‑139 (MQ=255)                                                                                                              aCCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAg                  >  1:1374412/1‑139 (MQ=255)                                                                                                              aCCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAg                  >  1:165345/1‑139 (MQ=255)                                                                                                                cGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGat                >  2:207522/1‑139 (MQ=255)                                                                                                                        cTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTa                                                                                      <  1:515595/61‑1 (MQ=255)                                                                                                                        cTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTa                                                                                      >  2:515595/1‑61 (MQ=255)                                                                                                                          ccGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTgagg                                                                                                <  1:685272/49‑1 (MQ=255)                                                                                                                          ccGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTgagg                                                                                                >  2:685272/1‑49 (MQ=255)                                                                                                                          ccGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAAt                                          <  2:1400610/103‑1 (MQ=255)                                                                                                                          ccGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAAt                                          >  1:1400610/1‑103 (MQ=255)                                                                                                                          ccGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGata               >  1:1040787/1‑130 (MQ=255)                                                                                                                          ccGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGata               <  2:1040787/130‑1 (MQ=255)                                                                                                                          ccGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGttttt      <  2:165345/139‑1 (MQ=255)                                                                                                                           cGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGTTTTTa     <  1:249713/139‑1 (MQ=255)                                                                                                                             gATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTgagg                                                                                                <  2:417266/46‑1 (MQ=255)                                                                                                                             gATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTgagg                                                                                                >  1:417266/1‑46 (MQ=255)                                                                                                                             gATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGACGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGTTTTTACt   >  2:1010190/1‑139 (MQ=255)                                                                                                                              aTGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGTTTTTACTg  <  2:1304319/139‑1 (MQ=255)                                                                                                                                  aGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTg                      >  2:56595/1‑115 (MQ=255)                                                                                                                                  aGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTg                      <  1:56595/115‑1 (MQ=255)                                                                                                                                   gCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCg                                                               >  1:622268/1‑73 (MQ=255)                                                                                                                                   gCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCg                                                               <  2:622268/73‑1 (MQ=255)                                                                                                                                        |                                                                                                                                 GCACCTGCATACTCTGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAGTTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGTTTTTACTG  >  NC_000913/3346934‑3347197

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