Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,798,3940CG40.0% 41.7 / 61.4 50F207L (TTC→TTG) bamA‑Iouter membrane protein assembly factor
Reads supporting (aligned to +/- strand):  ref base C (11/19);  new base G (0/20);  total (11/39)
Fisher's exact test for biased strand distribution p-value = 1.62e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCGCGCAACCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACAAGTA  >  NC_002947/1798245‑1798411
                                                                                                                                                     |                 
ccGCGCAACCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGttcttc                   <  2:1453348/150‑1 (MQ=32)
ccGCGCAACCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGttcttc                   <  1:9583447/150‑1 (MQ=32)
ccGCGCAACCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGttcttc                   <  1:5816485/150‑1 (MQ=32)
ccGCGCAACCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGttcttc                   >  2:7280168/1‑150 (MQ=32)
 cgcgCAACCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCa                  <  1:5732346/150‑1 (MQ=32)
  gcgcAACCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCTCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCaa                 <  2:1396910/150‑1 (MQ=17)
  gcgcAACCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCaa                 <  1:4720446/150‑1 (MQ=32)
  gcgcAACCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCaa                 <  1:48903/150‑1 (MQ=32)
  gcgcAACCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCaa                 >  2:3942729/1‑150 (MQ=32)
    gcAACCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGa               <  1:232682/150‑1 (MQ=32)
     cAACCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCCGTCGTTCTTgaagaa              <  2:1942952/150‑1 (MQ=17)
      aaCCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAac             >  2:4640796/1‑150 (MQ=32)
      aaCCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAac             <  1:9393678/150‑1 (MQ=32)
      aaCCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAac             <  2:417172/150‑1 (MQ=32)
      aaCCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAac             >  2:4641002/1‑150 (MQ=32)
      aaCCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAac             >  1:4962575/1‑150 (MQ=32)
      aaCCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAac             <  2:7847624/150‑1 (MQ=32)
        ccGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacga           <  2:7689773/150‑1 (MQ=17)
         cGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAacgac          <  1:4792777/150‑1 (MQ=32)
          gtgtGGCCCTGAAGATCAAGATCAACGACGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca         <  2:2582640/150‑2 (MQ=17)
          gtgtGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCACCATCGTTGGCAACAACGTATTCGATGATGAGACCCTGTGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca         <  2:2582583/150‑2 (MQ=12)
          gtgtGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca         <  1:2569256/150‑2 (MQ=17)
          gtgtGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca         <  2:4118912/150‑2 (MQ=17)
          gtgtGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca         <  2:6560798/150‑2 (MQ=17)
          gtgtGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca         <  2:3675323/150‑2 (MQ=17)
          gtgtGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca         <  2:3675191/150‑2 (MQ=17)
          gtgtGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca         <  2:3675180/150‑2 (MQ=17)
          gtgtGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca         <  2:8557345/150‑2 (MQ=17)
          gtgtGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca         <  2:2582491/150‑2 (MQ=17)
          gtgtGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca         <  1:7973376/150‑2 (MQ=17)
          gtgtGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca         <  1:7574839/150‑2 (MQ=17)
          gtgtGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca         <  1:7188858/150‑2 (MQ=17)
          gtgtGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCCGTCGTTCTTGAAGAacgaca         <  2:95834/150‑2 (MQ=17)
          gtgtGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCCGTCGTTCTTGAAGAacgaca         <  2:5655174/150‑2 (MQ=17)
          gtgtGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCCGTCGTTCTTGAAGAacgaca         <  2:5655318/150‑2 (MQ=17)
           tgtgTCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACCACGTATTCGATGTTGAGACCCTGGGGCTGCTGTTCGAGATGAAGACCTCCAACTGGATTTCGTTCTTCATGAacgtcta        <  2:8873184/150‑5 (MQ=11)
           tgtgGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgacag        <  1:7989409/150‑3 (MQ=17)
           tgtgGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAAcgacga        <  2:8872790/150‑1 (MQ=32)
           tgtgGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAAcgacga        <  1:6665253/150‑1 (MQ=32)
           tgtgGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAAcgacga        <  2:7110287/150‑1 (MQ=32)
           tgtgGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAAcgacga        <  1:6665156/150‑1 (MQ=32)
            gtgGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACCTCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACgacgac       <  2:5844373/150‑1 (MQ=17)
            gtgGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgacagc       <  1:8847774/150‑4 (MQ=17)
            gtgGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgacagc       <  1:4042177/150‑4 (MQ=17)
            gtgGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACgacgac       >  2:5984082/1‑150 (MQ=32)
            gtgGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAAAATCGTTGGCAACAACGTATTCGATGATGATACCCTGCGGAAGCTGTTCGAGCTGAAGACCACAAACTGGCTGTCGTTCTTCAAGAACgacgac       >  2:5984376/1‑150 (MQ=11)
              ggCCCTGATGTTCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTTAATACCACCAACGGGGTGTCGTTCTTCAAGAACGACGACaa     <  1:4123356/150‑1 (MQ=11)
              ggCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACaa     <  2:3137649/150‑1 (MQ=32)
               gCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAacgac          <  2:8699074/144‑1 (MQ=32)
               gCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAacgac          >  1:8699074/1‑144 (MQ=32)
               gCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACAAg    >  1:6492702/1‑150 (MQ=32)
               gCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACAAg    <  1:4911038/150‑1 (MQ=32)
                 ccTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACAAGTa  >  2:574232/1‑150 (MQ=11)
                 ccTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACAAGTa  >  1:7201349/1‑150 (MQ=11)
                 ccTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACAAGTa  >  1:7201391/1‑150 (MQ=11)
                                                                                                                                                     |                 
CCGCGCAACCGTGTGGCCCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACAAGTA  >  NC_002947/1798245‑1798411

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