Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,015,5160GT25.6% 199.7 / 31.0 82intergenic (+57/‑199)pqiC/rmfintermembrane transport lipoprotein PqiC/ribosome modulation factor
Reads supporting (aligned to +/- strand):  ref base G (24/37);  new base T (21/0);  total (45/37)
Fisher's exact test for biased strand distribution p-value = 2.46e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACC  >  NC_000913/1015468‑1015565
                                                |                                                 
aTTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGt                                                  <  1:7398721/50‑1 (MQ=255)
aTTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGt                                                  >  1:1783858/1‑50 (MQ=255)
 tttGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGtc                                                 <  1:3598448/50‑1 (MQ=255)
   tGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTtctc                                               >  1:3760317/1‑50 (MQ=255)
   tGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTtctc                                               >  1:3321415/1‑50 (MQ=255)
   tGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGtct                                                <  1:2536504/49‑1 (MQ=255)
    gTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTtctc                                               >  1:1891695/1‑49 (MQ=255)
    gTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGtctc                                               >  1:3455721/1‑49 (MQ=255)
     tAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTtctcc                                              >  1:3009924/1‑48 (MQ=255)
     tAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTttct                                              >  1:4200461/1‑49 (MQ=38)
     tAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTccct                                              >  1:3197581/1‑49 (MQ=255)
      aaaTATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTTTCTaa                                            >  1:4280231/1‑50 (MQ=39)
      aaaTATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTaa                                            >  1:3437978/1‑50 (MQ=255)
      aaaTATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTa                                             >  1:4200567/1‑49 (MQ=255)
       aaTATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTaa                                            <  1:7821523/49‑1 (MQ=255)
        atatAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACg                                          <  1:6038050/50‑1 (MQ=255)
         tataACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGt                                         >  1:6169093/1‑50 (MQ=255)
         tataACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGt                                         <  1:2067753/50‑1 (MQ=255)
         tataACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGt                                         <  1:6665543/50‑1 (MQ=255)
         tataACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGt                                         >  1:5841006/1‑50 (MQ=255)
         tataACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGt                                         <  1:8555644/50‑1 (MQ=255)
         tataACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGt                                         <  1:5814960/50‑1 (MQ=255)
         tataACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGt                                         <  1:3557542/50‑1 (MQ=255)
         tataACCGTCTCCGGTATGGTGCCTGAGGCGGTTTTTTTGTCTCTAACGt                                         <  1:797034/50‑1 (MQ=255)
           taACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTTACGTGc                                       >  1:5878052/1‑50 (MQ=39)
           taACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGc                                       >  1:658886/1‑50 (MQ=255)
           taACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTg                                        <  1:4803370/49‑1 (MQ=255)
           taACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTg                                        <  1:2254233/49‑1 (MQ=255)
           taACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGc                                       <  1:7890087/50‑1 (MQ=255)
           taACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGc                                       <  1:4236062/50‑1 (MQ=255)
            aaCCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCg                                      >  1:3626473/1‑50 (MQ=255)
            aaCCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCg                                      >  1:4919720/1‑50 (MQ=255)
            aaCCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTTTCTAACGTGCg                                      >  1:8230408/1‑50 (MQ=255)
            aaCCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTTTCCAACGGGCg                                      >  1:4410045/1‑50 (MQ=38)
            aaCCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGGGCg                                      >  1:1696061/1‑50 (MQ=255)
             aCCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTTTCTTACGTGCgg                                     >  1:2089395/1‑50 (MQ=39)
             aCCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCgg                                     >  1:6152172/1‑50 (MQ=255)
             aCCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCgg                                     <  1:3809666/50‑1 (MQ=255)
              ccGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCCAACGTGCgg                                     >  1:785399/1‑49 (MQ=255)
               cGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACCTGCggg                                    >  1:8535986/1‑48 (MQ=38)
               agTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTTACGTGCggg                                    >  1:177816/2‑48 (MQ=39)
                gTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCggga                                   >  1:6780555/1‑47 (MQ=255)
                gTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGGCTCTAACGTGCGGaaa                                  <  1:5467209/50‑1 (MQ=255)
                 tctcCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGaaaa                                 <  1:5009262/50‑1 (MQ=255)
                 tctcCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGaaaa                                 <  1:5435160/50‑1 (MQ=255)
                  ctcCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCGGaaaaa                                >  1:8326823/1‑50 (MQ=255)
                   tcCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCGGGAAAAt                               >  1:700869/1‑50 (MQ=39)
                    ccGGTATGTTGCCTGAGGCGGTTTTTTTTTCTTCAACGGCCGGaaaaaca                              >  1:3362799/1‑48 (MQ=25)
                    ccGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCGGGAAAAtt                              >  1:1383315/1‑50 (MQ=39)
                    ccGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCGGAAAAAtt                              >  1:3055019/1‑50 (MQ=255)
                     cGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGaaaaa                                <  1:6120157/47‑1 (MQ=255)
                     cGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAAttt                             <  1:7858694/50‑1 (MQ=255)
                     cGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAAttt                             <  1:1386945/50‑1 (MQ=255)
                     cGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAAttt                             <  1:3723232/50‑1 (MQ=255)
                        tATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCGGGAAAAtttttt                          >  1:3930387/1‑50 (MQ=38)
                         aTGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCGGGAAAAATTGttt                         >  1:4569082/1‑49 (MQ=38)
                         aTGTTGCCTGAGGCGGTTTTTTTGTCTCTAAAGTGCGGAAAAATTTGTTc                         >  1:2030597/1‑50 (MQ=255)
                          tGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTcc                        <  1:238345/50‑1 (MQ=255)
                          tGTTGCCTGAGGCGGTTTCTTTGTCTCTAACGTGCGGAAAAATTTGTTcc                        <  1:5637094/50‑1 (MQ=255)
                            ttGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCctc                      <  1:3116241/50‑1 (MQ=255)
                              gCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCctct                     <  1:2196382/49‑1 (MQ=255)
                              gCCTGAGGCGGTTTTTTTGTCTCTAACCTGCGGAAAAATTTGTTCCTCtt                    >  1:4460431/1‑50 (MQ=255)
                               ccTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCtt                    <  1:1027788/49‑1 (MQ=255)
                                   aGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTcaca                <  1:6520365/49‑1 (MQ=255)
                                   aGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTcaca                <  1:7885433/49‑1 (MQ=255)
                                   aGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAt               <  1:8515594/50‑1 (MQ=255)
                                   aGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAt               <  1:6594330/50‑1 (MQ=255)
                                    ggCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAtt              >  1:5559410/1‑50 (MQ=255)
                                    ggCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAtt              >  1:3871258/1‑50 (MQ=255)
                                     gCGGTTTTTTTGTCTCTAACGGGCGGAAAAATTTGTTCCTCTTCACAttt             <  1:7296376/50‑1 (MQ=255)
                                      cGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAtttt            >  1:6050695/1‑50 (MQ=255)
                                      cGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAtttt            <  1:3365621/50‑1 (MQ=255)
                                       ggTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAttttt           >  1:1043625/1‑50 (MQ=255)
                                       ggTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAttttt           >  1:1481393/1‑50 (MQ=255)
                                        gTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAtttttt          >  1:2456415/1‑50 (MQ=255)
                                         tttttttGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTg         >  1:3431863/1‑50 (MQ=255)
                                          ttttttGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTg         >  1:5893594/1‑49 (MQ=255)
                                           tttttGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGt        <  1:4619211/49‑1 (MQ=255)
                                           tttttGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGt        <  1:6059817/49‑1 (MQ=255)
                                             tttGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTAc      >  1:4392542/1‑49 (MQ=255)
                                              ttGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACaa    >  1:5914334/1‑50 (MQ=255)
                                              ttGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACa     <  1:2093078/49‑1 (MQ=255)
                                                gTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAAcc  >  1:7730911/1‑50 (MQ=255)
                                                gTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAAc   <  1:6615980/49‑1 (MQ=255)
                                                |                                                 
ATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACC  >  NC_000913/1015468‑1015565

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