Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,027,2330TA41.4% 26.1 / 11.7 29intergenic (+5/+376)fre/fadANAD(P)H‑flavin reductase/3‑ketoacyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base T (14/3);  new base A (1/11);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 1.16e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCGCTTTCGTTGAAGCATCCTGGTCTGCAAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACC‑GT‑GTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAG‑ATGGCGAAAA‑TTGCCCGCGATCTGTTTTGC‑AGTGAGCGTA‑ATGCGCGGGAAGA‑TCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTA  >  NC_000913/4026985‑4027472
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
ggCGTTTTCGTTGAAGCATCCTGGTCTGCAAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACC‑GT‑GTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAG‑ATGGCGAAAA‑TTGCCCGCGATCTGTTTTGC‑AGTGAGCGTA‑ATGCGCGGGAAGA‑TCGCCTGTTTTGCGATGCGTTTGCATTTATCTGAGATATaa                                                                                                                                                                                                                                               <  2:325919/251‑1 (MQ=255)
      ttCGTTGAAGCATCCTGGTCTGCAAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACC‑GT‑GTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAG‑ATGGCGAAAA‑TTGCCCGCGATCTGTTTTGC‑AGTGAGCGTA‑ATGCGCGGGAAGA‑TCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATaaaaaa                                                                                                                                                                                                                                           <  2:168315/249‑1 (MQ=255)
                    cTGGTCTGCAAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACC‑GT‑GTTAACGGCGGTATTGGAGGGTCACAGGAAGCTGGGAGAGCAAGATATCTTTATTTCCGGAAGGTTTTAG‑ATGGGGGAAA‑ATGCCCGCGGTCTGTTTTGG‑AGTGAGCGTA‑ATGCGCGGGGAGAATC‑CCTGTTTTGCGGTGCGGTTGGATTTTTTTGGGTTATAAAAAAACCCCCCCCCGAccgg                                                                                                                                                                                                                           >  1:404400/1‑248 (MQ=255)
                      ggTCTGCAAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACC‑GT‑GTTAACGGCGGGATTGCACGATCACGGTACGCTGGCCGAGAAAGGAATTTCTATTTCCGGGCGCTTTTTG‑ATGGCGGAAA‑TTGCCCGCCGCCTGTTTTGC‑GGTGAGAGTA‑ATGCGCGGGAAGA‑TCGCCTGTTTTGCGATGTGTTTGGCTTTTTCTAGGATATAAAAAAAACCGCCCCTGACaaggg                                                                                                                                                                                                                         >  1:468306/1‑247 (MQ=255)
                        tCTGCAAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACC‑GT‑GTTAACGGCGGTATTTCAGGAGAACGGTACGCTCGCAGAGCAAGATATCTTTATTTCCCGGCGTTTTTAG‑ATGGCGAAAAATTGCC‑GCCGTCTGTTTTGC‑AGTGAGCGGA‑ATGGGAGGGGAGA‑TAGCCCGTTTTGGGGTGCGCTTGCATTTATCTTAGATATAAAAAAACCCGCCCCCGTCAGGCggg                                                                                                                                                                                                                       >  1:335305/1‑251 (MQ=255)
                                    gcgggtggtTGAACAACGGGAAGCGGGGGGGGGTGGGCGT‑CTGGCACCCGT‑GTTAGCGGCGGTTTTGCAGGATCACGGTACGGTGGCAGAGCATTATATCTATTTTGCCGGCCGTTT‑GAGGAGGGCGAAAA‑TGGCCCGCGTTCGTTTTTGC‑AGTGGGTGTA‑AGGGGGGGAAAGA‑TCGCCTGTTTGGCGTTGCTTTTGTATTTATCGGAAATAAAAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAAc                                                                                                                                                                                                           <  1:320394/248‑1 (MQ=255)
                                      cGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACC‑GT‑GTTAACGGCGGCATTACAGGGAGACGGGACGCCGGCAGAGACTGATATCTTTATTTCCGGACGTTTTTAG‑ATGGCGGAAA‑TTGTCCGCGGTCTGTTTTTT‑AGAGAGGGGA‑ATGCGCGGGAAGA‑AAGCCTGCTTTGCGATGCGTTTGCGTTTTTCTGGGATAAAAAAAAACACGCCCCTGACCGGCGGGAAGAACgagaact                                                                                                                                                                                                          >  1:250837/1‑244 (MQ=255)
                                          ggttgacaaccggggaGCGGGCTGGCGTGGGCGTACGGGC‑CC‑GTTGTTAGCGGCGGTTTGGCAGGAGCACGTTCCGCTGGCAGAGAATGATTTCTTTATTCCCGGACTTTTTGAG‑ATGGCAAAAA‑TTCCCCGGGTTCTTTTTTGC‑AGTGGGCTTA‑AGGCGCGGAAAGA‑TCGCCTTTTTGGCGATGCTTTTGCATTTATCTGAGATAAAAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACt                                                                                                                                                                                                     <  2:314485/236‑1 (MQ=255)
                                                           ggggcggggggggtggtgggggcccc‑GT‑GTTAAGGGGGGTTTTGAAGGATCAGGTTCCGTTGGCAGGGCATAATTTGTATATGGCGGGAGTTTTTGAG‑GGGGCAAAAA‑TTCCCCGCGATCTGTTTTGC‑AGTGGGGATA‑AGGCGGGGGAATA‑CCCCTTGTTGGGGGAGGCGTTTGAATTTATTTGAGATAAAAAAAACCCCGCCCCGGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCAttt                                                                                                                                                                                    <  2:281509/227‑1 (MQ=255)
                                                                     ggggCGTCCGGGCACC‑GT‑TTTAGCGGCGTTATTGCGGGATCACGGTACGCGGGCAGAGCATGATTTCTATATTGCCGGACGTTTTGAG‑ATGGCAAAAT‑TTGCCCGGGTTTTGTTTTGC‑AGTGAGCTTA‑ATGCGGGGGAAGA‑CCGCCTGTTGGGCGTTGCGTTTGCATTTATCTGAGATAAAAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGa                                                                                                                                                                          <  1:153613/250‑1 (MQ=255)
                                                                                              cggcggTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAG‑ATGGCGAAAA‑TTGCCCGCGATCTGTTTTGC‑AGTGAGCGTA‑ATGCGCGGGAAGA‑TCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTTTTATTCAGTGGCATTTAGATTTATGACGTATCTGGCAAACCATGCCCg                                                                                                                                                    >  1:322455/1‑250 (MQ=255)
                                                                                                   gtctttttggATCCCGGTACGCTGGCAAAGCATGATATCTATATTGCCGGACGTTTTGAG‑ATGGCGAAAA‑TTGCCCGCGATCTGTTTTGC‑AGTGAGCGTA‑ATGCGCGGGAAGA‑TCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCgg                                                                                                                                              >  1:365575/9‑250 (MQ=255)
                                                                                                            gATCACGGTACGCTGGCAGAGCATGATATCTATATTGCAGGACGTTTTGAG‑ATGGCGAAAA‑TTGCCCGCGATCTGTTTTGC‑AGTGAGCGTA‑ATGCGCGGGAAGA‑TCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTcgc                                                                                                                                     >  1:413155/1‑251 (MQ=255)
                                                                                                                     aCGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAG‑CTGGCGACAA‑TTGCCCGCGATCTGATTTGC‑AGTGAGAGTA‑ATGCGCGGGAAGA‑TCGACTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAAATGTTATTCAGGGGCATTTAGATTAATGAAGTATCTGGCAAAACATGCCCGATGCGGCGCTGTCGCGGCTTATCg                                                                                                                            >  1:230705/1‑251 (MQ=255)
                                                                                                                     aCGCTGGCAAATCATGATATCTATATGGCGGGCCGTTTTGAG‑GTGGAGAAAA‑TTCCCCGCGATCTTTTTTGC‑AGTGAGCGTA‑AGGCGCGGAAAGA‑CCGCCTGTTCGGCGATGCTTTGGCATTTATCGAAGATAAAAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATc                                                                                                                             <  1:263886/250‑1 (MQ=255)
                                                                                                                          ggCAGAGCATGATATTTATATTGCGGGCCGTTTGGAG‑ATGGCGAAAA‑TTGCCCGGGATCTGTTTTGA‑AGGGAGCGTA‑ATGCGCGGGAAGA‑CCGCCTTTTTGGCGAGGCGTTGGAATTTATCTGAGAAAAAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCt                                                                                                                       <  2:32517/251‑1 (MQ=255)
                                                                                                                                           attttccgggACGTTTGGAG‑AGGGCGAAAA‑TTCCCCGCG‑TTTTTTTTGCGAGTGAGGGTA‑ATGCGGGGAAAGA‑TCCCGTGTTGGCCGATGTGTTGGCTTTTCTTGGAAATAAAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAcc                                                                                                       <  1:356903/242‑1 (MQ=255)
                                                                                                                                                       gTTTTGAG‑ATGGCGAAAA‑TTGCCCGCGATCTGTTTTGC‑AGTGAGCGTA‑ATGCGCGGGAAGA‑TCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAACCCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATa                                                                                           >  1:401067/1‑250 (MQ=255)
                                                                                                                                                              g‑aTGGAGAAAA‑TTGCCCGGGTTCTTTTTTGC‑AGTGAGC‑TACATGCGCGGAAAGA‑TCCCTTTTTTGCCGATGCGTTTCCATTTATCTGAGATAAAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGTATAAGGCGttt                                                                                   <  2:75134/251‑1 (MQ=37)
                                                                                                                                                                aTGGCGAAAA‑TTGCCCGGGTTTTGTTGTGC‑GGTGAGCGTA‑AGGGGCGGGAAGA‑TCTCCTGTTTGGCGGTGCGTTTGCATTTATCTGAGATAAAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTa                                                                                  <  2:335305/251‑1 (MQ=38)
                                                                                                                                                                  ggCGAAAA‑TTGCCCGCGATCTGTTTTGC‑AGTGAGCGTA‑ATGCGCGGGAAGA‑TCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCACGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGACCGAACCGTAGGCCGGGTAAGGCGTTTAcg                                                                                >  2:341503/1‑251 (MQ=255)
                                                                                                                                                                            ggCCCGCGATCTGTTTTGC‑AGTGAGCGTA‑ATGCGCGGAAAGA‑TCGCCTGTTGGGCGATGCGTTTGCATTTATCTGAGATAAAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                       <  1:138845/250‑1 (MQ=255)
                                                                                                                                                                               cccgcattctgtttgggggtgggCGTA‑ATGCGGGGGAATA‑TTGCGTGTTTGGGTATGCGTTGACATTAATCTTGGAAATAAAAAATGCCGCCCCAGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCAGTTAGATCTATGACGGATCGGGCAAACCAAGCCGGAGGCGACGCGGTCGTGGCTTATCGTGCCTTCATATAGTGCGAGCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                     <  1:454612/229‑1 (MQ=255)
                                                                                                                                                                                              c‑AGTGAGCGTA‑ATGCGCGGGAAGA‑TCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGc                                                     >  1:285022/1‑251 (MQ=255)
                                                                                                                                                                                                   gAGCGTA‑ATGCGCGGGAAGA‑TCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATTGGGCCTCCAAAGAGTCCGAACCGTAGGCCGGTGAAGGCGTTTACGCCGCCTACGGCACTTGGTGCATGATGCCTgg                                                 >  1:384433/1‑250 (MQ=255)
                                                                                                                                                                                                     gCGTA‑ATGCGCGGGAAGA‑TCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATa                                                                                                                >  2:249867/1‑186 (MQ=255)
                                                                                                                                                                                                       gTA‑ATGCGGGGGAAGA‑TCGCTTGTTTGGCGAGGCTTTTGCTTTTATCGGAGATAAAAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCg                                             <  1:19801/251‑1 (MQ=25)
                                                                                                                                                                                                                 gggAAGA‑TCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTAcgcc                                                                              >  1:138500/1‑209 (MQ=255)
                                                                                                                                                                                                                                                    tATCTGAGATATAAAAAACCCCGCCCCTGACAGGCGGGAAGCCCGGCAAATAAAATGTTATTCAGTGGCATTTACATCTATGAAGCATCTGGCAAACCATGCCAGCTGCGACGCTGTCGAGGCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTTCGCCGCATCCGGCCCTTGGTGGATGATGCCTGGTGCGACGCTGTCCCGTCTTATTCTGCCTACAAATAGGCCGAACCGtt  >  2:440556/1‑250 (MQ=21)
                                                                                                                                                                                                                                                    atgcgagaaagtaAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAATCGGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGAGGCGACGCTGTCGCGTCTTATCGTGCCTGCAAAAAGTCCGAACCGTAGGCCGTATAAGGGGTTTACGCCGCATTCGGCAATGGGTGCATGATGCCTGATGGTACGCTGTCGCTGCTTATCGTGCCTACAAATAGTCCGACCCGt   <  2:436797/239‑1 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
GGCGCTTTCGTTGAAGCATCCTGGTCTGCAAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACC‑GT‑GTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAG‑ATGGCGAAAA‑TTGCCCGCGATCTGTTTTGC‑AGTGAGCGTA‑ATGCGCGGGAAGA‑TCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTA  >  NC_000913/4026985‑4027472

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