Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,117	0	T	C	43.5%	50.6 / 15.2	46	N217D (AAT→GAT)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base T (19/7);  new base C (0/20);  total (19/27)
Fisher's exact test for biased strand distribution p-value = 2.49e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTAT  >  NC_000913/4537972‑4538264                                                                                                                                                      |                                                                                                                                                     ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑caca                                                                                                                                                 >  2:561308/1‑150 (MQ=255)         tACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACcgccgc                                                                                                                                         >  2:19988/1‑150 (MQ=255)          aCGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCAc                                                                                                                                       >  2:476105/1‑151 (MQ=255)                        tgaGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTAc                                                                                                                         >  2:259963/1‑151 (MQ=255)                         gagaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                                        >  1:563502/1‑151 (MQ=255)                              gggAGCAAGTTTTTCCCTTTTAAAATTATCCCCGGTAAAAC‑AAATTAAAAATACGCCA‑CCCTTCCCAACTCTTGTTTTGC‑CTCCCCCATAAAAAGCCCAGT‑TTTTATCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAggg                                                                                                                    <  2:55223/150‑1 (MQ=255)                              gggAGCAAGTTTATTCCTTTTAAAATTATTCCCGGTAAAAT‑AAAGTTAAAATCCGGCATCCGTTCGCAA‑CCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGT‑TTTTATCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAg                                                                                                                      <  1:119450/148‑1 (MQ=255)                              gggAGCAAGTTTATTCCTTTTAAAATTATTACCGGTGAAAT‑AAAGTTAAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATAAAAAACCCAGT‑TTTTATCCCCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAgg                                                                                                                     <  1:586003/149‑1 (MQ=255)                              gggAGCAACTTTATTCCTTTTAAATTTTTTCCCGTTAAAAT‑CAAGTTAAAATCCGGCATCCTTCCGC‑A‑TCCCTGTTTGG‑GCTCCCCAATAAATAGCCAGGTTTTTATC‑CCTTTATTCCACACCCGCCGCACCAGCCGTTCCGTACCAggg                                                                                                                    <  2:624051/150‑1 (MQ=255)                               ggAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                                  >  2:315543/1‑151 (MQ=255)                               ggAGCAA‑CTTTTTTCATTTTAAATTTTTACCGGGAAAATCCAAGTTAAAAGCGGGCACCCGTTCCCAA‑‑CCTTGGTTGG‑GTTCCCCCTTAAAAAACCCGGTTTTTACCCCTTTTTTT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                                  <  2:7765/151‑1 (MQ=255)                                gAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                                   >  2:488147/1‑149 (MQ=255)                                 aGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAt                                                                                                                >  1:433282/1‑151 (MQ=255)                                            tCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                                      >  2:15136/1‑150 (MQ=255)                                                 tttAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                                >  2:651618/1‑151 (MQ=255)                                                             ccGGTAAAAT‑CAAGTTAAAATACGGCATCCGTTCGCAA‑TCCGGTTTTGG‑CTTCCCCATTAAAAAGCCAGG‑TTTTATCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                                     <  1:109589/150‑1 (MQ=255)                                                                 tGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                                >  1:426669/1‑151 (MQ=255)                                                                 gttaaaaaaaattCAAAACCCGCCTCCCTTCCAAA‑ACCTGGTTTGGCCTGCGCCATAAAAAAGCCGGG‑TGTTATCCCTTTAATC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                                     <  2:339532/137‑1 (MQ=255)                                                                   aaaT‑CAAGTTAAAATCCGGCACCCTTTCGCAA‑TCCTGTTTTGG‑TTTCGCCATAAAAAACCCAGT‑TTTTACCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                              <  1:539650/151‑1 (MQ=255)                                                                           gtacAATTACGGCATCCGTCCGAAT‑TCGTGGTTGGG‑CTCCGCCATTAAAAAGCCGGT‑TTTTCCCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAg                                                                        <  1:212667/147‑1 (MQ=255)                                                                                   ccgcatccctttCGCAT‑CCTTGTTTTGC‑CTCCCCCATAAATAACCCGGG‑TTTTACCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                               <  2:528100/141‑1 (MQ=255)                                                                                   actgtaTCCCTTCGAAA‑TCTTGTTTTGG‑TTGCCCCATAAAAACGCCAGT‑TTTTACCCCTTTTATC‑CCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAg                                                                 <  1:573871/144‑1 (MQ=255)                                                                                          ccGTCCGCAA‑TCCTGTTGTGG‑CTCCGCCATAAAAAAGCCGGG‑TTTTCCCCCCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAt                                                         <  2:543384/150‑1 (MQ=255)                                                                                            gctCGCAC‑TCC‑GGTTGGGCCTTCCCCATTAAAAAGCCAGT‑TTTAATCCCTTTTTTC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                                      <  2:100224/149‑1 (MQ=255)                                                                                                 cAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaagaa                                                  >  1:348421/1‑150 (MQ=255)                                                                                                  aa‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa                                                 >  2:248544/1‑150 (MQ=255)                                                                                                   aacccggttttGG‑CTCCCCCATAAAAAAGCCGGT‑TTTTATCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                                   <  1:709403/140‑1 (MQ=255)                                                                                                       cTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtc                                            >  1:140674/1‑151 (MQ=255)                                                                                                       cTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtc                                            >  2:530934/1‑151 (MQ=255)                                                                                                          ttttGC‑CTTCGCCATAAAAAAGCCGGT‑TTTTACCACCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt                                          <  1:343076/149‑1 (MQ=255)                                                                                                              gg‑CTTCCCCATTAATAGGCCGGT‑TTTTATCCCTTTTATC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                       <  1:316061/149‑1 (MQ=255)                                                                                                                    cGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                                >  1:491488/1‑151 (MQ=255)                                                                                                                    cGCCATAAAAACGCCAGT‑TTTAACCCCCTTTTTC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                                <  2:428070/151‑1 (MQ=255)                                                                                                                           aaaaaGCCGGT‑TTTTATCCCCTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa                           <  1:299644/146‑1 (MQ=255)                                                                                                                            aaaaggccggtttttAGCCCTTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                          <  2:210635/139‑1 (MQ=255)                                                                                                                            aaaaGCCAGT‑TTTTATCACTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                        <  2:499143/149‑1 (MQ=255)                                                                                                                             aaccccAGT‑TTTTCCCCCCTTTATC‑CCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                        <  2:70576/146‑1 (MQ=255)                                                                                                                              aaGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                        >  1:574799/1‑149 (MQ=255)                                                                                                                                 ccAGT‑TTTTATCCCCTTAATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATa                       <  2:246137/147‑1 (MQ=255)                                                                                                                                   aGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc                 >  2:367066/1‑151 (MQ=255)                                                                                                                                   aGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc                 >  1:397627/1‑151 (MQ=255)                                                                                                                                     gtttttATCTCTTTTATT‑CCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc                 <  1:669663/148‑1 (MQ=255)                                                                                                                                        tttATCCCCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAtt              <  1:215474/150‑1 (MQ=255)                                                                                                                                         ttctCACTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg             <  2:55831/147‑1 (MQ=255)                                                                                                                                             cccttttATT‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAt           <  1:417234/144‑1 (MQ=255)                                                                                                                                                  ttATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCt    <  1:412388/150‑1 (MQ=255)                                                                                                                                                     tt‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCtat  >  1:8002/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTAT  >  NC_000913/4537972‑4538264

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