Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,662	0	T	A	41.7%	72.0 / 15.4	49	E384D (GAA→GAT)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (7/21);  new base A (20/0);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 8.13e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGC  >  NC_000913/3508515‑3508805                                                                                                                                                       |                                                                                                                                                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)                                                                                                                                                ggttataTTCATACACGTTAGCCCCCACATCAGAAAGCGCAGCGCGCACACCGAGCTGCGTAGCGAACATCCCTTGCCGATAAACCAGCACTGGCACCACTTCGCGCTGATGGATGCGGTGACCAATGGTCGGATCTTTGCCGTGa       >  2:972752/7‑146 (MQ=255)                                                                                                                                                       |                                                                                                                                                CGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGT‑GGCGCG‑CGGA‑AAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGC  >  NC_000913/3508515‑3508805

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