Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 882,839 T→C 21.2% intergenic (‑44/‑40) ybjJ ← / → rcdA putative transporter/transcriptional regulator of csgD and ybiJI; autoregulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273882,8390TC21.2% 23.3 / 5.0 19intergenic (‑44/‑40)ybjJ/rcdAputative transporter/transcriptional regulator of csgD and ybiJI; autoregulator
Reads supporting (aligned to +/- strand):  ref base T (11/4);  new base C (2/2);  total (13/6)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-02

TTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAAT  >  CP009273/882699‑882979
                                                                                                                                            |                                                                                                                                            
ttGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGCAGACAACATTAGCAAGGTTGTGTACAt                                                                                                                                     >  1:154208/1‑150 (MQ=255)
    ccAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGt                                                                                                                                 >  1:84931/1‑150 (MQ=255)
                             aaaaaGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACCATTGCCGGGAGAGGACATGa                                                                                                        >  1:249656/1‑150 (MQ=255)
                             aaaaaGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGCTGTGTACATTTGTACACAATTGCAGAGAGAGGaa                                                                                                            >  2:216973/1‑146 (MQ=255)
                             aaaaaGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGCTGTGTACATTTGTACACAATTGCAGAGAGAGGaa                                                                                                            >  2:216925/1‑146 (MQ=255)
                                     gaaCATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGAATACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTAGACAATTGCAGAGAGAGGAAATGACATGCGTc                                                                                                >  1:176628/1‑150 (MQ=255)
                                        cATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTcgcg                                                                                             >  1:214674/1‑150 (MQ=255)
                                                    ccAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGGGAGAGTAAATCACGTGCGTAGTGCTAACGATCCGc                                                                                 >  1:80206/1‑150 (MQ=255)
                                                             gCGTTTCAATGCATTACGTAAAAAATTTACGGTCATAGATCACCTCAAAATAAGAAGAGGGAAAACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCgag                                                                        <  1:2464/150‑1 (MQ=255)
                                                               gTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCg                                                                          >  2:223395/1‑146 (MQ=255)
                                                                tttCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCga                                                                         >  2:117603/1‑146 (MQ=255)
                                                                            aCGTGAAGAATTGACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATc                                                             <  2:155355/146‑1 (MQ=255)
                                                                                gaagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCATCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCAc                                                     >  1:115051/1‑150 (MQ=255)
                                                                                                       cGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGCTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACg                              <  1:243616/150‑1 (MQ=255)
                                                                                                       cGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGCTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACg                              <  1:243566/150‑1 (MQ=255)
                                                                                                                  gcagagGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCt                   <  1:220105/148‑1 (MQ=255)
                                                                                                                        ggAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGt                 >  2:214806/1‑146 (MQ=255)
                                                                                                                                 cttcaGCAAGTTTTGGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCg        <  2:218795/142‑1 (MQ=255)
                                                                                                                                       cAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAAt  >  2:85039/1‑146 (MQ=255)
                                                                                                                                            |                                                                                                                                            
TTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAAT  >  CP009273/882699‑882979

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