Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,657	0	A	G	42.9%	64.0 / 11.0	49	F386S (TTC→TCC)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (6/22);  new base G (21/0);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 7.74e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGC  >  NC_000913/3508507‑3508805                                                                                                                                                         |                                                                                                                                                      gcAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAAACCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGCCCATTTTGTGGCGGT‑GGCGCG‑CAGa                                                                                                                                                       <  2:345762/151‑1 (MQ=255)        aCGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGGGCG‑CGGa‑aaaaat                                                                                                                                                >  1:846021/1‑148 (MQ=255)         cGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGg‑aaaaat                                                                                                                                                >  1:707865/1‑147 (MQ=255)         cGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTca                                                                                                                                              <  2:856092/151‑1 (MQ=255)         cGTTGTTGTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTGTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑GGGA‑AAAAtt                                                                                                                                                <  2:861647/149‑1 (MQ=255)             gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGG‑AAAAATcaca                                                                                                                                            >  1:818659/1‑149 (MQ=255)                tttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACg                                                                                                                                        <  2:880452/150‑1 (MQ=255)                 ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGG‑AAAAATcacaaa                                                                                                                                          >  1:160591/1‑145 (MQ=255)                    aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGc                                                                                                                                   <  1:1045404/151‑1 (MQ=255)                        cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGGGCG‑CGGG‑AAAAATCACACACGGGAGccccc                                                                                                                               >  1:440641/1‑151 (MQ=255)                        cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGccccc                                                                                                                               <  2:440465/151‑1 (MQ=255)                            gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGGGCG‑GGGG‑AAAAATCAAAAACgggggacccc                                                                                                                               >  2:55587/1‑139 (MQ=255)                            gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACAt                                                                                                                           <  1:163178/151‑1 (MQ=255)                               aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGGGCG‑CGGG‑AAAAATCACACACGGGAGCCCCcaca                                                                                                                            >  1:731669/1‑147 (MQ=255)                                    tatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑TGGGCG‑CGGG‑AAAAATCAAACCCCGGAGGCCCCCCAACAaaaa                                                                                                                     >  1:886861/1‑149 (MQ=255)                                        cAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGG‑AAAAATTAAAAACAGGAGCCCCCACAACAGAAAGCgggg                                                                                                               >  2:272858/1‑151 (MQ=255)                                        cAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGG‑AAAAATAACACACGGGAGCCCCCACATCAGAAAAAgtgg                                                                                                               >  2:412902/1‑151 (MQ=255)                                        cAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCgt                                                                                                                 <  2:542109/149‑1 (MQ=255)                                           aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGG‑AAAAATAAAACAAAGTAGCCCCCCCATCCGCAAgagggg                                                                                                               >  1:824998/1‑143 (MQ=255)                                                cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGca                                                                                                         <  1:1028482/149‑1 (MQ=255)                                                 aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGAG‑CGGAGAAAATTAAAAAACGGACGCCCCC‑CATCTGaaaaagaggggcggcaa                                                                                                       >  1:875073/1‑136 (MQ=255)                                                       gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGG‑AAAAATTACACACCGTAGCCCCCCCATCAGAAAGCGGGGTGCGCACACCCAGc                                                                                                 >  1:52310/1‑150 (MQ=255)                                                          ggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCg                                                                                             <  2:969192/151‑1 (MQ=255)                                                               aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGc                                                                                         <  1:166693/150‑1 (MQ=255)                                                                  ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatc                                                                                     <  1:457309/151‑1 (MQ=255)                                                                  ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatc                                                                                     <  1:321092/151‑1 (MQ=255)                                                                    gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatcat                                                                                   <  2:8858/151‑1 (MQ=255)                                                                           cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                                                            <  1:720823/151‑1 (MQ=255)                                                                             aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCC‑CGGG‑AAAAATCACACACGGTAGGCCCCCCATCCGAAAACGGGGGGCGCCCACCGAGCTGCGTAGCGATGCTcccct                                                                              >  2:919274/1‑145 (MQ=255)                                                                              aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGG‑AAAAATCACACACGGTAGCCCCCCCATCAGAAAACGTGGTGCGCACACCGCGCTGCGGCACGATCATCCCTTTCt                                                                           >  1:640569/1‑149 (MQ=255)                                                                                  cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATaa                                                                      <  1:328307/150‑1 (MQ=255)                                                                                  cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATaa                                                                      <  1:325031/150‑1 (MQ=255)                                                                                  cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATa                                                                       <  2:900205/149‑1 (MQ=255)                                                                                  cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCC‑CGGG‑AAAATTTACCCAATGTAACCCCCCCATCCGAAAACGGGGTGCGGCCCCCCAGCAGCGTAGAGATTATCCCCTTCTGATaa                                                                      >  1:362633/1‑150 (MQ=255)                                                                                          gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGG‑AAAAATCACACACGGGAGCCCCCACATCAGAAAGCGGGGTGCGCACACCCAGCCGCGTAGCGATCATCCCCTTCTGATAAAACAGCa                                                               >  1:1049233/1‑149 (MQ=255)                                                                                                tGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGG‑AAAAATCACACACGGGAGCCCCCCCCTCAGAAAACGGGGGGGGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACTGGcacc                                                       >  1:15120/1‑151 (MQ=255)                                                                                                tGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGG‑GGCGCG‑CGGG‑AAAAATTTCACACCGGGGGCCCCACATCAGAAAACGCGGGGGGGACAACCACCCGCCTAGCGATTATCCCTTTCTGTTaaaaccggaccgccac                                                        >  2:300107/1‑137 (MQ=255)                                                                                                   aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACt                                                    <  2:222768/151‑1 (MQ=255)                                                                                                       cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGG‑AAAAATCACACACGGTAGCCCCCACATCAGAAAGCGTGGGGGGCACACCGCGCGGCGTAGCGATCATCCCCTGCTGATTAAACACCACTGGccccac                                                     >  1:278164/1‑146 (MQ=255)                                                                                                        ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGG‑AAAAATCACACACGGGAGCCCCCACCTCAGAAAACGTGGGGCGCACACCGAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAACACTGGgaccacct                                                   >  1:688324/1‑145 (MQ=255)                                                                                                            ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGG‑CGGA‑AAAA‑TCAAAAACGGTAACCCCCCCAACAGAAAAACTGGTGCGCACACCCAGCTGCGGAGCGAACATCCCCTTCTGTTAAAACCGCACCGGCACCCCTTCGCggg                                             >  2:985543/1‑149 (MQ=255)                                                                                                              aggaAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAAATCACAAAAGGGAGCCCCCACATCCGAAAACGGGGGGGGCCCACCCAGCTGCGGAGCGATCATCCCCTGCTGATAAAACAAAACCGGCACCACTTCGCGGGGAtg                                         >  1:769615/1‑151 (MQ=255)                                                                                                                 aaagaaaagaGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggtg                                      <  1:771656/151‑1 (MQ=255)                                                                                                                    gaaaagaGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGG‑AAAAATCACACACGGGAGCCCCCACATCAGAAAGAGGGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCCCGGGGAtgg                                        >  2:987948/1‑146 (MQ=255)                                                                                                                       aagaGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGG‑AAAATTCACACACCGCAGCCCCCCCCTCAGAAAGCGTGGGGGGCACACCGAGCTGCGGAGCGATCATCCCCTGCTGATAAACCAGCACTGGCACCACCTCGCGGGGATGGGGGGtgtg                                >  1:665393/1‑151 (MQ=255)                                                                                                                         gagCGACCATTTTGTGGCGGG‑GGCGCGCCGGA‑AAAATT‑AAACACCGGAGCCCCCACAACCGAAAACGGGGGGCGCCCACCCAGCCGCGCGAGGAATATCCCCTGCTGGTAAACCACCACCCGCACCACCTCCCGGGGAtggtggg                                    >  2:704694/1‑144 (MQ=255)                                                                                                                          agCGACCATTTTGTGGCGGT‑GGCGCG‑CGGG‑AAAAATTACAAAAGGGAGCCCCCCCATCAGAAAGCGGGGGGGGCACACCGCGCCGCGGAGCGATCATCCCCTTCTGATAAACCAGCACCGGGCCCCCTTCGCGGGGGTGGGGGCTGTGAcc                             >  2:837142/1‑151 (MQ=255)                                                                                                                                cATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGt                       <  2:479546/151‑1 (MQ=255)                                                                                                                                cATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGt                       <  1:298245/151‑1 (MQ=255)                                                                                                                                         gcggt‑ggcgCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt              >  1:1003306/1‑151 (MQ=255)                                                                                                                                               ggcgCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCCCCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg          <  2:886852/150‑1 (MQ=255)                                                                                                                                                       ggA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGCTGCCGTGATCCGc  >  1:914546/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                      GCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGC  >  NC_000913/3508507‑3508805

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