Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,852,080 T→C 19.2% V343A (GTC→GCC)  PP_3399 → minor curlin subunit CsgB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,852,0800TC19.2% 34.3 / 4.3 21V343A (GTC→GCC) PP_3399minor curlin subunit CsgB
Reads supporting (aligned to +/- strand):  ref base T (7/10);  new base C (2/2);  total (9/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

AAAACCTCAACCTGACCATCAACACCCAAGGCAATGGTAACAGCATCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAGGCAACCAGATGTCGACCAAGCAGGATGATGTCAGCAACGA  >  NC_002947/3851942‑3852176
                                                                                                                                          |                                                                                                
aaaaCCTCAACCTGACCATCAACACCCAAGGCAATGGTAACAGCATCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGAccc                                                                                       >  2:70112/1‑150 (MQ=255)
                 aTCAACACCCAAGGCAATGGTAACAGCATCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCaa                                                                      <  2:192541/150‑1 (MQ=255)
                           aaGGCAATGGTAACAGCATCACCGCCAGTAACTGCAAAGGCGCCAAGATGGGTGGGCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAAGGGCAACGGCGTCAATCTCACGCAGAAGGGAGCTGACCAACTGGCCACGc                                                            <  2:57389/150‑1 (MQ=255)
                             ggCAATGGTAACAGCTTCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGCCAATCTGa                                                                                          >  2:576898/1‑118 (MQ=255)
                             ggCAATGGTAACAGCTTCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGCCAATCTGa                                                                                          <  1:576898/118‑1 (MQ=255)
                             ggCAATGGTAACAGCATCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGCCAATCTGa                                                                                          <  1:576726/118‑1 (MQ=255)
                             ggCAATGGTAACAGCATCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGCCAATCTGa                                                                                          >  2:576726/1‑118 (MQ=255)
                                     tAACAGCATCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGAc                                                        <  2:710986/144‑1 (MQ=255)
                                     tAACAGCATCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGAc                                                        >  1:710986/1‑144 (MQ=255)
                                                 cGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAGGCaa                                      <  2:258657/150‑1 (MQ=255)
                                                  gCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAGGCAAc                                     <  2:194848/150‑1 (MQ=255)
                                                  gCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAGGCAAc                                     >  1:194848/1‑150 (MQ=255)
                                                           tGGAAAGGCGCCAAGATGGGTGGCCAGTTAGGCACCAACCAGACCGCAACCACCAACCAGGCTGGCAATGGCCCCGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCAAGACGGATAAAGGCAACCGGACGTcc                            >  2:170862/1‑149 (MQ=255)
                                                                      cAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCt                                                                            <  2:268601/91‑1 (MQ=255)
                                                                      cAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCt                                                                            >  1:268601/1‑91 (MQ=255)
                                                                       aacaTGGGTGGCCAGTTCAGCACAAACCAAAACGCCACCATCATGCTGGCTGGCAATTGCAACGGCGTCAATCTGTATCAGAAGGGCGCTAACCGACTGGCCACGCTGACCCAGACGGGTAAAAGCAACCAGATGTCGACCAAGCAGgat                <  2:144464/147‑1 (MQ=255)
                                                                                   cAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAg                                          >  1:587047/1‑112 (MQ=255)
                                                                                   cAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAg                                          <  2:587047/112‑1 (MQ=255)
                                                                                     gTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAGGCAACCAGATGTCGACCAAGCAGGATGATGTCAGCAACGa  <  1:70112/150‑1 (MQ=255)
                                                                                              caaccaGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAGGCAACCAGATGTCGa                           >  2:414304/1‑116 (MQ=255)
                                                                                              caaccaGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAGGCAACCAGATGTCGa                           <  1:414304/116‑1 (MQ=255)
                                                                                                                                          |                                                                                                
AAAACCTCAACCTGACCATCAACACCCAAGGCAATGGTAACAGCATCACCGCCAGTAACTGGAAAGGCGCCAAGATGGGTGGCCAGTTCGGCACCAACCAGAACGCCACCATCAACCAGGCTGGCAATGGCAACGGCGTCAATCTGACCCAGAAGGGCGCTGACCAACTGGCCACGCTGACCCAGACGGGTAAAGGCAACCAGATGTCGACCAAGCAGGATGATGTCAGCAACGA  >  NC_002947/3851942‑3852176

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