Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,372,680	A→C	intergenic (‑106/+3)	nanT ← / ← nanA	sialic acid transporter/N‑acetylneuraminate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,372,680	0	A	C	100.0%	168.4 / NA	45	intergenic (‑106/+3)	nanT/nanA	sialic acid transporter/N‑acetylneuraminate lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (23/22);  total (23/22)

GCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGA  >  NC_000913/3372539‑3372822                                                                                                                                              |                                                                                                                                               gCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCcg                                                                                                                                       <  1:36514/151‑1 (MQ=255)  cGATACCACGGGATATTCTGGGTTGTAATACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCcgc                                                                                                                                      >  2:365407/1‑151 (MQ=255)   gATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCcgcg                                                                                                                                     <  1:518226/151‑1 (MQ=255)    aTACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCcgcg                                                                                                                                     >  1:86590/1‑150 (MQ=255)        cACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCtt                                                                                                                                <  2:674563/151‑1 (MQ=255)         aCGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTg                                                                                                                               >  1:713875/1‑151 (MQ=255)         aCGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTg                                                                                                                               >  1:730453/1‑151 (MQ=255)           gggATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCa                                                                                                                             <  2:762858/151‑1 (MQ=255)             gATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATc                                                                                                                           <  1:524945/151‑1 (MQ=255)                     gggTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAActg                                                                                                                      <  2:854632/148‑1 (MQ=255)                      ggTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTgg                                                                                                                  >  2:845839/1‑151 (MQ=255)                       gTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGACGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTggg                                                                                                                 <  2:345031/151‑1 (MQ=255)                        ttGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGc                                                                                                                >  1:425797/1‑151 (MQ=255)                           tagtaCTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGcc                                                                                                               >  1:175290/1‑149 (MQ=255)                                       tgAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTc                                                                                                 >  1:323275/1‑151 (MQ=255)                                        gAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCt                                                                                                <  1:570180/151‑1 (MQ=255)                                              cccGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGc                                                                                             <  2:932745/148‑1 (MQ=255)                                              cccGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGa                                                                                          <  1:441464/151‑1 (MQ=255)                                               ccGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGc                                                                                             >  1:368096/1‑147 (MQ=255)                                                cgcgGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGat                                                                                         >  1:723837/1‑150 (MQ=255)                                                   ggTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATAttt                                                                                     >  1:685897/1‑151 (MQ=255)                                                      gggCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTc                                                                                  <  2:51537/151‑1 (MQ=255)                                                          gATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCt                                                                              <  1:704315/151‑1 (MQ=255)                                                                  cacGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTc                                                                       <  2:735921/150‑1 (MQ=255)                                                                      cTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAAc                                                                  >  1:623365/1‑151 (MQ=255)                                                                       tttGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACg                                                                 >  1:333552/1‑151 (MQ=255)                                                                        ttGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACgg                                                                >  2:112993/1‑151 (MQ=255)                                                                             aTGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTgc                                                           <  1:558809/151‑1 (MQ=255)                                                                             aTGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTg                                                            <  2:139870/150‑1 (MQ=255)                                                                              tGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTgc                                                           >  2:313200/1‑150 (MQ=255)                                                                                 aaaTTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTgcggc                                                        >  2:373709/1‑150 (MQ=255)                                                                                     tGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGTAGATATTTTTCATCTACCGGTCCAAACGGTTTgcggcacagc                                                   >  2:884634/1‑151 (MQ=255)                                                                                      gTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTgcggcacagc                                                   <  1:130435/150‑1 (MQ=255)                                                                                         gggTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCGCCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTgcggcacagcggca                                               <  1:91655/151‑1 (MQ=255)                                                                                           gTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAGCGGTTTgcggcacagcggcaca                                             <  1:316155/151‑1 (MQ=255)                                                                                                     gTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACAt                                   >  2:684985/1‑151 (MQ=255)                                                                                                      tttATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTTCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTCGCGGCACAGCGGCACAGAAACGACATc                                  >  1:228884/1‑151 (MQ=255)                                                                                                         aTTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCat                               >  2:391170/1‑151 (MQ=255)                                                                                                            tCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCatat                             <  1:282227/150‑1 (MQ=255)                                                                                                                   aCCCGGTAGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATggagg                     >  1:629953/1‑151 (MQ=255)                                                                                                                          aGGGGCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGttt              >  1:31021/1‑151 (MQ=255)                                                                                                                            gggCGAGCGAGGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGTTTTc            <  1:894805/151‑1 (MQ=255)                                                                                                                                gagcgagGGGAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGTTTTCAGGc        <  2:737763/151‑1 (MQ=255)                                                                                                                                      ggggAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGa  <  1:281851/151‑1 (MQ=255)                                                                                                                                       gggAAACCACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGa  >  2:129880/1‑150 (MQ=255)                                                                                                                                              |                                                                                                                                               GCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGATGCCCACACGCTTTGGTATGAAAATTGTAGGGTACAGATGCGTTTATTTCCCCTCACCCGGTAGGGGCGAGCGAGGGGAAACAACTCACCCGCGCTCTTGCATCAACTGCTGGGCCAGCGCCTTCAGTTCTGGCAGATATTTTTCATCTACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGA  >  NC_000913/3372539‑3372822

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