Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,372,254	G→A	G107G (GGC→GGT)	nanT ←	sialic acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,372,254	0	G	A	100.0%	114.4 / NA	37	G107G (GGC→GGT)	nanT	sialic acid transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/17);  total (20/17)

CCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCA  >  NC_000913/3372120‑3372389                                                                                                                                       |                                                                                                                                        ccTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTg                                                                                                                                     >  1:464285/1‑139 (MQ=255)       aTCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCaa                                                                                                                               <  1:677249/139‑1 (MQ=255)            aaaaCCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCg                                                                                                                          <  1:855849/139‑1 (MQ=255)                    tGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCa                                                                                                                        <  2:906773/133‑1 (MQ=255)                    tGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCa                                                                                                                        >  1:906773/1‑133 (MQ=255)                          tGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCaa                                                                                                                               >  1:1180437/1‑120 (MQ=255)                          tGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCaa                                                                                                                               <  2:1180437/120‑1 (MQ=255)                                   gATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCgagaa                                                                                                   <  2:464285/139‑1 (MQ=255)                                                                gTGGCACTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCg                                                                                                                          <  2:707574/87‑1 (MQ=255)                                                                gTGGCACTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCg                                                                                                                          >  1:707574/1‑87 (MQ=255)                                                                    cgcTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCgagaa                                                                                                   <  1:226858/106‑1 (MQ=255)                                                                    cgcTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCgagaa                                                                                                   >  2:226858/1‑106 (MQ=255)                                                                    cgcTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGacgacg                                                                  >  2:352996/1‑139 (MQ=255)                                                                         gAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGt                                                             <  2:962267/139‑1 (MQ=255)                                                                          aaCCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTa                                                            >  2:443872/1‑139 (MQ=255)                                                                           aCCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAg                                                           <  2:1093520/139‑1 (MQ=255)                                                                                aTTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTc                                                      >  2:115163/1‑139 (MQ=255)                                                                                    acgcGCCATCCCCATGCCGATGACCAGACGAGCGATAAACCTGGTGACGTAACCCCGCGCAACGCCGCAGGCCCGAGTACCGGACGAACAGAGAACGATGCAGGCGACCATCCCCGGACGACGCCGCTCGCGGtgcccc                                                  >  2:914811/4‑134 (MQ=255)                                                                                      ccGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCAt                                                >  2:595407/1‑139 (MQ=255)                                                                                         ccATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCga                                                                                                      <  1:1192398/82‑1 (MQ=255)                                                                                         ccATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCga                                                                                                      >  2:1192398/1‑82 (MQ=255)                                                                                                  tGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAg                                              >  1:1026885/1‑129 (MQ=255)                                                                                                  tGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCACAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCAt                                    >  2:478594/1‑139 (MQ=255)                                                                                                  tGCCGATGACCAGACGAGCGATAAACATGGAGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAg                                              <  2:1026885/129‑1 (MQ=255)                                                                                                       aTGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGTCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGgc                               >  1:143572/1‑139 (MQ=255)                                                                                                       aTGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGcc                                                                                                        <  1:920792/66‑1 (MQ=255)                                                                                                       aTGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGcc                                                                                                        >  2:920792/1‑66 (MQ=255)                                                                                                       aTGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGgc                               >  2:370300/1‑139 (MQ=255)                                                                                                        tGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGgcc                              >  1:950762/1‑139 (MQ=255)                                                                                                        tGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGgcc                              >  1:102598/1‑139 (MQ=255)                                                                                                              gACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGaa                        <  1:115163/139‑1 (MQ=255)                                                                                                                    cGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCg                  >  2:485213/1‑139 (MQ=255)                                                                                                                             tGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATaaa         <  1:352996/139‑1 (MQ=255)                                                                                                                                   tGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGc   >  2:441370/1‑139 (MQ=255)                                                                                                                                   tGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGc   <  2:950762/139‑1 (MQ=255)                                                                                                                                   tGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGc   <  2:102598/139‑1 (MQ=255)                                                                                                                                    gTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCa  <  1:370300/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                        CCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCA  >  NC_000913/3372120‑3372389

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