Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,222,8160TG50.0% 16.2 / 10.5 32V205G (GTC→GGC) oprJouter membrane protein OprJ
Reads supporting (aligned to +/- strand):  ref base T (10/6);  new base G (0/16);  total (10/22)
Fisher's exact test for biased strand distribution p-value = 2.48e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAGCGCCAGGCCTACAACGCGCTGGTGCTGT  >  NC_002947/3222681‑3222950
                                                                                                                                       |                                                                                                                                      
aGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGGGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGt                                                                                                                                        >  2:188778/1‑136 (MQ=255)
aGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGCCCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGt                                                                                                                                        >  2:182557/1‑136 (MQ=255)
aGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGt                                                                                                                                        >  1:68738/1‑136 (MQ=255)
aGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGt                                                                                                                                        >  2:86561/1‑136 (MQ=255)
aGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGt                                                                                                                                        >  2:230032/1‑136 (MQ=255)
aGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGt                                                                                                                                        >  2:73761/1‑136 (MQ=255)
        gggCGTGCCGCACGAATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCCGGCGCAGACCCGACCGACCCAGGCCCGCCGCGACGATACCCAGGCGATGGTCGGCCAgc                                                                                                                                >  1:288518/1‑136 (MQ=255)
                cgcacgcaTCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGcgcc                                                                                                                        >  2:129124/1‑136 (MQ=255)
                           cccTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGcacggc                                                                                                             >  2:360290/1‑136 (MQ=255)
                                  cgccgcGGTGAGCCTGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGgcg                                                                                                      <  2:38280/136‑1 (MQ=255)
                                        ggTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGac                                                                                                >  2:186129/1‑136 (MQ=255)
                                            aGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGGCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGaccacc                                                                                            <  2:304174/136‑1 (MQ=255)
                                             cccctgCCTCCCTGGCCTGGGCCGCCGCCGCGCGGCGCCAGGCGCGGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGGCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGAccacca                                                                                           <  1:189782/131‑1 (MQ=255)
                                               cAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCCGCCCGGCCACGGCGCGGGTCCACCaaa                                                                                         >  1:326003/1‑135 (MQ=255)
                                                      aCCGGGCCTGGGCCGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCGGGCGTTGGGCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCt                                                                                  <  1:254911/136‑1 (MQ=255)
                                                       ccTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGCCAGCCCCTGGCCAGCCGCGAGGTTTCCCTGGCGTGGGGCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTg                                                                                 <  1:27268/136‑1 (MQ=255)
                                                        cTGGCGTTGCACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCCGACCCTGGCCAGCCGCGAGGAGTCCCTGGCGTTGGGCGGCCAGCGGCGCGCGGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTgg                                                                                <  1:336625/136‑1 (MQ=255)
                                                        cTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACCGCCACGGCGCTGGGCCACCAAGAGGCGCAgg                                                                                >  1:179077/1‑136 (MQ=255)
                                                          ggCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCGGCGGCGCGCCCCCGGCCCGGCCACGGCGCGGGACCCCCAACAGGCGCTgggg                                                                              >  2:38380/1‑136 (MQ=255)
                                                           gCCTTGGACGGCGCCGAGCGGCGCCAGGCGCGGACCCGGCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGGCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGc                                                                             <  1:251128/136‑1 (MQ=255)
                                                                ggcgggCGCCGCGCGGCGCCTGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGGCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTc                                                                        <  1:214240/132‑1 (MQ=255)
                                                                ggACGGCGCCGTGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTc                                                                        <  2:8818/136‑1 (MQ=255)
                                                                ggACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGGCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTc                                                                        <  1:86561/136‑1 (MQ=255)
                                                                    ggCGCCGCGCGGCGCCAGGCGCGGACCCGTCAGACCCGGGCCAGCCGCGAGGATTCCCTGGCGTTGGGCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGc                                                                    <  1:194941/136‑1 (MQ=255)
                                                                                 cccAGGCGCAGACCCGTCTGACCCGGGCCTGCCGCGAGGATTCCCTGGCGTTGGTCGGCCTGCGGCGGGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGa                                                       <  1:242382/135‑1 (MQ=255)
                                                                                    gggCGCTGCCCCGCCACCCCCGGGCCCGCCGCGAGGATGCCCGGGCGTTGGGCGGCCCGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCt                                                    <  2:289182/135‑1 (MQ=255)
                                                                                                    caCCCTGGCCAGCCGCGAGGCTGCCCGGGCGGGGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGTAGTCCCGTGCCGAGCTCGAAGTCACCGGgcgc                                    <  2:274635/135‑1 (MQ=255)
                                                                                                         tGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCCCGGCGCTGGACCCCCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAAcc                               >  2:333439/1‑135 (MQ=255)
                                                                                                         tGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCCCCGCCCCGCCCGCCAAGGCCCGGGGCCCCCAAAAGGCGCAGGGGGGGGTCGAGCACTCCCGATCCGCCCGCGAAGTGACCGGGCGCCcaca                               >  2:127035/1‑136 (MQ=255)
                                                                                                          ggCCAGCCGCGAGGATTCCCTGGCGTTGGGCGGCCAGCGGCGCGCCGCCGGCATGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAg                              <  1:236865/136‑1 (MQ=255)
                                                                                                          ggCCAGCCGCGAGGATGCCCTGGCGTGGGGCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAg                              <  1:73761/136‑1 (MQ=255)
                                                                                                          ggCCAGCCGCGAGGATCCCCTGGCGTGGGGCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAg                              <  1:112523/136‑1 (MQ=255)
                                                                                                               gCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAGCGCCa                         <  1:170142/136‑1 (MQ=255)
                                                                                                                        aTTCCCTGGCGTTGGGCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAGCGCCAGGCCTACaa                <  1:191797/136‑1 (MQ=255)
                                                                                                                                gCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAGCGCCAGGCCTACAACGCGCTgg        <  2:242274/136‑1 (MQ=255)
                                                                                                                                gCGTTGGGCGGCCTGCGGCGCGCCGCCGGCACGGCCTCGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAGCGCCAGGCCTACAACGCGCTgg        <  2:29824/136‑1 (MQ=255)
                                                                                                                                  gTTGGGCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAGCGCCAGGCCTACAACGCGCTGGTg      <  2:79626/136‑1 (MQ=255)
                                                                                                                                     ggTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGCGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCCTGCGCCAACACCGCC‑GGCCTACAACGCGCTGGTGCGGt  >  1:330999/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                                      
AGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAGCGCCAGGCCTACAACGCGCTGGTGCTGT  >  NC_002947/3222681‑3222950

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