Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	535,920	0	C	T	59.5%	87.4 / ‑2.7	42	R71C (CGC→TGC)	hyi	hydroxypyruvate isomerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (4/21);  minor base G (17/0);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 4.70e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACG  >  NC_000913/535780‑536067                                                                                                                                             |                                                                                                                                                     cTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGGGCGGGAtt                                                                                                                                              >  2:708987/1‑149 (MQ=255)        gtgtGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTAtt                                                                                                                                     <  2:295124/151‑1 (MQ=255)                   gcggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGt                                                                                                                                        <  2:321878/137‑1 (MQ=255)                   gcggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATGGGGGGGGGGGGGGGGCGGTGATGCCTTTAtatcgggagagcg                                                                                                                          >  1:321845/1‑139 (MQ=255)                            aaTTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGGGGGGGAGGGGGGGATTTCCTGTTTTTCTTGGCGGggag‑aaagggt                                                                                                                >  1:315215/1‑146 (MQ=255)                             aTTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGGGCGGTGTTGCCTGTATTTCTGGCCGGGTAG‑AAGAGttt                                                                                                               >  2:734003/1‑151 (MQ=255)                                    tttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTGGGGGGCGGGGGGGGGGGGGATTGCCTGTTTTCCTGGCCGTGGAG‑AAGGGGTTTggg                                                                                                           >  1:650911/1‑148 (MQ=255)                                       ccTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGGGGGGGGGGGCGGTATTGTCTGTATTTCTTGCCGCGAAG‑GAAGGGTTCGGGGTgg                                                                                                       >  2:312440/1‑149 (MQ=255)                                        cTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGGGGCGGGGGGGGGCGGGATTGCCTGTATTCCTGGGCGCGGAG‑AAGAGATTTggggggggg                                                                                                     >  1:566665/1‑144 (MQ=255)                                          tATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGAGGACGGGCGGGGGGGGGGGGCGGGAATGGCCGGTTTCCCGGCCCGGTAGAAAGAGTTTTggggggtg                                                                                                      >  2:260264/1‑143 (MQ=255)                                             gACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTagca                                                                                                <  2:645357/150‑1 (MQ=255)                                                 aCGACATTTAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTagcagca                                                                                             <  2:542686/149‑1 (MQ=255)                                                   gACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGCGGGGGGGGGCGGTATTGCCTGTATTTCTGGCCGTGGAG‑AAGAGGTTTGGGGTgggggtg                                                                                                  >  2:511918/1‑137 (MQ=255)                                                      aTTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGGGGGGGGGGGCGGGATTGCCTGTAATTCTGGCCGTGGAG‑AAGAGTTTTGGGGGGGCGTagcagcagc                                                                                           >  2:223447/1‑146 (MQ=255)                                                         gaagaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGt                                                                                    <  1:620303/150‑1 (MQ=255)                                                          aagaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTa                                                                                  <  2:12602/151‑1 (MQ=255)                                                              aTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGGGCGGGATTGCCTGTATTCCTGGCCGGGAAG‑AAAAGTTTTGGGGTGGGGTGGCAGGAGGGATTTGTTATGcc                                                                              >  2:517000/1‑151 (MQ=255)                                                                 aaaaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCg                                                                            <  1:593566/150‑1 (MQ=255)                                                                  aaaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCCGGGGGGCGGGGGGGGGCGGGATTTCCTGTGTTTCCCGCCGTGGAG‑AAGAGATTTGGGGTGGCGTAGCAGCAACGATTTGTTTTGCCCgcg                                                                          >  2:557390/1‑149 (MQ=255)                                                                   aaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGAGGGCGGGATTGCCTGTATTTCTGGCCGTGAAG‑GAGAGGTTTGGGATGGGGTAGCAGCCGCGATTTGTTATGGCCCTg                                                                          >  2:176804/1‑150 (MQ=255)                                                                    aCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGGGGGGGGGGGGGGGATTGCCTGTATTTCTGGCCGTGAAG‑AAGAGGTTTGGGGTGGGGTAGCAGCAGCGGTTTGTTATGCCCGTg                                                                          >  2:179161/1‑149 (MQ=255)                                                                     cATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTgcgc                                                                       <  1:423246/151‑1 (MQ=255)                                                                          gCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTggg                                                                   <  2:771164/150‑1 (MQ=255)                                                                          gCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGt                                                                  <  1:365256/151‑1 (MQ=255)                                                                             ggCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGta                                                                 <  1:491399/149‑1 (MQ=255)                                                                             ggCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGta                                                                 <  1:415230/149‑1 (MQ=255)                                                                                     ataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAAt                                                       <  2:356536/151‑1 (MQ=255)                                                                                        aaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTaa                                                    <  1:700963/151‑1 (MQ=255)                                                                                          cTCGAACATACGCTGCACAATTTACCGGCGGGGGACTTGGCGGCGGGGGGGGGCGGGATTGACTGTATTTCTGGCCGGGAAG‑AAGAGGTTTGGGGTGGGGGAGGAGCAGCGATTCGTTATGCCCGTGCGCTGGGGGATAAAAAAATTAACt                                                  >  2:26174/1‑151 (MQ=255)                                                                                                   aCGCTGCACAATTTACCGGCGGGTGAATGGGGGGGGGGGGGGGGCGGGGTTGCCTTTATTCCTGGCCGGGAAG‑AAGAGGTTCGGGAGGGGGGAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGGAATAAAAAAATTAACTGTCTGGt                                           >  2:180400/1‑149 (MQ=255)                                                                                                    cGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGTGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTc                                          <  2:570491/149‑1 (MQ=255)                                                                                                        gCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTat                                     <  1:435114/150‑2 (MQ=255)                                                                                                        gCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTa                                      <  1:101325/149‑1 (MQ=255)                                                                                                         cacaATTTACCGGCGGGTGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTaaaa                                   <  2:551093/151‑1 (MQ=255)                                                                                                               ttACCGGCGGGTGACTGGGGGGGGGGGGGGGGCGGGATTTCCTTTATTTCTGGCCGTGGAG‑AAGAGGTTCGGGGGGGGGGAGCAGCAGCGATTCGTTTTGGCCCGGGGCTGGGGAATAAAAAAATTAACTGTCTTGGCGGTtaaaaacccg                             >  1:730909/1‑145 (MQ=255)                                                                                                                          tGACTGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTcagc                  >  2:654243/1‑151 (MQ=255)                                                                                                                              tGGGCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTg              >  1:300489/1‑151 (MQ=255)                                                                                                                               gggCGGCGGGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGa             <  1:521810/151‑1 (MQ=255)                                                                                                                                     cggGGGAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAAc           >  1:180921/1‑147 (MQ=255)                                                                                                                                       ggggAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTc     <  1:516948/151‑1 (MQ=255)                                                                                                                                          gAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACg  <  2:650969/151‑1 (MQ=255)                                                                                                                                          gAGTGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACg  >  1:202156/1‑151 (MQ=255)                                                                                                                                             |                                                                                                                                                     CTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACG  >  NC_000913/535780‑536067

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