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

CCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537962‑4538258                                                                                                                                                      |                                                                                                                                                      ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCt                                                                                                                                                    >  2:237196/1‑151 (MQ=255)   ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCttt                                                                                                                                                  >  1:760857/1‑151 (MQ=255)     gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTAt                                                                                                                                                >  2:53181/1‑151 (MQ=255)                         aTCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGcc                                                                                                                            >  2:339768/1‑151 (MQ=255)                                      aCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAgg                                                                                                               >  2:860869/1‑151 (MQ=255)                                          gAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                            >  2:781699/1‑150 (MQ=255)                                           aGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAt                                                                                                          >  2:290985/1‑151 (MQ=255)                                               aGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                                      >  1:792509/1‑151 (MQ=255)                                               aGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                                      >  1:515703/1‑151 (MQ=255)                                                gCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACAGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                                      >  2:29970/1‑150 (MQ=255)                                                 cTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     >  2:702902/1‑150 (MQ=255)                                                  ttATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     >  1:824052/1‑149 (MQ=255)                                                   tATTCCATTTTAAATTATTCCCGGTAAAATAAATTTCAAATACGGCATCCGTTCGCAATCCTGTTT‑T‑GGCTTCGCCATAAATAAGCCAGGTTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAg                                                                                                   <  1:636789/150‑1 (MQ=255)                                                   tATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                                  >  1:166313/1‑151 (MQ=255)                                                     tcccATTTAAAATTATTCCCGGTAAAATCAAGTTAAAATCCGCCAT‑CGTTTCCCATCCTTGTT‑TGGGCTTCCCCATAAATAACCCGGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg                                                                                                 <  2:174896/148‑1 (MQ=255)                                                         aTTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt                                                                                            >  1:16904/1‑151 (MQ=255)                                                              aaaTTTTTCCCGGTGAAATCAATTTCAAATACGGCACCCGTTCGCAATCCGGTTT‑T‑GGTTCCCCCATTAATAACCCGGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                       <  1:57978/151‑1 (MQ=255)                                                              aaaTTTTTCCCGGTAAAACCAAGTTAAAATCGGGCACCCGTTCCCACCCCTGTTT‑T‑GTTTCCGCCATAAAAAAGCCGGTTTTTATCCCCTTTTTCCCCAACCGCCGCTCCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                       <  2:234423/151‑1 (MQ=255)                                                                       cgggTAAAATAAAGTTAAAATGCGGCTTCCGTTCGCAATCTTGTTT‑T‑GGCTTCCCAATTAATAAGCCGGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                              <  2:678528/149‑1 (MQ=255)                                                                       ccGGTAAAATCAATTTCAAATACGGCACCCTTTCGCAACCCTGTTT‑G‑GGCTCCCCCATTAATAAGCCAGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                               <  1:702/150‑1 (MQ=255)                                                                           tGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                          >  2:392004/1‑151 (MQ=255)                                                                             aaaccaattaaaaaTACGGCCTCTG‑TCCCAACCCCGGTTGG‑GGCTGCGCCATTAAAAACCCGGTTTTTATCCCTTTTATACAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGa                                                                         <  2:205495/139‑1 (MQ=255)                                                                             aaaTCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                        >  2:426865/1‑151 (MQ=255)                                                                              aaacAATTTAAAATACTGCCTCCTTTCGCAATCTGGTTT‑T‑GGTTCCTCCATCAATAAGCCAGTTTTTATCCCCTTTATTCAAATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                                       <  2:853383/148‑1 (MQ=255)                                                                              aaTAAAGTTAAAATCCGGCCCCCTTTCCCATCCCTGTTT‑T‑GGCTCCCCCTTTAAAAACCCAGTTTTTATCCCTTTTTTCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                                       <  2:54631/151‑1 (MQ=255)                                                                               aTCAAGTTAAAATACGGCATCCGTTCGCAATCCTGTTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                      <  1:568186/151‑1 (MQ=255)                                                                                           tACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                             >  1:184990/1‑148 (MQ=255)                                                                                              ggCATCCGTCCGCAACCCTGTTT‑T‑GGCTTCGCCATTAATAAGCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                                        <  1:152107/150‑1 (MQ=255)                                                                                                cttccttccccAATCCGGTTT‑T‑GGCTCCCCAATTAATAACCCAGTTTTTACCCCCTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                                        <  2:605640/138‑1 (MQ=255)                                                                                                    cGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                                 >  1:460988/1‑151 (MQ=255)                                                                                                      ttCGCAACCCTGTTT‑T‑GGTTTCGCAATTAATAAGCCGGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                                 <  1:3697/149‑1 (MQ=255)                                                                                                      tccccattctggttt‑t‑GGTTCCCCCATTAAAAAGCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                          <  2:646924/127‑1 (MQ=255)                                                                                                      tccGCAATCCTGTTT‑T‑GGCTCCGCCATTAAAAAGCCAGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                                 <  1:658994/147‑1 (MQ=255)                                                                                                       tCGCAACCTTGTTT‑T‑GGCTTCGCCATTAAAAAGCCAGGTTTTACCCCCTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaa                                               <  2:309503/150‑1 (MQ=255)                                                                                                        cGCAACCCTGGTT‑T‑GGCTTCGCCATTAAGAAGCCAGGTTTTACCCCTTTAATGCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                             <  1:630231/151‑1 (MQ=255)                                                                                                          cactccggttt‑t‑GCCTGCGCAATAAATAACCCAGTTTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                             <  2:201965/140‑1 (MQ=255)                                                                                                          cAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaagaa                                            >  1:125079/1‑150 (MQ=255)                                                                                                           aaTCCTGTTT‑T‑GGCTTCGCCATTAATAGCCCAGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaagaa                                            <  1:678545/149‑1 (MQ=255)                                                                                                                 ggTT‑T‑GGCTTCGCCATAAATAAGCCAGTTTTTACCCCCTTTTTTCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct                                     <  1:310787/150‑1 (MQ=255)                                                                                                                   tt‑t‑GCCTCCCCCATAAATAACCCAGTTTTTACCCCCTTTATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtc                                      <  2:231374/147‑1 (MQ=255)                                                                                                                    t‑t‑GGTTTCCCCATAAATAACCCAGTTTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                                  <  1:194861/150‑1 (MQ=255)                                                                                                                            tCGCCATTAATAACCCAGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                            <  1:784694/150‑1 (MQ=255)                                                                                                                                      taaGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                  >  1:839332/1‑150 (MQ=255)                                                                                                                                       acgCCGGGTTTTATCGCCTTTATCCCCATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                      <  2:267998/143‑1 (MQ=255)                                                                                                                                           cAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg              >  2:476405/1‑149 (MQ=255)                                                                                                                                           cAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg              >  2:822520/1‑149 (MQ=255)                                                                                                                                               ttttATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc           <  2:121694/148‑1 (MQ=255)                                                                                                                                               ttttACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc           <  1:382943/148‑1 (MQ=255)                                                                                                                                                   aTCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtt    >  1:694153/1‑151 (MQ=255)                                                                                                                                                      aCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  2:91052/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTT‑T‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537962‑4538258

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: ‑