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	94.1%	53.4 / ‑3.5	17	G147G (GGT→GGC)	nanX	sialic acid transporter NanX
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (4/12);  total (4/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

ATGGCTGATAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGTTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCTTTTTTTATAGGCCTGTTACCAGTTTTACTAGTTCTTTGGATCAGAAAAAGTGCTCCAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTT  >  NC_000913/4504274‑4504720                                                                                                                                                                                                                                 |                                                                                                                                                                                                                               aTGGCTGATAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGc                                                                                                                                                                                                      <  1:162547/251‑1 (MQ=255)          aaaTATTGTCGTAAGCCAATGCGGATGTGGTCAATTTTCATTTTCTCATGCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGTTAGCAGTTTGCCGATTTATTGTTGGCTTATGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGTCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCTCAtataat                                                                                                                                                                                             <  2:215660/251‑1 (MQ=255)           aaTATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAaataata                                                                                                                                                                                            <  2:214696/251‑2 (MQ=255)                ggTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCa                                                                                                                                                                                       <  1:129723/251‑1 (MQ=255)                        gCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTg                                                                                                                                                                               <  1:268949/251‑1 (MQ=255)                                        gCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGgag                                                                                                                                                               <  1:129024/251‑1 (MQ=255)                                                cATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAAGCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCtt                                                                                                                                                       <  1:143691/251‑1 (MQ=255)                                                                aCAGGCCTTAGCGGTATTGCTACAAACTTATATATGCGCGCAGTTTGCCGGTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCTTTTTTTATAGGCCTGtt                                                                                                                                       <  1:153215/251‑1 (MQ=255)                                                                         aGCGGTATTGCTACAAGCTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATTTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGGGCATTTTTGGGAAGTTGTTTTTCTGTTGGAAATATGATGGCGGCACAAATAATCCCTCAGTTGGCTGAAGTATATGGATGTAGAAACTCTCTTTTTATAGGCCTGTTAGCAGtttt                                                                                                                              <  2:18302/251‑1 (MQ=255)                                                                         aGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCTTTTTTTATAGGCCTGTTACCAGttt                                                                                                                               <  1:211370/250‑1 (MQ=255)                                                                               aTTGCTACAAACTTATATATGCGCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTTCTGAAGTATATGGATGGAGAAACTCTTTTTTTATAGGCCTGTTACCAGTTTTACTAGt                                                                                                                        <  2:87856/251‑1 (MQ=255)                                                                                                             tGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAAtatta                                                                                                                                                                                                             >  1:230097/1‑136 (MQ=255)                                                                                                             tGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAAtatta                                                                                                                                                                                                             <  2:230097/136‑1 (MQ=255)                                                                                                                         gttgGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCTTTTTTTATAGGCCTGTTACCAGTTTTACTAGTTCTTTGGATCAGAAAAAGTGCTCCAGAAAGTCAGAAGTGGa                                                                               <  1:74390/250‑1 (MQ=255)                                                                                                                                                    gCATGTGCTTAAACTGATGCGGACGAACGTTGGCCTAAACATCTTCAATCTCAACCTAGTGCTTTTTACGAAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCACTCAGTTTGCTGAAGTCTATGGATGGAGAAACTCTTTTTTTATAGGCCTGTTACCAGTTTTACTAGTTCTTTGGATCAGAAAAAGTGCTCCAGAAAGCCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAAc                                                   >  2:164637/1‑251 (MQ=255)                                                                                                                                                                                                 cAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCTTTTTTTATAGGCCTGTTACCAGTTTTACTAGTTCTTTGGATCAGAAAAAGTGCTCCAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGAt      >  2:94310/1‑251 (MQ=255)                                                                                                                                                                                                     cTAAAGCTAGTGCTTTTTTGGTAAGTGGCTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCTTTTTTTATAGGCCTGTTACCAGTTTTACTAGTTCTTTGGATCAGAAAAAGTGCTCCAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGtt  >  1:193094/1‑251 (MQ=255)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                               ATGGCTGATAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGGGATGTCTGGTGAATATGCATGTGCTTCAACTTATGCGGTAGAAAGTTGGCCTAAAAATCTTCAATCTAAAGCTAGTGCTTTTTTGGTAAGTGGTTTTTCTGTTGGAAATATTATTGCGGCACAAATAATCCCTCAGTTTGCTGAAGTATATGGATGGAGAAACTCTTTTTTTATAGGCCTGTTACCAGTTTTACTAGTTCTTTGGATCAGAAAAAGTGCTCCAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTT  >  NC_000913/4504274‑4504720

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