Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,005,0000GC23.4% 81.4 / 27.9 47intergenic (+38/‑138)PP_5259/ydcJLysR family transcriptional regulator/metalloprotein
Reads supporting (aligned to +/- strand):  ref base G (16/20);  new base C (7/4);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.33e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGC  >  NC_002947/6004853‑6005132
                                                                                                                                                   |                                                                                                                                     
aTGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCGTCTTCGCGGCTAAAGCCGCTCCTACGGGGACcgc                                                                                                                                      >  3:47362/1‑149 (MQ=255)
 tGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCcccg                                                                                                                                     <  7:273498/149‑1 (MQ=255)
   cAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTa                                                                                                                                   <  7:75444/149‑1 (MQ=255)
           aGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGc                                                                                                                           <  6:197206/149‑1 (MQ=255)
                 ggAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGc                                                                                                                     >  2:245317/1‑149 (MQ=255)
                    gcgcATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCgg                                                                                                                            <  3:200298/139‑1 (MQ=37)
                    gcgcATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCgg                                                                                                                            >  4:200298/1‑139 (MQ=37)
                      gcATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCt                                                                                                                          <  5:74250/139‑1 (MQ=255)
                      gcATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCt                                                                                                                          >  6:74250/1‑139 (MQ=255)
                        aTTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCCCTACAACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGaa                                                                                                               >  6:222986/1‑148 (MQ=255)
                           aCCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAgg                                                                                                            <  4:154250/148‑1 (MQ=255)
                                   tATGCCGAGCATGCGGCAGAGGTGCCATAGGTGCGGGCGCTGGTGTAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCTCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAg                                                                                                   <  2:337535/149‑1 (MQ=255)
                                     tGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGt                                                                                                 <  1:271595/149‑1 (MQ=255)
                                       ccGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTc                                                                                                >  1:286917/1‑148 (MQ=255)
                                        cGAGCATGCGGCAGAGGGGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGTAAAGGGCCGTCGAGGTCGa                                                                                              <  8:193257/149‑1 (MQ=255)
                                             aTGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCaa                                                                                          >  4:121320/1‑148 (MQ=255)
                                                cGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCaaaaaa                                                                                      <  1:245317/149‑1 (MQ=255)
                                                         tGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGCCGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGaca                                                                            <  8:312994/149‑1 (MQ=255)
                                                           ccAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGgagc                                                                                                                              <  8:270121/98‑1 (MQ=25)
                                                           ccAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGgagc                                                                                                                              >  7:270121/1‑98 (MQ=25)
                                                           ccAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGc                                                                                                                     <  4:57114/107‑1 (MQ=18)
                                                                    cGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACTGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACtttt                                                                 <  5:222986/149‑1 (MQ=255)
                                                                      ggCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCGGGCGAGGTCGAGCAAAAAAGGA‑ATGACACACGCAGtttttt                                                               >  6:362816/1‑149 (MQ=255)
                                                                                gAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCTGCTAAAGCCGCTCCTACGGGGACCGCTTATGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGa                                                     <  5:226689/149‑1 (MQ=255)
                                                                                  aTGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGACATGACACACTCACTTTTTTTCGCTTGTGAg                                                    <  1:213900/149‑1 (MQ=255)
                                                                                            gAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCAtt                                         <  3:121320/149‑1 (MQ=255)
                                                                                               cAGTTGCAGGTTTGACGGCCTCTTCGCGGCTTAAGCCGCACCAACGGGGACCGCGTAGGCGCGGAAATATCCCGGAAACGGCCCTACAGGTAGATCAAAAAAAGA‑AATACACACGCAAATTTTTTACCTTGTTATCAAAGTTCATAcc                                       >  2:67405/1‑148 (MQ=255)
                                                                                                  ttGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGGTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGGCACAGTGACTTTTTTTCGCTTGTGAGCAAAGTGGATTCCCAGc                                   >  5:294238/1‑149 (MQ=255)
                                                                                                  ttGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGc                                   >  5:289279/1‑149 (MQ=255)
                                                                                                    gCAGGTTTGTCGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGGAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCt                                  <  2:286917/148‑1 (MQ=255)
                                                                                                        gTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTTGCCGGGTAAGGGCCGTCGAGGTCGATGAAAAAAGGT‑ATGACACACTCACTTTTTTTCGCTTGTGAGGAAAGTGGATTCCCAGCTTTGAt                             >  5:5132/1‑149 (MQ=255)
                                                                                                        gTTTGACGGCCTCTTCGCGGCTAAAGCCGATCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCAt                             <  1:67405/149‑1 (MQ=255)
                                                                                                             aCGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAg                        >  8:237291/1‑149 (MQ=255)
                                                                                                                 ccTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTc                    >  5:156988/1‑149 (MQ=255)
                                                                                                                     ttCGCGGCTAAAGCCGCTCCTACTCGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCg                  <  6:294238/147‑1 (MQ=255)
                                                                                                                     ttCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCTATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCa                                     <  7:68066/128‑1 (MQ=255)
                                                                                                                     ttCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCg                  <  6:289279/147‑1 (MQ=255)
                                                                                                                     ttCGCGGCAAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTATCCGGGAAATGGCCGTCGAGGTCGATCAAAAAAGCA‑ATGAAACACACACATTTTTTAGCAAGTGAGCAAAGTGCATTcaca                                     >  8:68066/1‑128 (MQ=255)
                                                                                                                        gcgGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTACATTCCCAGCTTTCATGTAAGCGTCCGAAcc              >  6:112232/1‑148 (MQ=255)
                                                                                                                             tAAAGCCGCTCCTACGGGGACAGAGTAGGCGCGGCTTTAGCAGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCACAGCTATCATGTAAGCGTACGAACCCaaa          >  8:317446/1‑146 (MQ=255)
                                                                                                                               aaGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCa        <  7:237291/147‑1 (MQ=255)
                                                                                                                                 tccGCGCCTACGCGGTCACCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGATTTCATGTAAGCGTACGAACCCAACCAgg      >  7:12131/2‑147 (MQ=38)
                                                                                                                                 gCCGCTCCTCCGGGGACCGCGGAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGaa    <  4:47362/149‑1 (MQ=255)
                                                                                                                                 gCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACGTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAgg      >  7:17243/1‑147 (MQ=38)
                                                                                                                                  ccGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAg   >  2:200894/1‑149 (MQ=39)
                                                                                                                                    gcgccTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCTAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCAGTCCCAGCTTTGATGTAAGCGTCCGAACCCAACCAGGAAGc  >  6:220922/4‑148 (MQ=38)
                                                                                                                                    gCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGc  >  3:184784/1‑148 (MQ=255)
                                                                                                                                      gccTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTc                                                             >  3:143498/2‑87 (MQ=25)
                                                                                                                                      gccTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTc                                                             <  4:143498/86‑1 (MQ=25)
                                                                                                                                                   |                                                                                                                                     
ATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGC  >  NC_002947/6004853‑6005132

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