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%	167.0 / NA	52	*239* (TAG→TAA)	nanC	N‑acetylnuraminic acid outer membrane channel protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/29);  total (23/29)

GTCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTACTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCA  >  NC_000913/4538649‑4538903                                                                                                                                         |                                                                                                                       gtCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTctatta                                                                                                                      <  2:411047/139‑1 (MQ=255)      tATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTctatta                                                                                                                      <  1:939199/134‑1 (MQ=255)              gCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGac                                                                                                         >  1:71927/1‑139 (MQ=255)               cTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGaca                                                                                                        >  2:470578/1‑139 (MQ=255)                  cGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGacacac                                                                                                     >  2:553484/1‑139 (MQ=255)                   ggTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACAcc                                                                                                    <  1:470578/139‑1 (MQ=255)                    gTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGTAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCa                                                                                                   >  1:910284/1‑139 (MQ=255)                    gTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCa                                                                                                   >  1:125150/1‑139 (MQ=255)                          cTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGa                                                                                             <  1:986253/139‑1 (MQ=255)                                     cacaGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGa                                                                                                          >  1:1246387/1‑115 (MQ=255)                                     cacaGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGa                                                                                                          <  2:1246387/115‑1 (MQ=255)                                     cacaGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTc                                                                                  >  2:715554/1‑139 (MQ=255)                                                    tAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATaaa                                                                           >  2:1025505/1‑131 (MQ=255)                                                    tAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATaaa                                                                           <  1:1025505/131‑1 (MQ=255)                                                     aaTACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                  <  2:142824/139‑1 (MQ=255)                                                         cAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGt                                                              <  1:243921/139‑1 (MQ=255)                                                            aCGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGt                                                           <  2:910284/139‑1 (MQ=255)                                                              gcgTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGt                                                                                                                  >  2:493946/1‑82 (MQ=255)                                                              gcgTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGt                                                                                                                  <  1:493946/82‑1 (MQ=255)                                                              gcgTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTaa                                                         >  1:962096/1‑139 (MQ=255)                                                                                gCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAgg                                       <  2:540118/139‑1 (MQ=255)                                                                                gCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCCATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAgg                                       <  1:715554/139‑1 (MQ=255)                                                                                  catcatGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTa                                     <  2:125150/139‑1 (MQ=255)                                                                                       tGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCat                                <  1:1197655/139‑1 (MQ=255)                                                                                       tGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAAGGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCat                                <  2:962096/139‑1 (MQ=255)                                                                                         aTAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCatat                              <  1:473566/139‑1 (MQ=255)                                                                                           aGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTc                            >  2:967199/1‑139 (MQ=255)                                                                                                gcgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatgt                       >  2:1231066/1‑139 (MQ=255)                                                                                                 cgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatgta                      <  2:520628/139‑1 (MQ=255)                                                                                                    cGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATgg                   >  1:179558/1‑139 (MQ=255)                                                                                                      ttATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCg                 <  1:967199/139‑1 (MQ=255)                                                                                                                aTTCATTGTTACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTAGCTCTGCCGTTGTAAACACCCTGACGGTc                                           <  2:1088066/103‑1 (MQ=255)                                                                                                                aTTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTc                                           >  1:1088066/1‑103 (MQ=255)                                                                                                                aTTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCggg       <  2:71927/139‑1 (MQ=255)                                                                                                                aTTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCggg       <  1:1231066/139‑1 (MQ=255)                                                                                                                 ttCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                  >  2:239046/1‑79 (MQ=255)                                                                                                                 ttCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                  <  1:239046/79‑1 (MQ=255)                                                                                                                 ttCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTa                                                          >  1:575171/1‑87 (MQ=255)                                                                                                                 ttCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTa                                                          <  2:575171/87‑1 (MQ=255)                                                                                                                    aTTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatgtatg                    >  1:1059668/1‑122 (MQ=255)                                                                                                                    aTTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatgtatg                    <  2:1059668/122‑1 (MQ=255)                                                                                                                     ttGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCa  <  1:1043032/139‑1 (MQ=255)                                                                                                                      tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGTGATAACTGTTTTCCGATAAATTATCTCTGcc                                                                >  2:951253/1‑76 (MQ=255)                                                                                                                      tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGa                                                                               <  2:283378/61‑1 (MQ=255)                                                                                                                      tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGa                                                                               >  1:283378/1‑61 (MQ=255)                                                                                                                      tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGcc                                                                <  1:951253/76‑1 (MQ=255)                                                                                                                              tctcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                  >  2:438873/1‑66 (MQ=255)                                                                                                                              tctcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                  <  1:438873/66‑1 (MQ=255)                                                                                                                              tctcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGt                                   <  1:207817/97‑1 (MQ=255)                                                                                                                              tctcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGt                                   >  2:207817/1‑97 (MQ=255)                                                                                                                              tctcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCggg       <  1:371279/125‑1 (MQ=255)                                                                                                                              tctcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCggg       >  2:371279/1‑125 (MQ=255)                                                                                                                                     tCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATgg                   <  2:932165/106‑1 (MQ=255)                                                                                                                                     tCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATgg                   >  1:932165/1‑106 (MQ=255)                                                                                                                                         |                                                                                                                       GTCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTACTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCA  >  NC_000913/4538649‑4538903

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: ‑