Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,915,7270GA32.8% 158.3 / 128.9 122R224Q (CGG→CAG) barAsensory histidine kinase BarA
Reads supporting (aligned to +/- strand):  ref base G (39/43);  new base A (18/22);  total (57/65)
Fisher's exact test for biased strand distribution p-value = 8.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-03
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 129 total aligned reads displayed.

AACATGGTGAATACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGAAAAACGGTATC  >  NC_000913/2915678‑2915776
                                                 |                                                 
aaCATGGTGAATACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCg                                                   >  1:4778580/1‑50 (MQ=255)
aaCATGGTGAATACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCg                                                   >  1:2175792/1‑50 (MQ=255)
aaCATGGTGAATACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCg                                                   <  2:3683790/50‑1 (MQ=255)
  cATGGTGAATACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCggg                                                 >  2:1358280/1‑50 (MQ=255)
  cATGGTGAATACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCagg                                                 <  1:4931185/50‑1 (MQ=255)
   aTGGTGAATACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCAGGt                                                >  1:4975449/1‑50 (MQ=255)
   aTGGTGAATACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCAGGt                                                <  1:2589749/50‑1 (MQ=255)
   aTGGTGAATACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCAGGt                                                <  2:2551716/50‑1 (MQ=255)
   aTGGTGAATACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCAGGt                                                >  1:180050/1‑50 (MQ=255)
    tGGTGAATACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCGGGTg                                               >  2:3889836/1‑50 (MQ=255)
    tGGTGAATACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCGGGTg                                               >  1:353300/1‑50 (MQ=255)
      gTGAATACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGa                                             <  2:3949952/50‑1 (MQ=255)
      gTGAATACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCAGGTgg                                              <  1:5096818/49‑1 (MQ=255)
       tGAATACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGaa                                            <  1:427382/50‑1 (MQ=255)
       tGAATACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCAGGTgg                                              >  1:1011815/1‑48 (MQ=255)
        gAATACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAg                                           >  1:4056409/1‑50 (MQ=255)
         aaTACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAgg                                          <  1:27923/50‑1 (MQ=255)
         aaTACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAg                                           <  1:844712/49‑1 (MQ=255)
          gtACCGTCGACCGCAGGCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGa                                         <  2:1586893/49‑1 (MQ=255)
          aTACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGa                                         <  1:3351793/50‑1 (MQ=255)
          aTACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGa                                         >  2:1939108/1‑50 (MQ=255)
           tACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGAt                                        >  1:1154748/1‑50 (MQ=255)
           tACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGAt                                        <  1:1650450/50‑1 (MQ=255)
           tACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGAt                                        >  2:569863/1‑50 (MQ=255)
            aCCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCAGGTGGAAGGAtt                                       >  1:2109367/1‑50 (MQ=255)
             ccGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCAGGTGGAAGGAtt                                       <  1:4532465/49‑1 (MQ=255)
              cGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCAGGTGGAAGGATTTa                                     >  1:4167739/1‑50 (MQ=255)
               gTCGACCGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTag                                    <  2:3469357/50‑2 (MQ=255)
               gTCGACCGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTAt                                    >  1:3618017/1‑50 (MQ=255)
               gTCGACCGCATCCGTCGCGGGCAACTCGACAGCCAGGTGGAAGGATTTAt                                    <  2:2526604/50‑1 (MQ=255)
                tCGACCGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTATg                                   >  1:4906414/1‑50 (MQ=255)
                tCGACCGCATCCGTCGCGGGCAACTCGACAGCCAGGTGGAAGGATTTATg                                   <  1:718878/50‑1 (MQ=255)
                 cGACCGCATCCGTCGCGGGCAACTCGACAGCCAGGTGGAAGGATTTATGc                                  >  1:4159190/1‑50 (MQ=255)
                  gACCGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTATGCt                                 >  2:1666174/1‑50 (MQ=255)
                  gACCGCATCCGTCGCGGGCAACTCGACAGCCAGGTGGAAGGATTTATGCt                                 >  2:1549948/1‑50 (MQ=255)
                   aCCGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTc                                <  2:3610997/50‑1 (MQ=255)
                   aCCGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTc                                >  2:152189/1‑50 (MQ=255)
                   aCCGCATCCGTCGCGGGCAACTCGAAAGCCGGGTGGAAGGAGTTATGCTc                                <  2:3738659/50‑1 (MQ=255)
                    ccGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCg                               <  1:39241/50‑1 (MQ=255)
                    ccGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCg                               <  2:2875654/50‑1 (MQ=255)
                    ccGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCg                               >  2:410520/1‑50 (MQ=255)
                    ccGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCg                               <  2:4409429/50‑1 (MQ=255)
                    ccGCATCCGTCGCGGGCAACTCGACAGCCAGGTGGAAGGATTTATGCTCg                               <  1:873649/50‑1 (MQ=255)
                    ccGCATCCGTCGCGGGCAACTCGACAGCCAGGTGGAAGGATTTATGCTCg                               <  2:1706199/50‑1 (MQ=255)
                      gCATCCGTCGCGGGCAACTCGACAGCCAGGTGGAAGGATTTATGCTCGGc                             <  2:606510/50‑1 (MQ=255)
                      gCATCCGTCGCGGGCAACTCGACAGCCAGGTGGAAGGATTTATGCTCGGc                             <  2:4819773/50‑1 (MQ=255)
                      gCATCCGTCGCGGGCAACTCGACAGCCAGGTGGAAGGATTTATGCTCGGc                             <  2:2861143/50‑1 (MQ=255)
                       cATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCg                            <  2:1753960/50‑1 (MQ=255)
                       cATCCGTCGCGGGCAACTCGACAGCCAGGTGGAAGGATTTATGCTCGGCg                            <  2:266233/50‑1 (MQ=255)
                        aTCCGTCGCGGGCAACTCGACAGCCGGGCGGAAGGATTTATGCTCGGc                             >  2:250996/1‑48 (MQ=255)
                         tCCGTCGCGGGCAACTCGACAGCCAGGTGGAAGGATTTATGCTCGGCGAg                          >  1:681451/1‑50 (MQ=255)
                           cGTCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGc                         <  2:746498/49‑1 (MQ=255)
                           cGTCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGc                         <  1:4370922/49‑1 (MQ=255)
                           cGTCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCt                        <  2:3939052/50‑1 (MQ=255)
                           cGTCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCt                        >  2:4044931/1‑50 (MQ=255)
                           cGTCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCt                        >  1:340120/1‑50 (MQ=255)
                            gTCGCGGGCAACTCGACAGCCAGGTGGAAGGATTTATGCTCGGCGAGCTg                       >  2:642499/1‑50 (MQ=255)
                            gTCGCGGGCAACTCGACAGCCAGGTGGAAGGATTTATGCTCGGCGAGCTg                       >  2:55421/1‑50 (MQ=255)
                             tCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTg                       >  2:187685/1‑49 (MQ=255)
                              cgcgGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGa                     >  2:4207096/1‑50 (MQ=255)
                              cgcgGGCAACTCGACAGCCAGGTGGAAGGATTTATGCTCGGCGAGCTGGa                     >  2:4208227/1‑50 (MQ=255)
                               gcgGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGat                    <  1:3856079/50‑1 (MQ=255)
                               gcgGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGa                     <  1:2252058/49‑1 (MQ=255)
                               gcgGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGa                     >  1:1915177/1‑49 (MQ=255)
                               gcgGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGa                     <  1:3242892/49‑1 (MQ=255)
                                cgGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGata                   >  2:1922574/1‑50 (MQ=255)
                                cgGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGata                   <  1:3493830/50‑1 (MQ=255)
                                cgGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGata                   <  2:3618111/50‑1 (MQ=255)
                                cgGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGata                   <  2:4928791/50‑1 (MQ=255)
                                 gggCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGatat                  >  2:529790/1‑50 (MQ=255)
                                 gggCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGatat                  <  2:3879002/50‑1 (MQ=255)
                                 gggCAACTCGACAGCCAGGTGGAAGGATTTATGCTCGGCGAGCTGGatat                  <  2:4071508/50‑1 (MQ=255)
                                 gggCAACTCGACAGCCAGGTGGAAGGATTTATGCTCGGCGAGCTGGatat                  <  2:884829/50‑1 (MQ=255)
                                  ggCAACTCGACAGCCGGGTGTAAGGATTTATGCTCGGCGAGCTGGATATg                 <  1:2788967/50‑1 (MQ=255)
                                  ggCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATg                 <  2:2649346/50‑1 (MQ=255)
                                  ggCAACTCGACAGCCAGGTGGAAGGATTTATGCTCGGCGAGCTGGATATg                 <  2:631646/50‑1 (MQ=255)
                                   gCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGc                <  1:5033002/50‑1 (MQ=255)
                                    cAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCt               <  1:4945646/50‑1 (MQ=255)
                                     aaCTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCt               >  2:3845301/1‑49 (MQ=255)
                                       cTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGa             >  2:5219655/1‑49 (MQ=255)
                                       cTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGa             >  1:5184058/1‑49 (MQ=255)
                                       cTCGACAGCCAGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGa             >  2:3657473/1‑49 (MQ=255)
                                        tCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGaa            >  2:3829258/1‑49 (MQ=255)
                                        tCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGaa            >  1:4045920/1‑49 (MQ=255)
                                        tCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGaa            >  1:126144/1‑49 (MQ=255)
                                         cGACAGCCAGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGaaa           >  1:1777126/1‑49 (MQ=255)
                                          gACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGaaaaa         <  1:2388681/50‑1 (MQ=255)
                                          gACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGaaaa          <  1:3945273/49‑1 (MQ=255)
                                          gACAGCCAGGTGGAAGGATTTATGCTCGGCGGGCTGGATATGCTGaaa           <  2:4177283/48‑1 (MQ=39)
                                          gACAGCCAGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGaaaa          <  2:899244/49‑1 (MQ=255)
                                           aCAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGAAAAAc        <  1:1409399/50‑1 (MQ=255)
                                           aCAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGAAAAAc        <  1:2364984/50‑1 (MQ=255)
                                           aCAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGAAAAAc        <  1:2881373/50‑1 (MQ=255)
                                            cAGCCGGGTGGACGGATTTATGCTCGGCGAGCTGGATATGCTGAACAACg       >  1:215045/1‑50 (MQ=255)
                                             aGCCAGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGAAAAACgg      >  1:4799027/1‑50 (MQ=255)
                                             aGCAGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGAAAAACgg      >  2:2329316/1‑50 (MQ=255)
                                              gCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGAAAAGCGGt     <  2:3425594/50‑1 (MQ=255)
                                               ccGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGAAAAACGGTa    >  1:4670590/1‑50 (MQ=255)
                                                cGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGAAAAACGGTAt   <  2:4788407/50‑1 (MQ=255)
                                                 gggTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGAAAAACGGTATc  <  1:3521790/50‑1 (MQ=255)
                                                 |                                                 
AACATGGTGAATACCGTCGACCGCATCCGTCGCGGGCAACTCGACAGCCGGGTGGAAGGATTTATGCTCGGCGAGCTGGATATGCTGAAAAACGGTATC  >  NC_000913/2915678‑2915776

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