Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,104,949 C→T intergenic (+47/‑12) csgA → / → csgC curlin subunit, amyloid curli fibers, cryptic/curli assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,104,9490CT97.8% 156.4 / ‑3.3 46intergenic (+47/‑12)csgA/csgCcurlin subunit, amyloid curli fibers, cryptic/curli assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/25);  minor base . (0/1);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

GCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGC  >  NC_000913/1104816‑1105077
                                                                                                                                      |                                                                                                                                
gctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGgag                                                                                                                             <  1:227553/139‑1 (MQ=255)
 ctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGgaga                                                                                                                            <  1:460480/139‑1 (MQ=255)
         gACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATg                                                                                                                    >  2:507250/1‑139 (MQ=255)
            cAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGgag                                                                                                                             >  1:226580/1‑127 (MQ=255)
            cAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGgag                                                                                                                             <  2:226580/127‑1 (MQ=255)
                             ctccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGc                                                                                                >  2:206095/1‑139 (MQ=255)
                              tccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCg                                                                                               <  2:358119/139‑1 (MQ=255)
                               ccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCgg                                                                                              >  2:457661/1‑139 (MQ=255)
                               ccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCgg                                                                                              <  1:461257/139‑1 (MQ=255)
                                          aCTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGt                                                                                   >  2:531608/1‑139 (MQ=255)
                                                          gTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTT‑GGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATAc                                                                  <  1:216964/139‑1 (MQ=255)
                                                              caacGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGAt                                                                              >  1:89381/1‑124 (MQ=255)
                                                              caacGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGAt                                                                              <  2:89381/124‑1 (MQ=255)
                                                               aacGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAcc                                                              <  1:531608/139‑1 (MQ=255)
                                                                        gCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAggg                                                     <  1:206095/139‑1 (MQ=255)
                                                                            agaaGTACTTATTCATCATTTGTATTACAGAAACTGGGCGCAAGCCCTGTTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCTTTTTCCAGTCAGATAACCTTTAATACGACCCTGCAAGGGGAt                                                  <  1:429581/136‑1 (MQ=255)
                                                                                    tAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTgg                                                                                                                                <  2:484416/52‑1 (MQ=255)
                                                                                    tAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATAcc                                         >  1:201475/1‑139 (MQ=255)
                                                                                      aTACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCat                                       >  1:175600/1‑139 (MQ=255)
                                                                                           tcatTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGgag                                                                                                                             <  1:524954/48‑1 (MQ=255)
                                                                                           tcatTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGgag                                                                                                                             >  2:524954/1‑48 (MQ=255)
                                                                                           tcatTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGa                                                   <  1:474073/122‑1 (MQ=255)
                                                                                           tcatTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGa                                                   >  2:474073/1‑122 (MQ=255)
                                                                                           tcatTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTc                                  >  2:69803/1‑139 (MQ=255)
                                                                                           tcatTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTc                                  >  1:423089/1‑139 (MQ=255)
                                                                                             atTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGgag                                                                                                                             <  1:264870/46‑1 (MQ=255)
                                                                                             atTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGgag                                                                                                                             >  2:264870/1‑46 (MQ=255)
                                                                                                tGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                             >  1:168151/1‑75 (MQ=255)
                                                                                                tGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                             <  2:168151/75‑1 (MQ=255)
                                                                                                  tATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGt                           <  2:175600/139‑1 (MQ=255)
                                                                                                    ttACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAgggg                                                    >  1:249093/1‑112 (MQ=255)
                                                                                                    ttACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAgggg                                                    <  2:249093/112‑1 (MQ=255)
                                                                                                         gAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAgggg                                                    >  2:284885/1‑107 (MQ=255)
                                                                                                         gAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAgggg                                                    <  1:284885/107‑1 (MQ=255)
                                                                                                         gAAACAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAAGGGATGTGTATACCATTATTCCTGAAGTCACTCtt                    <  1:507250/139‑1 (MQ=255)
                                                                                                           aaCAGGGCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTAc                  <  2:201475/139‑1 (MQ=255)
                                                                                                               gggCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATa                                                                             <  2:239938/76‑1 (MQ=255)
                                                                                                               gggCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATa                                                                             >  1:239938/1‑76 (MQ=255)
                                                                                                               gggCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAg                                                       >  1:81413/1‑98 (MQ=255)
                                                                                                               gggCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAg                                                       <  2:81413/98‑1 (MQ=255)
                                                                                                                ggCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACt                                                                                                      <  2:205578/50‑1 (MQ=255)
                                                                                                                ggCGCAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACt                                                                                                      >  1:205578/1‑50 (MQ=255)
                                                                                                                    cAAGCCCTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGGGGCACTTTCCAGTTAGATAACCTTTAATACGAACCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTg         >  2:273590/1‑139 (MQ=255)
                                                                                                                          cTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCtgtg   <  2:423089/139‑1 (MQ=255)
                                                                                                                          cTG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCtgtg   <  1:457661/139‑1 (MQ=255)
                                                                                                                           tG‑TTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGc  <  1:69803/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                
GCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGC  >  NC_000913/1104816‑1105077

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