Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,504,498	T→C	G147G (GGT→GGC)	nanX →	sialic acid transporter NanX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,504,498	0	T	C	100.0%	81.6 / NA	23	G147G (GGT→GGC)	nanX	sialic acid transporter NanX
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/15);  total (8/15)

GGTGGTTTTTTTGGTGCCATGGCTGATAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGTTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCTTTTTTTATAGGCCTGTTACCAGTTTTACTAGTTCTTTGGATCAGAAAAAGTGCTCCAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTT  >  NC_000913/4504256‑4504744                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       ggtggtTTTTTTGGTGCCATGGCTGATAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGt                                                                                                                                                                                                                                                <  1:214871/251‑1 (MQ=255)        tttttGGTGCCATGGCTGATAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGaaa                                                                                                                                                                                                                                          <  1:61154/250‑1 (MQ=255)                             aaTATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATc                                                                                                                                                                                                                    >  1:212644/1‑251 (MQ=255)                                            cAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGAATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGa                                                                                                                                                                                                      >  1:423000/1‑250 (MQ=255)                                                   gatgTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTatat                                                                                                                                                                                               <  1:102983/250‑1 (MQ=255)                                                   gatgTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATatg                                                                                                                                                                                              >  2:299076/1‑251 (MQ=255)                                                        gggCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATatggatg                                                                                                                                                                                          <  2:137829/250‑1 (MQ=255)                                                                     ttACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCtttt                                                                                                                                                                             <  2:230818/250‑1 (MQ=255)                                                                       aCTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAAGATTATTGCGGCACAAATAAGCCCTCAGTTTGCGGAAGTATATGTATGGAGAAACTCttttttt                                                                                                                                                                          <  2:57840/251‑1 (MQ=255)                                                                           aGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCTTTTTTTATAg                                                                                                                                                                      <  1:221039/251‑1 (MQ=255)                                                                                 aaCAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCTTTTTTTATAGGCCTGt                                                                                                                                                                <  1:248549/251‑1 (MQ=255)                                                                                                   tGCTACAAACTTATATATGCTCGCAGATTGCAGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGACGTATATGGATGGAGAAACTCTTTTTTTATAGGCCTGTTACCAGTTTTACTAGTTc                                                                                                                                              >  2:378038/1‑251 (MQ=255)                                                                                                       aCAAACTTATATATGCTCGCATGTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCTTTGTTTATAGGCCGGTTACCAGTTTTACTAGTTCTTTg                                                                                                                                          <  2:423000/251‑1 (MQ=255)                                                                                                               atatatGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCTTTTTTTATAGGCCTGTTACCAGTTTTACTAGTTCTTTGGATCAGa                                                                                                                                   >  2:10787/1‑250 (MQ=255)                                                                                                                                        aTTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCTTTTTTTATAGGCCTGTTACCAGTTTTACTAGTTCTTTGGATCAGAAAAAGTGCTCCAGAAAGTCAGGAGt                                                                                                          <  2:115244/250‑1 (MQ=255)                                                                                                                                          tgttgGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTGTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCTTTTTTTATAGGCCTGTTACCAGTTTTACTAGTTCTTTGGATCAGAAAAATTGCTCCAGAAAGTCAGGAGTgg                                                                                                        <  2:185724/250‑1 (MQ=255)                                                                                                                                                   aGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAACTAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAActct                                                                                                                                                                                >  1:435117/1‑169 (MQ=255)                                                                                                                                                                 aaTATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGTACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCTGTTTTTATAGGCCTGTTACCAGTTTTACTAGTTCTTTGGATCAGAAAAAGTGCTGCAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATaaa                                                                                <  1:22725/251‑1 (MQ=255)                                                                                                                                                                         tgtgCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAActc                                                                                                                                                                                 <  2:88784/146‑1 (MQ=255)                                                                                                                                                                         tgtgCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAActc                                                                                                                                                                                 >  1:88784/1‑146 (MQ=255)                                                                                                                                                                            gCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTAAGTTTGCTGAAGTATATGGATGGAGAAACTCTTTTTTTATAGGCCTGTTACCAGTTTTACTAGTTCTTTGGATCAGAAAAAGTGCTCCAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACAttttt                                                                     >  2:237026/1‑251 (MQ=255)                                                                                                                                                                                                                                  gCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCTTTTTGTATAGGCCTGTTACCAGTTTTACTAGTTCTTTGGATCAGAAAAAGTGCGCCAGAAAGTCAGGAGTGGAGTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGGCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCATCGTCt                <  2:271947/250‑1 (MQ=255)                                                                                                                                                                                                                                               gTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCTTTTTTTATAGGCCTGTTACCAGTTTTACTAGTTCTTTGGATCAGAAAAAGTGCTCCAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTAttt  <  1:201842/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       GGTGGTTTTTTTGGTGCCATGGCTGATAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGTTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCTTTTTTTATAGGCCTGTTACCAGTTTTACTAGTTCTTTGGATCAGAAAAAGTGCTCCAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTT  >  NC_000913/4504256‑4504744

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