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.8%	78.3 / 20.2	55	E384D (GAA→GAT)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (9/23);  new base A (23/0);  total (32/23)
Fisher's exact test for biased strand distribution p-value = 1.98e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAC  >  NC_000913/3508514‑3508811                                                                                                                                                     |                                                                                                                                                      aCGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAtt                                                                                                                                                      <  2:492952/150‑1 (MQ=255)     tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATcaca                                                                                                                                                  >  2:102653/1‑150 (MQ=255)               gCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGccc                                                                                                                                       <  2:561437/151‑1 (MQ=255)                  ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGTAGccccc                                                                                                                                     >  1:879908/1‑150 (MQ=255)                     gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACAt                                                                                                                                 <  1:282224/151‑1 (MQ=255)                       gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACCGTAGCCCCcac                                                                                                                                   >  1:1023455/1‑147 (MQ=255)                       gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCa                                                                                                                               <  1:158311/151‑1 (MQ=255)                        aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGCCCCCCCATc                                                                                                                                >  1:102840/1‑149 (MQ=255)                          cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGAAAAAATCACAAACGGGAGACCCCCCATCAGa                                                                                                                             >  1:456730/1‑150 (MQ=255)                          cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATAACAAAACGGAGCCCCCCCCTCAGa                                                                                                                             >  2:955384/1‑150 (MQ=255)                            atatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTAGCCAGCGGATAAAATTACAGACGGTGTGACCCCCATCAGa                                                                                                                             >  2:18781/1‑148 (MQ=255)                                 cAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGGCCCCCCATCAAAAAACgcgg                                                                                                                     >  1:534249/1‑151 (MQ=255)                                    aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTgc                                                                                                                  <  1:27749/151‑1 (MQ=255)                                         cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGca                                                                                                               <  1:271319/149‑1 (MQ=255)                                          aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGGGGGAAAAATCAAAAACGGGAGCCCCCACATCAAAAAACgggggg                                                                                                                   >  2:945350/1‑142 (MQ=255)                                          aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACAAACAGGAGCCCCCACATCAGGAAGCgggggggg                                                                                                                 >  2:295430/1‑142 (MQ=255)                                                 aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCCACATCAGAAAACGTGGGGCGCacac                                                                                                            >  1:773615/1‑144 (MQ=255)                                                      gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGGGGGAAAATTTACAAACGGGAGCCCCCCCCTCAGGAAGCGGGGTGCGCACACCGAGCCGc                                                                                                    >  2:54594/1‑147 (MQ=255)                                                      gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTa                                                                                                 <  2:106627/150‑1 (MQ=255)                                                         aCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGa                                                                                             <  2:368263/151‑1 (MQ=255)                                                          cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGat                                                                                            <  1:674871/151‑1 (MQ=255)                                                              tCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatcat                                                                                         <  2:312161/150‑1 (MQ=255)                                                                 gATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCt                                                                                     <  1:544203/151‑1 (MQ=255)                                                                    cAAAACTCATAGGGTGTCTCCTGCAAACCGGCGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTACACAAGGGATCCCCCACATCACAAAGGGTGGTGCGCCCACCGAGCATCGCTGCGATCATtcc                                                                                      >  1:919458/1‑147 (MQ=255)                                                                       aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGa                                                                               <  1:806859/151‑1 (MQ=255)                                                                        aCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGCCCCCCCCTCCGAAAACGGGGGGCGCACACCCAGCTGCGTAGCGATCCTCCCtt                                                                                    >  2:563312/1‑145 (MQ=255)                                                                           cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATaaa                                                                           <  2:921568/151‑1 (MQ=255)                                                                            aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATCAAACACCGGAGCCCCCCCAACAGAAAACGGGGGGGGGACACCCAGCAGCGTAGCGATCAACCCCTGCTGGTaaa                                                                           >  2:524653/1‑150 (MQ=255)                                                                               gggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATCCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAg                                                                       <  1:801171/151‑1 (MQ=255)                                                                                ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGc                                                                      <  2:943854/151‑1 (MQ=255)                                                                                 gtgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCAAACACCGGAACCCCCACCTCAGGAAGCGGGGGGGGCACCCCCCGCTGCGGAGCGATCATCCCTTGCTgtaaaaaca                                                                        >  2:710280/1‑140 (MQ=255)                                                                                             aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACtt                                                         <  1:216469/151‑1 (MQ=255)                                                                                                 ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACAAACGGGAGCCCCCCCATCAAAAAAAGGGGGGCGCACAACCAGCTGCCTAGCGGTAATCCCCTGCTGGTAAAACAGAACTGGgaccacct                                                         >  2:431287/1‑145 (MQ=255)                                                                                                 ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTACACACCGTAGCCCCCCCATCAGAAAACGTGGGGCGCACACCGAGCTGCGGAGCGGTCATCCCCTTCTGATAAAACAGCACCGGccccacct                                                         >  1:1065279/1‑145 (MQ=255)                                                                                                  gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGGGGGAAAAATTAAAAACGGGAGCCCCCACCACAGAAAGGGGGGGGGGGACAACCAGCTGCGGAGGGGTCATCCCCTGCTGATTAAACAGCACCGGCACCACCTCCCgg                                                    >  2:565270/1‑151 (MQ=255)                                                                                                  gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACAAACGGGAGCCCCCCCCTCAGAAAACGGGGGGGGCACACCCAGCCGCGGTGCGAACATCCCCTGCTTATAAAACAGCACTGGCACCCCCTcgccg                                                    >  2:125065/1‑149 (MQ=255)                                                                                                       aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGTAGCCCCCCCATCAGAAAGCGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTTCTGATAAAACAGCACTGGCACCCCCTCGCGGGGAt                                                >  1:997054/1‑150 (MQ=255)                                                                                                       aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtg                                               <  2:12835/151‑1 (MQ=255)                                                                                                        ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGTAGCCCCCACATCAGAAAGCGGGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTTCTGGTAAAACAACACCGGCCCCACCTCGCggg                                                   >  1:509327/1‑146 (MQ=255)                                                                                                          aaagaaaagaGCGACCATTTTGTGGCGGTGGGGCGGGGGAAAATTTACACACGGTAGCCCCCACATCACAAAGCGCGGGGGGCACACCCGGCTGCGTAGCGATTCCCCCCTGCTGGTAAAACAACACCTGGACCACTTCGCgggggttgtg                                            >  1:745805/1‑144 (MQ=255)                                                                                                          aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACAAACCGGCGCCCCCACATCAAGAAACGGGGGGGGCACACCCACCTTGGGCGCGACCATCCCTCGCTGCTAAAACCACACTGGTCCCACCTCGCGGGGAt                                                >  2:610356/1‑147 (MQ=255)                                                                                                           aagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACGGTAGCCCCCCCATCAGAAAGCGTGGGGCGCACACCGAGCTGCGTAGCGATCATCCCCTGCTGATAAACCAGCACTGGCACCACTTCGCGGGGAtggggg                                           >  1:821634/1‑151 (MQ=255)                                                                                                                 agagCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGTAGCCCCCACCTCAGAAAACGTGGGGCGCACACCGAGCTGCGTAGCGATCATCCCCTGCTGATAAACCAGCACTGGCACCACTTCGCGGGGATGGGGGGtgtg                                      >  1:626911/1‑150 (MQ=255)                                                                                                                 agagCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTAAACACCGGAACCCCCCCAACAGAAAGCGGGGGGCGCACACCCCGCTACGGAACGCTACTCCCCTTCTGATAAAACAACACAGGCACACCCTCGCGGGGATGGGGGGtgtg                                      >  2:881349/1‑150 (MQ=255)                                                                                                                   agCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGGAGCCCCCCCCTCCGAAAGCGGGGTGGGCACACCCAGCTGCGTAGCGATCATCCCTTTCTGATAAACCAGCACTGGCACCACCTCGCGGGGATGGTGGGTGTGAcc                                   >  2:52937/1‑151 (MQ=255)                                                                                                                        ccATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATgg                              <  2:397919/151‑1 (MQ=255)                                                                                                                        ccATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATg                               <  2:124462/150‑1 (MQ=255)                                                                                                                         cATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGt                             <  1:1038686/151‑1 (MQ=255)                                                                                                                            tttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCg                           <  1:431369/150‑1 (MQ=255)                                                                                                                              tgtgGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGa                         <  2:505850/150‑1 (MQ=255)                                                                                                                                   cggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg                   <  2:940066/151‑1 (MQ=255)                                                                                                                                          gcgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGAt            >  2:697475/1‑151 (MQ=255)                                                                                                                                            gcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc          >  1:832342/1‑151 (MQ=255)                                                                                                                                              ggAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc        >  1:1030/1‑151 (MQ=255)                                                                                                                                                   aaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCa   <  2:1007110/151‑1 (MQ=255)                                                                                                                                                    aTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAc  >  1:581173/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                      ACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAC  >  NC_000913/3508514‑3508811

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