Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,049,9100CG23.8% 70.5 / 19.5 42intergenic (+43/‑58)quiA/oprB‑IIIquinate dehydrogenase/carbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base C (13/19);  new base G (8/2);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 6.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGC  >  NC_002947/4049762‑4050054
                                                                                                                                                    |                                                                                                                                                 
cTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGc                                                                                                                                                   >  5:282996/1‑149 (MQ=255)
       tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCGACCAGGccccgg                                                                                                                                              >  2:304346/1‑145 (MQ=255)
          ccGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGa                                                                                                                                         <  5:166729/149‑1 (MQ=255)
            gAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGc                                                                                                                                       <  7:230220/149‑1 (MQ=255)
                    gCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCggg                                                                                                                                    >  5:347619/1‑144 (MQ=255)
                    gCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCggg                                                                                                                                    <  6:347619/144‑1 (MQ=255)
                    gCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGc                                                                                                                               <  3:159742/149‑1 (MQ=255)
                        gCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGccc                                                                                                                                                 >  6:151678/1‑127 (MQ=255)
                                  gtgcgcACTGCTGGCGGCGCCCGCCAGTCGACTCACCGCGGCGACTATGTGATGTCTTACGCGCTGCCGTAAATGCGCTTCGCCTGCTGCCCGGGCAAGCCCGCTCCTACCAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCt                                                                                                                 <  1:316733/144‑1 (MQ=37)
                                    gATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTATGAGc                                                                                                                                       >  4:72302/1‑125 (MQ=255)
                                       cACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCa                                                                                                               <  7:297996/146‑1 (MQ=255)
                                       cACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCa                                                                                                               >  8:297996/1‑146 (MQ=255)
                                        aCTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGa                                                                                                           <  6:282996/149‑1 (MQ=255)
                                        aCTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGa                                                                                                           >  7:102053/1‑149 (MQ=255)
                                         ctgctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCAACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGAc                                                                                                          <  3:359987/149‑1 (MQ=255)
                                          tgctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTGCTACAAGGCCCCTGTAGGGGCGGGCTTGCCCGGGAATGGCCCTCAAAGAc                                                                                                          >  7:167195/1‑148 (MQ=255)
                                          tgctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCCCGGCGACAAGGGCCCGGTAGGGGCGGGGGTGCCCCCGGAAGGCCCGCAAAGACa                                                                                                         >  8:58687/1‑149 (MQ=255)
                                                  gcgcCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCgaga                                                                                                 <  3:129389/149‑1 (MQ=255)
                                                      ctttcAGTCGACTGACCGTCGCGACTATGTGATTTCTTACCCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCt                                                                                             <  1:338136/145‑1 (MQ=255)
                                                       cgccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTg                                                                                            <  5:357363/149‑1 (MQ=255)
                                                                    aCCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCg                                                                                >  8:213105/1‑148 (MQ=255)
                                                                            gACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATg                                                                                       <  5:320077/133‑1 (MQ=255)
                                                                            gACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCACAAAGAAAACCCGACATTCTGTAAtt                                                                                       >  6:320077/1‑132 (MQ=255)
                                                                                        tCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTGCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCAcc                                                           >  3:276980/1‑148 (MQ=255)
                                                                                          ttACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCAGTAGGAGCGGGCTTGCCCGCGAATGCCCCAAAAAGCAAACCCGACAATCTGTAATGCCATCCGCAATTCGCCTCA‑CCCCCACCctg                                                        >  4:291287/1‑149 (MQ=38)
                                                                                                 ctGCCGTAAAACCGCGACGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGctgg                                                 <  5:233147/149‑1 (MQ=255)
                                                                                                      gTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGctggccag                                             <  3:385225/148‑3 (MQ=255)
                                                                                                                         ttCGCGGGCAAGCCCGCTCCAACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAAACCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCCAGGCCAGCACCACTGCCCTTg                         >  4:250226/1‑149 (MQ=255)
                                                                                                                           cgcgGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc                        <  8:167195/148‑1 (MQ=255)
                                                                                                                           cgcgGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc                        <  7:213105/148‑1 (MQ=255)
                                                                                                                                    gCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACACCCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAt              >  7:62476/1‑149 (MQ=255)
                                                                                                                                     cccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCCTGACCAGCACCACTGCCCTTGCCAGCGACGAtc             <  7:58687/149‑1 (MQ=255)
                                                                                                                                      ccGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct            >  2:321844/1‑149 (MQ=255)
                                                                                                                                        gCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa          >  6:86936/1‑149 (MQ=255)
                                                                                                                                        gCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa          >  3:293403/1‑149 (MQ=255)
                                                                                                                                        gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCATACCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc           <  1:304346/149‑1 (MQ=255)
                                                                                                                                            cTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAc      >  7:151998/1‑149 (MQ=255)
                                                                                                                                            cTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCACACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAc      >  7:341843/1‑149 (MQ=255)
                                                                                                                                            cTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAc      >  4:84290/1‑149 (MQ=255)
                                                                                                                                              ccAAGGCCCCTGTCGGAGCGGGCTTGCCCGCTAATGCCCCTCAACCAAAACCCGACATTCTGAAATGCCATCCGCAATTCGCATCC‑CCCCCACCCTGCGGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgacga               <  2:369325/137‑1 (MQ=255)
                                                                                                                                              aCAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgacga               >  1:369325/1‑138 (MQ=255)
                                                                                                                                              aCAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATg        <  4:117816/145‑1 (MQ=255)
                                                                                                                                              aCAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATg        >  3:117816/1‑145 (MQ=255)
                                                                                                                                                acggCCCCTGTAGGTGCGGGGTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgc  <  3:10803/147‑1 (MQ=255)
                                                                                                                                                aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCg   >  1:283267/1‑148 (MQ=255)
                                                                                                                                                    ccccTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCAGGGCCAGcagtag                                 >  6:369230/1‑110 (MQ=255)
                                                                                                                                                    ccccTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCAGGGCCAGcagtag                                 <  5:369230/114‑5 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCA‑CCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGC  >  NC_002947/4049762‑4050054

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