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

CCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCA  >  NC_000913/3372120‑3372389                                                                                                                                       |                                                                                                                                        ccTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTg                                                                                                                                     >  2:902355/1‑139 (MQ=255)   tGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGgc                                                                                                                                   >  1:533406/1‑139 (MQ=255)      aaTCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCa                                                                                                                                <  1:902355/139‑1 (MQ=255)         cAAAAAACCACTGGCTTTGTTACGCAGATGTCTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAg                                                                                                                             <  1:689513/139‑1 (MQ=255)                              aCGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGcc                                                                                                        <  1:1193138/139‑1 (MQ=255)                                gCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCga                                                                                                      <  1:856/139‑1 (MQ=255)                                            gCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGAt                                                                                          >  2:690958/1‑139 (MQ=255)                                                gCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCgagaagagaa                                                                                              <  2:615785/131‑1 (MQ=255)                                                gCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCgagaagagaa                                                                                              >  1:615785/1‑131 (MQ=255)                                                     cAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGTCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGAc                                                                                 <  1:1270673/139‑1 (MQ=255)                                                              aGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAg                                                                        >  2:184846/1‑139 (MQ=255)                                                                   gcgcTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAgacgac                                                                   >  2:1032869/1‑139 (MQ=255)                                                                    cgcTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCg                                                                                                           >  1:1255417/1‑98 (MQ=255)                                                                    cgcTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCg                                                                                                           <  2:1255417/98‑1 (MQ=255)                                                                    cgcTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGacgacg                                                                  >  1:612145/1‑139 (MQ=255)                                                                    cgcTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGacgacg                                                                  >  2:768587/1‑139 (MQ=255)                                                                      cTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGa                                                                                           <  2:991150/112‑1 (MQ=255)                                                                      cTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGa                                                                                           >  1:991150/1‑112 (MQ=255)                                                                      cTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGcc                                                                >  2:891226/1‑139 (MQ=255)                                                                        ggAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCg                                                              >  1:1155825/1‑139 (MQ=255)                                                                         gAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCgagaa                                                                                                   >  1:1065462/1‑101 (MQ=255)                                                                         gAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCgagaa                                                                                                   <  2:1065462/101‑1 (MQ=255)                                                                                aTTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTc                                                      <  1:891226/139‑1 (MQ=255)                                                                                    aCCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCAccc                                                  >  2:1256094/1‑139 (MQ=255)                                                                                             ccccATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCg                                         >  2:174852/1‑139 (MQ=255)                                                                                                  tGCCGATGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCAt                                    >  1:446543/1‑139 (MQ=255)                                                                                                       aTGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGgc                               >  1:84062/1‑139 (MQ=255)                                                                                                       aTGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGgc                               >  2:1146601/1‑139 (MQ=255)                                                                                                        tGACCAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGgcc                              >  1:964062/1‑139 (MQ=255)                                                                                                            cAGACGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGccgccg                          <  1:1146601/139‑1 (MQ=255)                                                                                                              gACGAGCGATAAACATGGTGATGTAACCTGGCGCaa                                                                                                                               >  2:480056/1‑36 (MQ=255)                                                                                                              gACGAGCGATAAACATGGTGATGTAACCTGGCGCaa                                                                                                                               <  1:480056/36‑1 (MQ=38)                                                                                                               aCGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGAc                                                                                 <  2:334802/81‑1 (MQ=255)                                                                                                               aCGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGAc                                                                                 >  1:334802/1‑81 (MQ=255)                                                                                                               aCGAGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAAc                       <  1:174852/139‑1 (MQ=255)                                                                                                                  aGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATg                                                                                         >  1:990817/1‑70 (MQ=255)                                                                                                                  aGCGATAAACATGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATg                                                                                         <  2:990817/70‑1 (MQ=255)                                                                                                                             tGGTGATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATaaa         <  1:962384/139‑1 (MQ=255)                                                                                                                                 gATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTcc                                                                                                             >  2:1142600/1‑35 (MQ=255)                                                                                                                                 gATGTAACCTGGCGCAAAGCCGCAGGCCAGCGTcc                                                                                                             <  1:1142600/35‑1 (MQ=255)                                                                                                                                    gTAACCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCa  <  2:533406/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                        CCTGAAATCAAAAAACCACTGGCTTTGTTACGCAGATGTTTTGGCCAGCTTTCAATGACATAGGTGGCGCTGGAACCGTATTCACCCGCCATCCCCATGCCGATGACCAGACGAGCGATAAACATGGTGATGTAGCCTGGCGCAAAGCCGCAGGCCAGCGTCCCGGCCGAGAAGAGAACGATGCTGGTGACCATTGCCAGACGACGCCCGTAGCGGTCACCCATAGCGCCGAGCATCAGGCCGCCGAACCAGCGAGAGATAAAGGCTGCA  >  NC_000913/3372120‑3372389

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