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%	149.2 / NA	48	G107G (GGC→GGT)	nanT	sialic acid transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/23);  total (25/23)

CCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAA  >  NC_000913/3372120‑3372382                                                                                                                                       |                                                                                                                                 ccTGAAATCAAAAAACCATTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTg                                                                                                                              >  1:644981/1‑139 (MQ=255)           aaaaaCCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGcc                                                                                                                    >  1:521333/1‑139 (MQ=255)                  aCTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGcc                                                                                                             >  2:572163/1‑139 (MQ=255)                                gCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCga                                                                                               >  2:1332390/1‑139 (MQ=255)                                gCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCga                                                                                               >  2:177378/1‑139 (MQ=255)                                  aGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCgaga                                                                                             <  1:807873/139‑1 (MQ=255)                                      gTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCgagaagag                                                                                         >  1:277977/1‑139 (MQ=255)                                         ttGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAAc                                                                                      <  1:572163/139‑1 (MQ=255)                                          tGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACg                                                                                     <  2:644981/139‑1 (MQ=255)                                                 cTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTgg                                                                              >  2:741202/1‑139 (MQ=255)                                                 cTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTgg                                                                              <  2:277977/139‑1 (MQ=255)                                                             tAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCa                                                                  >  2:983193/1‑139 (MQ=255)                                                                    cgcTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGacgacg                                                           >  1:789561/1‑139 (MQ=255)                                                                          aaCCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCa                                                                                                                 >  2:1380009/1‑79 (MQ=255)                                                                          aaCCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCa                                                                                                                 <  1:1380009/79‑1 (MQ=255)                                                                                   cACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGcc                                                                                                             <  1:390994/74‑1 (MQ=255)                                                                                   cACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGcc                                                                                                             >  2:390994/1‑74 (MQ=255)                                                                                   cACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCAcc                                            >  2:1429898/1‑139 (MQ=255)                                                                                       cGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATa                                        <  2:789561/139‑1 (MQ=255)                                                                                         ccATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGAc                                                                          >  2:1170369/1‑103 (MQ=255)                                                                                         ccATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGAc                                                                          <  1:1170369/103‑1 (MQ=255)                                                                                           aTCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTc                                               >  2:477147/1‑128 (MQ=255)                                                                                           aTCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTc                                               <  1:477147/128‑1 (MQ=255)                                                                                            tCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGcc                                   <  2:1257330/139‑1 (MQ=255)                                                                                             ccccATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCg                                  >  2:811768/1‑139 (MQ=255)                                                                                                cATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATg                                                                                  <  1:649735/88‑1 (MQ=255)                                                                                                cATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATg                                                                                  >  2:649735/1‑88 (MQ=255)                                                                                                       aTGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGCGAAGAGAACGATGCCGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATccgac                        >  2:1407423/1‑135 (MQ=255)                                                                                                       aTGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCg                                                  <  2:1165225/113‑1 (MQ=255)                                                                                                       aTGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCg                                                  >  1:1165225/1‑113 (MQ=255)                                                                                                       aTGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAg                                       <  2:473318/124‑1 (MQ=255)                                                                                                       aTGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAg                                       >  1:473318/1‑124 (MQ=255)                                                                                                       aTGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGgc                        >  2:572326/1‑139 (MQ=255)                                                                                                             aGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGa                  <  1:1429898/139‑1 (MQ=255)                                                                                                               aCGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAAc                <  1:741202/139‑1 (MQ=255)                                                                                                                     gATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCga          <  1:1332390/139‑1 (MQ=255)                                                                                                                     gATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCga          <  1:811768/139‑1 (MQ=255)                                                                                                                             tGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGCTGCTgg                                                                              >  2:951053/1‑63 (MQ=255)                                                                                                                             tGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTgg                                                                              <  1:951053/63‑1 (MQ=255)                                                                                                                             tGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAAc                >  1:1311860/1‑125 (MQ=255)                                                                                                                             tGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAACATCAGGCCGCCGAACCAGCGAGAGATaaa  <  1:17938/139‑1 (MQ=255)                                                                                                                             tGGTGATGTAACCTGGCGCAAAGCCGCAGGACAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAAc                <  2:1311860/125‑1 (MQ=255)                                                                                                                              ggTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCGGGCCGCCGAACCAGCGAGAGATaaa  >  2:629096/1‑138 (MQ=255)                                                                                                                              ggTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCGGGCCGCCGAACCAGCGAGAGATaaa  <  1:629096/138‑1 (MQ=255)                                                                                                                               gTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCg                                                       <  2:279046/84‑1 (MQ=255)                                                                                                                               gTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCg                                                       >  1:279046/1‑84 (MQ=255)                                                                                                                                 gATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAg                                                    >  1:1395341/1‑85 (MQ=255)                                                                                                                                 gATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAg                                                    <  2:1395341/85‑1 (MQ=255)                                                                                                                                       |                                                                                                                                 CCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAA  >  NC_000913/3372120‑3372382

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