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%	152.4 / NA	48	*164* (TAG→TAA)	ptsN	sugar‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/22);  total (26/22)

ACCTGCATACTCTGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAGTTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGTTTTTACTGCG  >  NC_000913/3346936‑3347199                                                                                                                                      |                                                                                                                                   aCCTGCATACTCTGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAattatt                                                                                                                               <  2:1520405/139‑1 (MQ=255) aCCTGCATACTCTGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAattatt                                                                                                                               <  2:1294726/139‑1 (MQ=255)   cTGCATACTCTGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCg                                                                                                                             >  2:696097/1‑139 (MQ=255)          ctctGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                      <  2:175404/139‑1 (MQ=255)          ctctGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                      <  2:214247/139‑1 (MQ=255)          ctctGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                      <  1:696097/139‑1 (MQ=255)             tGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGtctc                                                                                                                   >  2:1243672/1‑139 (MQ=255)             tGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGtct                                                                                                                    >  1:674104/1‑138 (MQ=255)             tGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGtct                                                                                                                    <  2:674104/138‑1 (MQ=255)               tCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAattatt                                                                                                                               <  2:1571943/125‑1 (MQ=255)               tCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAattatt                                                                                                                               >  1:1571943/1‑125 (MQ=255)               tCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTAt                                                                                                                                >  2:649129/1‑124 (MQ=255)               tCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTAt                                                                                                                                <  1:649129/124‑1 (MQ=255)                cGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCttt                                                                                                                >  1:715643/1‑139 (MQ=255)                        gCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGtt                                                                                                        >  1:326613/1‑139 (MQ=255)                         cGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTtg                                                                                                       <  1:1243672/139‑1 (MQ=255)                           aaaCGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTtgtg                                                                                                     <  2:1590835/139‑1 (MQ=255)                                   ggCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGaaa                                                                                             <  2:715643/139‑1 (MQ=255)                                          aaaaCCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGCATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                      >  1:947704/1‑139 (MQ=255)                                                  cTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACa                                                                                           >  2:702596/1‑126 (MQ=255)                                                  cTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACa                                                                                           <  1:702596/126‑1 (MQ=255)                                                                      cccAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACgttc                                                          >  2:370718/1‑139 (MQ=255)                                                                      cccAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACgtt                                                           >  1:586313/1‑138 (MQ=255)                                                                      cccAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACgtt                                                           <  2:586313/138‑1 (MQ=255)                                                                                  aaGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTaa                                              <  2:947704/139‑1 (MQ=255)                                                                                       cTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                      >  2:1117591/1‑94 (MQ=255)                                                                                       cTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                      <  1:1117591/94‑1 (MQ=255)                                                                                          tATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAgg                                                                                                            >  2:1001836/1‑68 (MQ=255)                                                                                          tATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGa                                                                                                            <  1:1001836/69‑1 (MQ=255)                                                                                            tCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGcc                                    >  2:1272969/1‑139 (MQ=255)                                                                                               aaTCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                      <  2:244978/86‑1 (MQ=255)                                                                                               aaTCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                      >  1:244978/1‑86 (MQ=255)                                                                                                 tCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTa                                                                                                                          >  2:735705/1‑48 (MQ=255)                                                                                                 tCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTa                                                                                                                          <  1:735705/48‑1 (MQ=255)                                                                                                    ttACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGaaa                                                                                             <  2:1294938/74‑1 (MQ=255)                                                                                                    ttACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGaaa                                                                                             >  1:1294938/1‑74 (MQ=255)                                                                                                    ttACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCtgcgtgc                            >  1:1491536/1‑139 (MQ=255)                                                                                                          aTACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                      <  2:186690/75‑1 (MQ=255)                                                                                                          aTACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                      >  1:186690/1‑75 (MQ=255)                                                                                                          aTACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGa                      >  2:579491/1‑139 (MQ=255)                                                                                                          aTACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATG‑TCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGc                                     <  1:484584/123‑1 (MQ=255)                                                                                                          aTACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATG‑TCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGc                                     >  2:484584/1‑123 (MQ=255)                                                                                                                aaGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGatat                >  2:784177/1‑139 (MQ=255)                                                                                                                      cTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACa                                                                                           <  2:317396/58‑1 (MQ=255)                                                                                                                      cTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACa                                                                                           >  1:317396/1‑58 (MQ=255)                                                                                                                      cTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGttt          >  2:1274489/1‑139 (MQ=255)                                                                                                                        ccGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACgttcagg                                                       <  1:193960/92‑1 (MQ=255)                                                                                                                        ccGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACgttcagg                                                       >  2:193960/1‑92 (MQ=255)                                                                                                                        ccGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCtgcg                               <  1:997104/116‑1 (MQ=255)                                                                                                                        ccGGATGAAGCGTAATTATTCGGTAATGCCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCtgcg                               >  2:997104/1‑116 (MQ=255)                                                                                                                            aTGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGTTTTTACTg    <  1:1274489/139‑1 (MQ=255)                                                                                                                              gAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGTGTTTACTGCg  >  2:1031880/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                   ACCTGCATACTCTGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAGTTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGTTTTTACTGCG  >  NC_000913/3346936‑3347199

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