Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,270,489	0	A	C	41.0%	45.7 / 17.5	40	intergenic (+887/+136)	yhaC/garK	pentapetide repeats‑related protein/glycerate kinase I
Reads supporting (aligned to +/- strand):  ref base A (9/14);  new base C (16/0);  total (26/14)
Fisher's exact test for biased strand distribution p-value = 8.56e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCG  >  NC_000913/3270340‑3270639                                                                                                                                                       |                                                                                                                                                        ggCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGT‑GGGTTTCCCCCCCAg                                                                                                                                                        <  1:654314/151‑1 (MQ=255) ggCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCAGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGT‑GGGTTTCCCCCCCAg                                                                                                                                                        <  1:93694/151‑1 (MQ=255)      tGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTTTTTCACTTGTCGT‑GGGGTTTCCCCCCCGGCg                                                                                                                                                     >  1:648845/1‑149 (MQ=255)        cTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTTCACTGGTCGT‑GGGGTTTCCCCCCCGGCggt                                                                                                                                                   >  1:427236/1‑147 (MQ=255)          ccGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGCTTGCTCTCTCTTTCACCGGGCGT‑GGGGTTTCCCCCCCCGCGTTAcc                                                                                                                                                >  2:446670/1‑150 (MQ=255)             ggTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTTCACCGGGCGG‑GGGGTTCCCCCCCCCGCGGTAcc                                                                                                                                                >  1:658770/1‑147 (MQ=255)                       aCCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCt                                                                                                                                  <  2:616310/151‑1 (MQ=255)                            gCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGa                                                                                                                             <  1:119369/151‑1 (MQ=255)                            gCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTCTTTGACTCGTCGT‑GGGTTTTCCCCCCCCGCGGTACCTGGCccccccccc                                                                                                                                   >  2:424809/1‑140 (MQ=255)                                     tGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTTCACTCGTCGT‑GGGGTTCCCCCCCCGGCGTTACCTGGCACCCTGCCCCATGGAGccccg                                                                                                                       >  1:499883/1‑146 (MQ=255)                                            aGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCACTGGTCGT‑GGGGTTCCCCCCCCGGCGTTACCTGGCCCCCTGCCCTATGGAGCCCGG‑GCTTtcctcc                                                                                                            >  1:686331/1‑151 (MQ=255)                                            aGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGTTCTCTGTTTCACTGGTCGT‑GGGGTTTCCCCCCAGGCGTTTCCCGGCACCCTGCCCCATTGAGCCCGG‑GCCTtcctcc                                                                                                            >  1:559984/1‑151 (MQ=255)                                                           tctTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTTCACTGGGCGT‑GGGGTTCCCCCCCCGGCGTTACCCGGCACCCCGCCCCATGGAGCCCCG‑GCTTTTCTCCCCCCCCCCCCTCTcc                                                                                             >  1:652866/1‑151 (MQ=255)                                                             ttACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTCTTGCCCTTGGCGT‑GGGGTTTCCCCCCCGGCGTTACCCCGCACCCCCCCCCATGGGGCCCGC‑CCCTTTCTCTCCCCCCCCCCTcccccc                                                                                           >  2:400228/1‑151 (MQ=255)                                                              tACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTGCACTGGTCGT‑GGGGTTTCCCCCCCGGCGGTTCCTGGCCCCCCGCCCCATTGAGCCCCG‑GCTTTTCTCCCCCCCCCCCCTCcccccc                                                                                          >  1:147210/1‑151 (MQ=255)                                                               aCCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCg                                                                                         <  1:341419/151‑1 (MQ=255)                                                                   cACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTTTGTTGCCCTCGTTCT‑GGGGTTTCCCCCCCGGCGGTTCCTGGCCCCCCTCCCTTTTGAGCCCGGCACTTT‑CTCCCCCCCCCCCCCCTcccc                                                                                           >  1:111317/1‑145 (MQ=255)                                                                     cccTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTTCACTGGGCGT‑GGGGTTCCCCCCCCGGCGTTACCTGGCACCCTTCCCTATGGAGCCCGG‑ACTTTTCTCCCCCCCCCCCGTCTCCCCCCGAGAgg                                                                                   >  1:10836/1‑151 (MQ=255)                                                                     cccTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTGCACCGGTCGT‑GGGTTTTCCCCCCCCGCGGTACCTCGCCCCCCGCCCCCTGGGGCCCGG‑ACTTTCCCCCCCCCCCCCCCTCccccccccccagg                                                                                   >  2:1236/1‑144 (MQ=255)                                                                      ccTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGga                                                                                  <  2:385353/151‑1 (MQ=255)                                                                       cTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCCTTTGTACTCGGCTTGGGGTTTCCCCCCCGGGGGTTCCCCGG‑GCCCTTCCCCTTGGGGCCCCG‑ACTTTCCTCCCCCCCGCCCGCCTccccccgagggggc                                                                                 >  2:188544/1‑142 (MQ=255)                                                                           cACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGgacgac                                                                              <  2:575701/150‑1 (MQ=255)                                                                             cccTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGgacgacgac                                                                           <  1:58418/151‑1 (MQ=255)                                                                              ccTTACCGGTTCCCGCTTGCGCGGGCCATCGGCGGTTTGCCCTCTGTTGCGCTGGTCGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGTacgccgacg                                                                          <  2:562952/151‑1 (MQ=255)                                                                                      gATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCACTGGTCGT‑GGGGTTCCCCCCCAGGCGTTAACTGGCACCCCGCCCCATGGAGCCCCG‑GCCTTTCTCCCCCCCCCCCCGCCCCCCCCAAGAGGACGACGACGaa                                                                        >  1:603578/1‑145 (MQ=255)                                                                                      gATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACCGGTCGT‑GGGTTTCCCCCCCAGGCGTTACCTGGTACCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAgcgg                                                                    <  2:466698/149‑1 (MQ=255)                                                                                       aTCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTGCACTGGTCGT‑GGGGTTCCCCCCCCGGCGTTACCCGGCACCCCGCCCTATGGAGCCCCG‑GCTTTCCTCCCCTCCCCCCGTCTCCCCCCAAGAGGACGACGACGAAgcggc                                                                   >  1:226188/1‑149 (MQ=255)                                                                                           cGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTCTTGAACCGGGCGT‑GGGGTTCCCCCCCCCGCGGTACCCCGCACCCCGCCCTATTGGGCCCGG‑CCTTTCCTCCCCCCCGCCCCGCTCCCCCCCAGAGGgcgccgccgcggccgggcccg                                                              >  2:216552/1‑138 (MQ=255)                                                                                            gCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCACTGGGCGG‑GGGGTTCCCCCCCCGGCGTTACCTGGCACCCCGCCCCATTGAGCCCGG‑ACTTTCCTCCCCCCCCCCCCCCTCCCCCCAAGAGGACGACGACGAAGCGGCGACTGt                                                             >  1:553262/1‑150 (MQ=255)                                                                                                               ggTTTGCTCTCTGTTGCACTGGTCGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGgcgcgc                                          >  1:237166/1‑150 (MQ=255)                                                                                                                 tttGCTCTCTGTTGCACTGGTCGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAg                                        >  2:463878/1‑150 (MQ=255)                                                                                                                              gCACTGGTCGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTgc                          >  2:460327/1‑151 (MQ=255)                                                                                                                                   ggTCGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTgcgcg                       <  1:470136/149‑1 (MQ=255)                                                                                                                                    gTCGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGccc                    <  1:188591/151‑1 (MQ=255)                                                                                                                                     tCGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCt                   >  1:481888/1‑151 (MQ=255)                                                                                                                                       gT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTg                 >  2:311932/1‑151 (MQ=255)                                                                                                                                               tCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGG‑CCTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGCTTGCGCGCCCTTGTCAcccc          >  1:651253/1‑151 (MQ=255)                                                                                                                                                cccccccAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCAcccc          <  2:18047/150‑1 (MQ=255)                                                                                                                                                cccccccAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCcg         >  2:58303/1‑151 (MQ=255)                                                                                                                                                     ccAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTg    <  1:159586/151‑1 (MQ=255)                                                                                                                                                       aGGCGTTACCTGGCCCCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTgcg  <  1:65290/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                        GGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGTCGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGG‑ACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCG  >  NC_000913/3270340‑3270639

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