Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,538,785	C→T	*239* (TAG→TAA)	nanC ←	N‑acetylnuraminic acid outer membrane channel protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,785	0	C	T	100.0%	218.2 / NA	68	*239* (TAG→TAA)	nanC	N‑acetylnuraminic acid outer membrane channel protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/35);  total (33/35)

GTCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTACTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAA  >  NC_000913/4538649‑4538918                                                                                                                                         |                                                                                                                                      gtCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTctatta                                                                                                                                     <  1:504229/139‑1 (MQ=255)  tCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCtattac                                                                                                                                    >  2:874147/1‑139 (MQ=255)     ttATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGt                                                                                                                                 <  2:887726/139‑1 (MQ=255)            ttGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATg                                                                                                                          >  1:174538/1‑139 (MQ=255)            ttGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATg                                                                                                                          >  2:963933/1‑139 (MQ=255)            ttGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATg                                                                                                                          >  2:128558/1‑139 (MQ=255)                   ggTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACAcc                                                                                                                   <  2:988324/139‑1 (MQ=255)                                            gTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATATCTGTTTTCCGATaaa                                                                                          <  2:154720/139‑1 (MQ=255)                                               tCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTa                                                                                       >  2:772005/1‑139 (MQ=255)                                                    tAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTG‑TTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                                 <  2:34725/139‑1 (MQ=255)                                                     aaTACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                                 >  1:994321/1‑139 (MQ=255)                                                       tACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGcc                                                                               >  2:992046/1‑139 (MQ=255)                                                           gACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAAc                                                                                                        <  1:185951/110‑1 (MQ=255)                                                           gACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAAc                                                                                                        >  2:185951/1‑110 (MQ=255)                                                            aCGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGt                                                                          >  1:339702/1‑139 (MQ=255)                                                            aCGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGt                                                                          >  1:180667/1‑139 (MQ=255)                                                                 tttGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAaca                                                                     >  2:570667/1‑139 (MQ=255)                                                                     cGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGCTGTAAACACCCt                                                                 <  1:992046/139‑1 (MQ=255)                                                                      ggCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTg                                                                <  2:180667/139‑1 (MQ=255)                                                                                gCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAgg                                                      <  1:482573/139‑1 (MQ=255)                                                                                gCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAgg                                                      <  1:963933/139‑1 (MQ=255)                                                                                    tcatGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGt                                                  <  2:309809/139‑1 (MQ=255)                                                                                    tcatGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGt                                                  <  2:857489/139‑1 (MQ=255)                                                                                      atGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCa                                                <  2:339702/139‑1 (MQ=255)                                                                                           aGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTc                                           <  2:994321/139‑1 (MQ=255)                                                                                             cAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGAt                                                                                             >  1:1002507/1‑87 (MQ=255)                                                                                             cAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGAt                                                                                             <  2:1002507/87‑1 (MQ=255)                                                                                               aGCGCCGTTATTGTTTTATTCATTGTGACTTTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGcc                                                                               <  1:385552/99‑1 (MQ=255)                                                                                               aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatg                                       >  1:34907/1‑139 (MQ=255)                                                                                               aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatg                                       >  2:564773/1‑139 (MQ=255)                                                                                               aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAAATGTTTTCCGATAAATTATCTCTGcc                                                                               >  2:385552/1‑99 (MQ=255)                                                                                                gcgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGcc                                                                               <  1:444312/98‑1 (MQ=255)                                                                                                gcgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGcc                                                                               >  2:444312/1‑98 (MQ=255)                                                                                                 cgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAgg                                                      >  1:857212/1‑122 (MQ=255)                                                                                                 cgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAgg                                                      <  2:857212/122‑1 (MQ=255)                                                                                                 cgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatgta                                     <  1:564773/139‑1 (MQ=255)                                                                                                  gcCGTCATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatgtat                                    >  1:62915/1‑139 (MQ=255)                                                                                                   ccGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatgtatg                                   >  1:379218/1‑139 (MQ=255)                                                                                                    cGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATgg                                  >  1:549562/1‑139 (MQ=255)                                                                                                      ttATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGaca                                                                                                                       >  2:928999/1‑52 (MQ=255)                                                                                                      ttATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGaca                                                                                                                       <  1:928999/52‑1 (MQ=255)                                                                                                                 ttCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGc                                 >  2:366202/1‑127 (MQ=255)                                                                                                                 ttCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGc                                 <  1:366202/127‑1 (MQ=255)                                                                                                                 ttCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGc                     >  2:225895/1‑139 (MQ=255)                                                                                                                      tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATaa                                                                                                         <  2:605248/50‑1 (MQ=255)                                                                                                                      tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATaa                                                                                                         >  1:605248/1‑50 (MQ=255)                                                                                                                       gtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCata               <  2:174538/139‑1 (MQ=255)                                                                                                                       gtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCata               <  2:150112/139‑1 (MQ=255)                                                                                                                       gtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCata               <  1:570667/139‑1 (MQ=255)                                                                                                                         gACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATg             <  2:379218/139‑1 (MQ=255)                                                                                                                          aCTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAaca                                                                     >  1:75621/1‑82 (MQ=255)                                                                                                                          aCTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAaca                                                                     <  2:75621/82‑1 (MQ=255)                                                                                                                          aCTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATgg            <  2:62915/139‑1 (MQ=255)                                                                                                                           cTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAgg                                                      >  2:696005/1‑96 (MQ=255)                                                                                                                           cTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAgg                                                      <  1:696005/96‑1 (MQ=255)                                                                                                                              tctcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCg                    >  2:106171/1‑127 (MQ=255)                                                                                                                              tctcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCg                    <  1:106171/127‑1 (MQ=255)                                                                                                                                tcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTcc                                                                                                <  2:745906/49‑1 (MQ=255)                                                                                                                                tcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTcc                                                                                                >  1:745906/1‑49 (MQ=255)                                                                                                                                tcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCggg                      <  2:279313/123‑1 (MQ=255)                                                                                                                                tcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCggg                      >  1:279313/1‑123 (MQ=255)                                                                                                                                    gTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAgg                                                      >  2:852734/1‑87 (MQ=255)                                                                                                                                    gTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAgg                                                      <  1:852734/87‑1 (MQ=255)                                                                                                                                    gTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGt                                                  >  1:645734/1‑91 (MQ=255)                                                                                                                                    gTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGt                                                  <  2:645734/91‑1 (MQ=255)                                                                                                                                    gTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTc                                                 <  1:833497/92‑1 (MQ=255)                                                                                                                                    gTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTc                                                 >  2:833497/1‑92 (MQ=255)                                                                                                                                    gTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGaa  >  2:878153/1‑139 (MQ=255)                                                                                                                                     tCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCg                                                                              >  2:335606/1‑62 (MQ=255)                                                                                                                                     tCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCg                                                                              <  1:335606/62‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      GTCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTACTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAA  >  NC_000913/4538649‑4538918

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