Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,109	0	A	C	42.5%	81.5 / 12.1	48	G219G (GGT→GGG)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (23/4);  new base C (0/20);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.27e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTG  >  NC_000913/4537960‑4538253                                                                                                                                                          |                                                                                                                                                  aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCAc                                                                                                                                                  >  2:568107/1‑151 (MQ=255)       gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTAt                                                                                                                                           >  2:363136/1‑151 (MQ=255)         aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTAtt                                                                                                                                          <  2:232038/150‑1 (MQ=255)         aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTc                                                                                                                                         <  2:427983/151‑1 (MQ=255)            aCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTcaca                                                                                                                                      <  1:270291/151‑1 (MQ=255)                  aGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAAccgcc                                                                                                                                >  1:748152/1‑151 (MQ=255)                     tACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACcgccgc                                                                                                                              >  2:382175/1‑150 (MQ=255)                       cGCCATCTGTTGATGAGACGTGAGCAAGCTTCTTCCTTTTTAAATTATTACC‑GGTGAAATCAACTTCAAATACGGCATCCCTTCGCCATCCTGGTTC‑GG‑CCTCGCCTTTAAAAAGCCAGG‑TTTTATCCCC‑TTTATTCACAACCGCCGCAc                                                                                                                            <  1:49584/150‑1 (MQ=255)                       cGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCAcc                                                                                                                           >  1:264673/1‑151 (MQ=255)                        gCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCAcc                                                                                                                           >  1:663623/1‑150 (MQ=255)                                        aCGGGAGCAATCTTTTTCCATTTTAAATTATTACG‑GGTGAAATAAAGTTAAAATCCGGCACCCTTTCGCAATCCTGTTTT‑GG‑CTCCCCTATTAATAAGCCAGT‑TTTTATCCCC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAg                                                                                                           <  2:795796/150‑1 (MQ=255)                                                   cTTATTCCATTTTAAATTATTACC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                >  2:466876/1‑150 (MQ=255)                                                    tttttCCTTTTAAATTTTTTACC‑GGTGAAATAAAGTCAAAATCCGGC‑TCCGTTCCACATCCCTGTTT‑GGTTTTCGCCATTAATAA‑CCAGGTTTTTATCCCC‑TTTTTCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCa                                                                                               <  1:586497/147‑1 (MQ=255)                                                        ttccttttAAATTTTTTCCGGGTAAAATAAATTAAAAATACGGCA‑CCCTTCCCCACCCCGGTTTGGG‑CTTCCCCA‑TAAAAAGGCCGG‑TTTTAACCCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg                                                                                            <  1:363114/146‑1 (MQ=255)                                                         ccATTTTAAATTATTACC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                          >  2:182377/1‑150 (MQ=255)                                                            ttttAAATTATTACC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTc                                                                                      >  2:238865/1‑151 (MQ=255)                                                             tttAAATTATTACC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                     >  1:273835/1‑151 (MQ=255)                                                                aaaTTATTCCC‑GGTAAAAACAATTTCAAATGCGGAATCCGTTCGCATCCCTGGTTG‑GG‑TTTCGCAATTAAAAATCCTGT‑TTTTACCCCT‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                  <  2:734479/151‑1 (MQ=255)                                                                aaaTTATTACC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                  >  2:435445/1‑151 (MQ=255)                                                                aaaTTATTACC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                  >  2:546231/1‑151 (MQ=255)                                                                               gAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                     >  1:739219/1‑150 (MQ=255)                                                                                aaataaatTTAAAATACGACACCCTTTCCCATCCCTGTTTT‑GG‑CTCCCCCTTTAATAAGCCGGT‑TTTTACCCCT‑TTTATCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg                                                                      <  1:242739/148‑1 (MQ=255)                                                                                        ttCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                           >  2:258340/1‑151 (MQ=255)                                                                                           aaaTCCGGCTCCCTTTCGCAATCCTGGTTG‑GC‑CTTCCCCATTAATAAGCCGGT‑TTTTATCCCT‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                          <  2:553669/149‑1 (MQ=255)                                                                                                 gcctcccttccGCACTCCTGTTTT‑GG‑CTCCGCCATAAAAACCCCGGT‑TTTTACCCCC‑TTTTTCCAAACCCGCCGCACCAGCCGTTCCGTACCAGAGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                    <  1:679564/139‑1 (MQ=255)                                                                                                   cATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                >  2:205065/1‑151 (MQ=255)                                                                                                   cATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                >  2:155191/1‑151 (MQ=255)                                                                                                    acccGTCCGCACCCCTGGTTT‑GG‑TTTCGCATTTAAAAGGCCAGT‑TTTTACCCCT‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAg                                                      <  1:628136/142‑1 (MQ=255)                                                                                                        gTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                            >  2:366650/1‑150 (MQ=255)                                                                                                            gCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                        >  2:476003/1‑150 (MQ=255)                                                                                                             catccttgtttt‑GC‑CTTCGCCATTAAAAAGCCAGG‑TTTTATCCCC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaagaa                                       <  2:33890/141‑1 (MQ=255)                                                                                                              aaTCCTGGTTG‑GT‑TTTCCCTATAAAAAACCCGGT‑TTTTACCCCC‑TTTATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                         <  2:596034/147‑1 (MQ=255)                                                                                                                 cttGGTTT‑GG‑CTCCGCAATTAATAAGCTAGT‑TTTTTCCACC‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                     <  2:313752/146‑1 (MQ=255)                                                                                                                  cTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtc                                 >  2:688963/1‑151 (MQ=255)                                                                                                                    gtttt‑GG‑CTTCGCCATTAATAAGCAAGG‑TTTTATCTCT‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt                               <  2:15141/149‑1 (MQ=255)                                                                                                                         ttgggctcccctttAAAAAACCCGG‑TTTTTCCCCCTTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct                                <  1:738926/136‑1 (MQ=255)                                                                                                                              ttCGCCATAAATAAGCCAGG‑TTTTATCCCT‑TTTATACACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtt                         <  2:494119/149‑1 (MQ=255)                                                                                                                              ttCGCAATAAATAAGCCGGT‑TTTTACCCCT‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                            <  1:576653/146‑1 (MQ=255)                                                                                                                                cGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                     >  1:362499/1‑151 (MQ=255)                                                                                                                                  caataaaaagcccAGT‑TTTTCTCCCT‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                       <  1:99295/136‑1 (MQ=255)                                                                                                                                      taataaGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa                >  2:591925/1‑150 (MQ=255)                                                                                                                                       aaaaaGCCAGT‑TTTTACCCCT‑TTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                 <  2:359723/145‑1 (MQ=255)                                                                                                                                          aaGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaaa              >  2:812819/1‑148 (MQ=255)                                                                                                                                           aaccaagtttttACCCCC‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAg           <  2:255382/143‑1 (MQ=255)                                                                                                                                           aGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAg           >  1:728203/1‑150 (MQ=255)                                                                                                                                                 gtttttACCCCC‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg         <  1:372814/145‑1 (MQ=255)                                                                                                                                                   ttttATCCCC‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt          <  2:581228/144‑1 (MQ=255)                                                                                                                                                     ttATCCCT‑TTTATCAACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg  <  1:584627/150‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                                  ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTG  >  NC_000913/4537960‑4538253

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