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

CCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTG‑AAATCAAGTTCAAATACGGCATCCGTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATT  >  NC_000913/4537962‑4538252                                                                                                                                                        |                                                                                                                                                ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTG‑AAATCAAGTTCAAATACGGCATCCGTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCt                                                                                                                                              <  2:715856/151‑1 (MQ=255)           ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTG‑AAATCAAGTTCAAATACGGCATCCGTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACaa                                                                                                                                    <  2:517471/151‑1 (MQ=255)            cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTG‑AAATCAAGTTCAAATACGGCATCCGTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAc                                                                                                                                   <  2:25355/151‑1 (MQ=255)                      gCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTG‑AAATCAAGTTCAAATACGGCATCCGTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCAcc                                                                                                                          >  2:117971/1‑150 (MQ=255)                                 atgaGACGGGAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTG‑AAATCAAGTTCAAATACGGCATCCGTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGt                                                                                                               >  1:56899/1‑150 (MQ=255)                                        gggAGCAAGTTTTTTCCATTTAAAATTATTCCG‑GGTA‑AAATAAATTTCAAATCCGGCATCCGTCCGAAATCC‑TGTTTTGG‑CTCCGCCATTAAAAACCCGGGTTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAgg                                                                                                         <  1:237869/149‑1 (MQ=255)                                        gggAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTG‑AAATCAAGTTCAAATACGGCATCCGTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                       >  2:108318/1‑151 (MQ=255)                                        gggAGAAGGTTTATCCCATTTTAATTTATTCCC‑GGAA‑AAATAAAGTTCAAATCCGCCACCCGTTCCCAATCC‑TGTTTGGG‑CTCCGCCTTAAATAGCCCGGTTTTTACCCCCTTTTTTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                       <  1:777158/151‑1 (MQ=255)                                              aaGCTTATTCCATTTTAAATTATTACC‑GGTG‑AAATCAAGTTCAAATACGGCATCCGTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                 >  2:88427/1‑151 (MQ=255)                                                 cTTATTCCATTTTAAATTATTACC‑GGTG‑AAATCAAGTTCAAATACGGCATCCGTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                               >  2:768446/1‑150 (MQ=255)                                                   ttatcctatttAAATTTATCCC‑GGTGAAAAACAATTTAAAACACGCCAGCGTTTACCACCCC‑TTGTTTGG‑TTCCCCCTTAAATACGCAAGGTTTCACCCCCTTAATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCa                                                                                              <  2:599891/142‑1 (MQ=255)                                                   tATTCCATTTTAAATTATTACC‑GGTG‑AAATCAAGTTCAAATACGGCATCCGTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                            >  1:688497/1‑151 (MQ=255)                                                     ttCCTTTTTAACTT‑TTACCGGGTG‑ATAAAAAGTAAAAATACGGCACCCGCTCTACTTCC‑CGTTTAGA‑CTCCACAATTAATCACCTCGTTTTTACCACTTTAATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg                                                                                           <  1:130459/150‑1 (MQ=255)                                                     ttCCATTTTAAATTATTACC‑GGTG‑AAATCAAGTTCAAATACGGCATCCGTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                          >  1:535198/1‑151 (MQ=255)                                                     ttCCATTTTAAATTATTACC‑GGTG‑AAATCAAGTTCAAATACGGCATCCGTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                          >  2:5204/1‑151 (MQ=255)                                                          ttttAAATTATTACC‑GGTG‑AAATAAAGTTAAAATACGCCATCCGTTCCCAATCC‑TGGTTGGG‑CTTCGCCATTAATAAGCCAGTTTTTATCCCCTTTATTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                         <  1:540183/147‑1 (MQ=255)                                                                           gggaaaataaatTTAAAATAC‑GCCTCCCTTCCCAAACCTTGTTTTGT‑TTTCGCCTTAAAAAAGCCAGTTTTTATCCCCTTTATCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGa                                                                       <  2:431802/146‑1 (MQ=255)                                                                                 atcaattaaaaaTACGGCACCCTTTCGCACCCC‑TGTTTTGG‑CTTCGCCATTAATAAGCCGGTTTTTATCCCCTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                <  1:662552/142‑1 (MQ=255)                                                                                   cAAGTTCAAATACGGCATCCGATCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                               >  2:599844/1‑150 (MQ=255)                                                                                    aatttCAAATCCGGCATCCTTTCGCAACCC‑TGTTTGGG‑CTTCGCCATAAATAACCCGGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                               <  1:842848/146‑1 (MQ=255)                                                                                    aaGTTCAAATACGGCATCCGTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                              >  2:169335/1‑150 (MQ=255)                                                                                          aaaTGCGGCACCCGTTCGAAACCC‑TGTTTTGG‑TTCCGCCATAAATAAGCCAGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGa                                                           <  1:8038/147‑1 (MQ=255)                                                                                           aaTACGGCATCCGTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                       >  1:369034/1‑150 (MQ=255)                                                                                               cGGCATCCGTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                    >  2:639320/1‑149 (MQ=255)                                                                                                ggCATCCGTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                                 >  1:752362/1‑151 (MQ=255)                                                                                                ggCATCCGTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                                 >  2:207834/1‑151 (MQ=255)                                                                                                gccacccGTTCGCACCCT‑TGTTTTGG‑CTCCCCCATAAATAAGCCAGTTTTTATCCCTTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                                  <  1:201907/145‑1 (MQ=255)                                                                                                     ccGTTCGCAACCT‑TGTTTTGG‑CTTCCCAATAAAAAAGCCAGTTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                               <  1:338217/148‑1 (MQ=255)                                                                                                       gTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                           >  1:125379/1‑150 (MQ=255)                                                                                                          cGCAATCC‑TGGTTTGG‑CTCCGCCATTAATAAGCCAGTTTTTACCCCCTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                        <  1:478006/150‑1 (MQ=255)                                                                                                          cGCAACCC‑TGTTTTGC‑CTTCCCCATAAAAAACCCGGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                               <  1:571356/143‑1 (MQ=255)                                                                                                           cccatccttGGTTGGC‑CTCCGCCTTAAAAAGGCCAGTTTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                            <  2:433775/138‑1 (MQ=255)                                                                                                            cAACCC‑TGTTTTGC‑CTTCCCAATAAAAAAGCCTGTTTTTACCCCCTTTATGCCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAATGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaagaa                                      <  2:75618/150‑1 (MQ=255)                                                                                                             aaTCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                    >  1:613268/1‑151 (MQ=255)                                                                                                             aaTCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                    >  1:167782/1‑151 (MQ=255)                                                                                                              aTCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                    >  2:299295/1‑150 (MQ=255)                                                                                                               attccggTTTGGCCTCCGCAATTAAAAACCCAGTTTTTATCCCCTTTATCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                       <  2:131787/143‑1 (MQ=255)                                                                                                                  cgttgttggccTTCCCAATTAAAACGCCTGGTTTTATCCCCTTAATCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataa                                  <  2:58085/140‑1 (MQ=255)                                                                                                                     gttttgctttCGCCTTTAAAAAGCCGGGTTTTACCCCTTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct                               <  2:659331/142‑1 (MQ=255)                                                                                                                      tttGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                            >  1:698524/1‑151 (MQ=255)                                                                                                                      tttGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                            >  2:639579/1‑151 (MQ=255)                                                                                                                       ttGC‑CTCCGCAATAAATAAGCCAGTTTTTACCCCCTTTATACACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa                                     <  2:647008/141‑1 (MQ=255)                                                                                                                                 ccATAAAAAACCCAGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                  <  1:839995/151‑1 (MQ=255)                                                                                                                                    taaaaaacccAGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                <  1:86425/142‑1 (MQ=255)                                                                                                                                       aaaccccgttttttATCCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt            <  1:848063/141‑1 (MQ=255)                                                                                                                                         aacccgttttttatcccttttatcccAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATa           <  1:624847/142‑1 (MQ=255)                                                                                                                                              aGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc     >  2:361469/1‑151 (MQ=255)                                                                                                                                               ggTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc     >  2:700009/1‑150 (MQ=255)                                                                                                                                                gTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAt   >  2:199535/1‑151 (MQ=255)                                                                                                                                                 ttttATCCCTTTTATCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt         <  1:894861/144‑1 (MQ=255)                                                                                                                                                 ttttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCTATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAtt  >  2:477989/1‑151 (MQ=255)                                                                                                                                                 ttttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAtt  >  1:465420/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                CCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACC‑GGTG‑AAATCAAGTTCAAATACGGCATCCGTTCGCAATCC‑TGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATT  >  NC_000913/4537962‑4538252

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