Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,373,078	T→C	T167A (ACC→GCC)	nanA ←	N‑acetylneuraminate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,373,078	0	T	C	100.0%	159.3 / NA	43	T167A (ACC→GCC)	nanA	N‑acetylneuraminate lyase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/15);  total (28/15)

CGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCA  >  NC_000913/3372930‑3373222                                                                                                                                                     |                                                                                                                                                 cgcCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCt                                                                                                                                                >  2:428917/1‑151 (MQ=255) cgcCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCt                                                                                                                                                <  1:730348/151‑1 (MQ=255)        ccATGATGTTGTAGGTACTGCCGATCCCCCCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAg                                                                                                                                         <  1:298443/151‑1 (MQ=255)            gatgTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAgcgc                                                                                                                                      >  2:216330/1‑150 (MQ=255)              tgttgtAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCAcc                                                                                                                                   >  1:296004/1‑151 (MQ=255)                   tAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTTATAGATATCGCCAGAGGCCTGTTTCAGCGCACCTACg                                                                                                                               <  1:105117/150‑1 (MQ=255)                     ggTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGcc                                                                                                                             >  1:374241/1‑150 (MQ=255)                         cTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGc                                                                                                                         >  1:442867/1‑150 (MQ=255)                           gCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAAt                                                                                                                      <  2:503233/151‑1 (MQ=255)                            ccGATACCCCCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATg                                                                                                                     <  2:95595/151‑1 (MQ=255)                                        tCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCAAGGCAATGTAACAAgtgtgt                                                                                                         >  2:377455/1‑151 (MQ=255)                                         cAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAgtgtgt                                                                                                         <  2:516240/150‑1 (MQ=255)                                              cccgccAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTtgatct                                                                                                   <  1:216338/151‑1 (MQ=255)                                               ccgccAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTtgatctg                                                                                                  <  2:747937/151‑1 (MQ=255)                                               ccgccAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTtgatctg                                                                                                  >  1:135935/1‑151 (MQ=255)                                                           cAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTc                                                                                        >  2:108001/1‑149 (MQ=255)                                                                 cGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTAc                                                                                >  2:704269/1‑151 (MQ=255)                                                                  gAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTAcc                                                                               >  1:697979/1‑151 (MQ=255)                                                                    aGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAg                                                                                      <  1:84342/142‑1 (MQ=255)                                                                    aGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAg                                                                                      >  2:84341/1‑142 (MQ=255)                                                                     gATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTAcccc                                                                             >  2:322339/1‑150 (MQ=255)                                                                           gtcgtAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAgg                                                                      >  1:314429/1‑151 (MQ=255)                                                                                aaCCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTgg                                                                 >  1:63317/1‑151 (MQ=255)                                                                                aaCCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTgg                                                                 >  2:275079/1‑151 (MQ=255)                                                                                   cGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAAt                                                              >  2:623817/1‑151 (MQ=255)                                                                                     ttATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAAtgt                                                            >  2:163813/1‑151 (MQ=255)                                                                                                  gATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCAt                                               <  2:535407/151‑1 (MQ=255)                                                                                                      aGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGc                                           >  1:517959/1‑151 (MQ=255)                                                                                                        gATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCa                                          >  1:688619/1‑150 (MQ=255)                                                                                                        gATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCa                                          >  1:658467/1‑150 (MQ=255)                                                                                                           gTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAAcc                                      >  2:47399/1‑151 (MQ=255)                                                                                                                    ggATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATccgccg                              >  1:115734/1‑150 (MQ=255)                                                                                                                          gCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGaa                            <  1:612642/146‑1 (MQ=255)                                                                                                                             ccATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCAAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaataat                    >  2:576978/1‑151 (MQ=255)                                                                                                                              cATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAAtt                   >  1:668586/1‑151 (MQ=255)                                                                                                                                tCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGc                 >  1:745729/1‑151 (MQ=255)                                                                                                                                   gATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGccc               >  1:497664/1‑150 (MQ=255)                                                                                                                                   gATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCg              <  1:47400/151‑1 (MQ=255)                                                                                                                                    aTAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCg              >  1:429851/1‑150 (MQ=255)                                                                                                                                         gaTCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTg        <  2:374225/151‑1 (MQ=255)                                                                                                                                             gCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCg    <  1:537513/151‑1 (MQ=255)                                                                                                                                               cAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCa  <  2:457016/151‑1 (MQ=255)                                                                                                                                                agagGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGc   >  2:233504/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCA  >  NC_000913/3372930‑3373222

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