Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,004,9970AT23.5% 78.7 / 25.9 51intergenic (+35/‑141)PP_5259/ydcJLysR family transcriptional regulator/metalloprotein
Reads supporting (aligned to +/- strand):  ref base A (17/22);  new base T (7/5);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.23e-02
Rejected as consensus: Frequency below/above cutoff threshold.

ACCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGC  >  NC_002947/6004850‑6005132
                                                                                                                                                   |                                                                                                                                        
aCCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGAc                                                                                                                                         <  7:143648/149‑1 (MQ=255)
 ccATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGAcc                                                                                                                                        <  2:8932/149‑1 (MQ=255)
   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                                                                                                                            >  4:200298/1‑139 (MQ=37)
                       gcgcATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCgg                                                                                                                            <  3:200298/139‑1 (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                                                                                                                              >  7:270121/1‑98 (MQ=25)
                                                              ccAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGgagc                                                                                                                              <  8:270121/98‑1 (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                                                             <  4:143498/86‑1 (MQ=25)
                                                                                                                                         gccTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTc                                                             >  3:143498/2‑87 (MQ=25)
                                                                                                                                                   |                                                                                                                                        
ACCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGA‑ATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGC  >  NC_002947/6004850‑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: ‑