Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,657	0	A	G	46.2%	66.7 / 15.3	52	F386S (TTC→TCC)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (9/19);  new base G (24/0);  total (33/19)
Fisher's exact test for biased strand distribution p-value = 1.40e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCG  >  NC_000913/3508519‑3508804                                                                                                                                           |                                                                                                                                                    gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGGGGGAAAAATtacac                                                                                                                                          >  2:1062127/1‑150 (MQ=255) gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTcacaca                                                                                                                                         <  1:855769/151‑1 (MQ=255)   ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCacacac                                                                                                                                        <  2:177181/150‑1 (MQ=255)         aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACCGTAGcc                                                                                                                                 >  1:856095/1‑151 (MQ=255)         aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACCGTAGcc                                                                                                                                 >  1:785444/1‑151 (MQ=255)                 agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATCACAAACGGGAGccccc                                                                                                                              >  1:1043852/1‑146 (MQ=255)                 agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGGGGGGGGGAAAAATTACAAACGGGAGccccccc                                                                                                                            >  1:992830/1‑146 (MQ=255)                     cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGGAGccccccc                                                                                                                            >  1:559273/1‑142 (MQ=255)                      aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATTACACACCGTAGCCCCCCCCCCAGaa                                                                                                                     >  2:583457/1‑150 (MQ=255)                      aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGaaa                                                                                                                    <  2:190607/151‑1 (MQ=255)                      aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGaaa                                                                                                                    <  2:1018822/151‑1 (MQ=255)                      aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGaaa                                                                                                                    <  2:907921/151‑1 (MQ=255)                       atatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCAAACACCGGAGCCCCCCCATCAGaaa                                                                                                                    >  2:52802/1‑150 (MQ=255)                         atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCg                                                                                                                 <  1:264504/151‑1 (MQ=255)                          tGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTAAACACCGGAGacccccccaccaaaa                                                                                                                    >  2:43571/1‑137 (MQ=255)                           gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtg                                                                                                               <  2:503353/151‑1 (MQ=255)                           gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtg                                                                                                               <  2:381530/151‑1 (MQ=255)                              taTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGTAGCCCCCCCACCAGAAAGCgtgggg                                                                                                            >  2:813980/2‑149 (MQ=255)                                            aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGGGGGCGGGAAAAATTACAAACGGTAGCCCCCCCATCAGGAAACGGGGGGCGGACACCCAGCTg                                                                                              >  1:688225/1‑151 (MQ=255)                                                 gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTa                                                                                          <  2:50974/150‑1 (MQ=255)                                                  tAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGACAGAAAAGAGCGACCATTTTGTGGCGGCGGGGCGCGGAAAAATTAACCCAACGTCGCCCCCCCCTCCGGAAGCGGGGGGGGGACAcctcgccgcggag                                                                                         >  2:300908/1‑139 (MQ=255)                                                         tCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAGCGGGGGGGGGACACCCAGCTGCGTAGCGgtc                                                                                    >  2:120514/1‑148 (MQ=255)                                                          cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGGAAATTCACACAAGGGAGCCCCCACCTCCGAAAGAGTGGGGCGGACACCGAGACGCGTAgagct                                                                                     >  1:511844/1‑142 (MQ=255)                                                          cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCACACACCGAGCTGCGTAGCGatcatc                                                                                 <  1:52803/150‑1 (MQ=255)                                                             aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAACCCCCCCATCAGAAAACGTGGGGCGCCCACCCAGCCGCGTAGCGATCCTccc                                                                               >  2:420635/1‑149 (MQ=255)                                                                aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACGGTAGCCCCCCCATCAGAAAACGTGGGGCGCACACCCAGCCGCGTAGCGATCATccc                                                                               >  2:1087995/1‑146 (MQ=255)                                                                  aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACAAACGGGAGCCCCCCCATCAGAAAACGGGGGGCGCCCACCCAGCTGCGGAGCGAACATccc                                                                               >  2:57806/1‑144 (MQ=255)                                                                         aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                                                                 <  1:404201/151‑1 (MQ=255)                                                                           ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGc                                                               <  1:330805/151‑1 (MQ=255)                                                                           ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGc                                                               <  1:991105/151‑1 (MQ=255)                                                                              gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCCAGCTGCGGAGCGATCATCCCCTTCTGATAAACCa                                                                 >  2:1055426/1‑146 (MQ=255)                                                                                      cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACCGTAGCCCCCCCCTCAGAAAACGGGGGGCGCACACCCAGCTGCGTAACGATCATCCCCTTCTGATAAAACAACACTGGcacccc                                                    >  1:769384/1‑149 (MQ=255)                                                                                         aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                 <  1:266701/151‑1 (MQ=255)                                                                                          ccGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcg                                                <  2:974218/151‑1 (MQ=255)                                                                                            ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACAAACGGGAGCCCCCCCATCAGAAAACGGGGGGGGCCCACCCCGCTGCGGAGCGGTCATCCCCTTCTGATAAACCAGCACTGGCACCACCTccc                                               >  2:318102/1‑148 (MQ=255)                                                                                                  aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCCACATCAGAAAGCGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTTCTGATAAAACAACACCGGcaccccctccccgggatg                                        >  1:326453/1‑137 (MQ=255)                                                                                                  aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGTAGCCCCCCCCTCAGAAAACGGGGGGCGCACACCCAGCTGCGGAGCGATCATCCCCTTCTGCTAAACCAGCACTGGCACCACCTCGCgggg                                           >  1:798723/1‑147 (MQ=255)                                                                                                     aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggtg                                     <  1:477013/151‑1 (MQ=255)                                                                                                      aagaaaagaGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATTACCAAAGGTAGCCCCCACCTCAACAAGCGGGGGGCGCCCACCCAGCTGCGGCGGGATCAACCCTTCGTGTTAAAACAACACCGACACCACCTCCCGGGGGtggtgg                                    >  1:194959/1‑151 (MQ=255)                                                                                                           aagaGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGGAGCCCCCCCCCCAGAAAAAGGGGGGGGCACCCCCAACCGCGTAGAGATCAAACCTTGCTGATAAAACCGCACTGGGCCCCCCTCGCGGGGAtggggggt                                  >  1:1087815/1‑146 (MQ=255)                                                                                                             gagCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCACATCAGAAAACGGGGGGGGCACACCCAGCTGCGGAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACCTCGCGGGGATGGTGGCtgtg                               >  1:852959/1‑149 (MQ=255)                                                                                                                 gACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACCGTAGCCCCCCCATCAGAAAAAGGGGGGCGCCCACCCAGCTGCGTAGCGATCATCCCCTTCTGATAAAACAACACTGGCCCCACCTCGCGGGGATGGGGGGTGTGAcc                            >  1:666091/1‑148 (MQ=255)                                                                                                                 gACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAAt                         <  2:786706/151‑1 (MQ=255)                                                                                                                       tttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCg                    <  2:59036/150‑1 (MQ=255)                                                                                                                            ggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGt              >  1:843151/1‑151 (MQ=255)                                                                                                                             gcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt             >  1:529491/1‑151 (MQ=255)                                                                                                                               ggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGc           >  1:522725/1‑151 (MQ=255)                                                                                                                                 tggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg         >  1:955979/1‑151 (MQ=255)                                                                                                                                      cgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGa      >  2:241597/1‑149 (MQ=255)                                                                                                                                      cgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc    >  2:1049255/1‑151 (MQ=255)                                                                                                                                        cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg  >  2:753110/1‑151 (MQ=255)                                                                                                                                        cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg  <  1:228264/151‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                    GTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCG  >  NC_000913/3508519‑3508804

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