Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,029,959	0	A	C	51.2%	59.8 / 79.9	44	C325G (TGT→GGT)	ygfT	putative oxidoreductase, Fe‑S subunit/nucleotide‑binding subunit
Reads supporting (aligned to +/- strand):  ref base A (18/3);  new base C (0/22);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 2.19e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-10

CCCATTGCGGTGAATATCTCTCGCCGCTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGCCCAG‑CGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGACTGCGCCCGAG  >  NC_000913/3029811‑3030101                                                                                                                                                     |                                                                                                                                                cccATTGCGGTGAATATCTCTCGCCGCTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACAcc                                                                                                                                               >  2:508118/1‑151 (MQ=255)       gCGGTGAATATCTCTCGCCGCTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAAcc                                                                                                                                         >  2:740783/1‑151 (MQ=255)         ggTGAATATCTCTCGCCGCTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCt                                                                                                                                       >  2:675135/1‑151 (MQ=255)            gAATATCTCTCGCCGCTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCa                                                                                                                                    <  2:267161/151‑1 (MQ=255)              atatCTCTCGCCGCTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAgg                                                                                                                                  <  1:883230/151‑1 (MQ=255)               tatCTCTCGCCGCTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGc                                                                                                                                 >  1:330269/1‑151 (MQ=255)               tatCTCTCGCCGCTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGc                                                                                                                                 >  1:588193/1‑151 (MQ=255)                atCTCTCGCCGCTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGcc                                                                                                                                >  2:250934/1‑151 (MQ=255)                  ctctcGCCGCTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGccc                                                                                                                               >  1:74993/1‑150 (MQ=255)                         cgcTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCGAACCCTGCAGGCCCAG‑CGCCaa                                                                                                                      <  2:632342/151‑1 (MQ=255)                          gcTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGCCCAG‑CGCCAAt                                                                                                                     >  1:763424/1‑151 (MQ=255)                                  aaTACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGCCCAG‑CGCCAATCACCGcc                                                                                                              >  1:891603/1‑150 (MQ=255)                                          tttATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGCCCAG‑CGCCAATCACCGCCACTTTTTCg                                                                                                     >  2:517973/1‑151 (MQ=255)                                           ttATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGACCTCCTGCGCGCGCAAGAATATCAGCACCCCCTACCCCTGCAGGCCCAG‑CGCCAATCACCGCCACTTTTTCg                                                                                                     <  1:989592/150‑1 (MQ=255)                                            tATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGCCCAG‑CGCCAATCACCGCCACTTTTTCGCt                                                                                                   >  1:937122/1‑151 (MQ=255)                                                  gTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGCCCAG‑CGCCAATCACCGCCACTTTTTCGCTACgggg                                                                                             >  1:405934/1‑151 (MQ=255)                                                    ttaaaaGGAGGAATCCAAAAAGCCACCAACCCACCATTTTCAGAATGGCGACAAAAGCCAGCGCCCTAACCTCCGGCGCGCGCCACAATATCAGCTCGCCATAACCCGGCAGGCCAAG‑CGCAAATCGCCGCCACTTTTTCGCTACggg                                                                                              <  2:87151/145‑1 (MQ=255)                                                           gaggaATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGCCCAG‑CGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACttt                                                                                    >  2:968567/1‑151 (MQ=255)                                                                    cAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGCCCAG‑CGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATc                                                                           >  1:251240/1‑151 (MQ=255)                                                                      aaaGTCAGCATCCCGCCATTTTCTGGATGGCGATAAAAGACATCGACCTGAACTCCTGCGCGCGCAAGAATATCAGCACCCCCTACCCCGGCAGGCCCAG‑CGCAAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCgg                                                                         <  1:77309/151‑1 (MQ=255)                                                                           cAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCCGCAGGCCCAG‑CGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGc                                                                    >  1:242467/1‑151 (MQ=255)                                                                                     ccAATTTCTGGAGGGCGATAAAAGACATCGACCTGAACTCCTGCGCGCGCCAAAATATCAGCACCCCCTACCCCGGCAGGCCCAG‑CGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATc                                                           <  1:602763/150‑1 (MQ=255)                                                                                       aaTTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGCCCAG‑CGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGcc                                                        >  1:557535/1‑151 (MQ=255)                                                                                       aaTTTCGGGATGGCGATCAAACCCATCCCCCAGAACCCCGGCGCGCGCGAAAATACCGGCCCCCCCTACCCCCGCGGCCCCAG‑CCCCAACCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGcc                                                        <  2:975429/151‑1 (MQ=255)                                                                                                    cGATCAAAGACACCGCCCTGACCTCCTGCGCGCGCCAAAATATCACCACCCCCTACCCCTGCAGCCCCAG‑CGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCg                                            <  1:86787/150‑1 (MQ=255)                                                                                                      accAAACACACCGCCCTGACCTCCGGCGCGCACAAGAATTCCAACCCCCCCAACCCCAGCAGGCCCAG‑CGCAAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGt                                          <  1:22374/148‑1 (MQ=255)                                                                                                      aTCAAAGACATCGACCTGACCTCCTGCGCGCGCAAGAATACCAGCACCCCCTACCCCTGCAGCCCCAG‑CGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCgg                                           <  1:371096/149‑1 (MQ=255)                                                                                                      aTCAAAGACACCGACCTGACCCCCTGCGCGCGCCAGAATACCAGCACCCCCTAACCCTGCAGGCCCAG‑CGCCAATCCCCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGt                                          <  2:333356/150‑1 (MQ=255)                                                                                                          aaGACATCTACCTGACCCCCGCCGCGCCCAAAAATACCAGCCCCCCCCACCCGGCCAGCCCCGG‑CCCCAACCGCCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGt                                     <  2:959714/151‑1 (MQ=255)                                                                                                                    cctaacccccTGCGCGCGCAAGAATCTCACCCCCCCCACCCCCGGCGGGCCCAG‑CGCCAATCGCCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTcc                             <  2:444447/141‑1 (MQ=255)                                                                                                                    ccgaacccccggcgcgcccAAGAATATCACCACCCCCTACCCCTGCAGGCCCGG‑CCCCAACCCCCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTcc                             <  2:655954/138‑1 (MQ=255)                                                                                                                     cTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGCCCAG‑CGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTcc                             >  1:607518/1‑148 (MQ=255)                                                                                                                        acccccTGCGCCCGCCAGAATATCCGCACCCCCTACCCCTGCGGGCCCAC‑CCCCAACCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGAtt                        <  1:918430/146‑1 (MQ=255)                                                                                                                        acccccTGCGCCCCCCAAAATACCACCACCCCCTACCCCGCCAGCCCCAG‑CCCCAACCCCCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGAt                         <  2:959106/145‑1 (MQ=255)                                                                                                                        accTCCGGCGCGCGCAAGAATATCAGCCCCCCCTACCCCGGCAGCCCCAG‑CCCAAACCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGAt                         <  2:135654/147‑1 (MQ=255)                                                                                                                            ccTGCGCGCGCAAGAATATCAGCACCCCCAACCCCTGCAGGCCCAG‑CGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTAcc                     <  1:688841/149‑1 (MQ=255)                                                                                                                                  cgcgcCAGAATATCAGCACACCCTAACCCTGCAGGCCCAG‑CGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATag                >  1:854926/1‑148 (MQ=255)                                                                                                                                  cgcgcAAAAATACCACCACCCCCAACCCCTGCAGCCCCGG‑CGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAgaga             <  2:74989/151‑1 (MQ=255)                                                                                                                                    cgctacaaactccgCACCCCCTCCCCCTGCA‑GCCCAGTCGCAATTCTCCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGAc            <  1:470612/137‑1 (MQ=255)                                                                                                                                    cgcAAGAATATCAGCCCCCCCTACCCCGGCAGGCCCAG‑CGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATaga               <  2:473341/147‑1 (MQ=255)                                                                                                                                        aaaaaaccaccaccCCCAACCCCTGAAGCCCCAG‑CCCAAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATa                 <  2:891516/131‑1 (MQ=255)                                                                                                                                        aaaaTATCAGCACCCCCCACCCCTGCAGGCCCAG‑CGCAATTCGCCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGACTgcg        <  2:663696/148‑1 (MQ=255)                                                                                                                                             atCAGCACCCCCTACCCCGGCAGCCCAAG‑CCCCAACCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGACTGCGCCCGa   <  2:733387/150‑1 (MQ=255)                                                                                                                                             accAGCACCCCCTACCCCGGCAGCCCCAG‑CGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGACTGCGCCCGAg  <  2:370604/149‑1 (MQ=255)                                                                                                                                               cAGCACCCCCTAACCCTGCAGGCCCAG‑CGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGACTGCGCCCGAg  <  1:567605/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                CCCATTGCGGTGAATATCTCTCGCCGCTGGCTTAATACCGTTTTATCGAGTTTGAAAGGAGGAATGCCAAAAGTCAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGCCCAG‑CGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGCGCGGTATCGGTGATGTAGCGTTCCAGATTACCGATAGAGACTGCGCCCGAG  >  NC_000913/3029811‑3030101

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