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

GCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGC  >  NC_000913/1104816‑1105077
                                                                                                                                      |                                                                                                                                
gctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGgag                                                                                                                             <  1:227553/139‑1 (MQ=255)
 ctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGgaga                                                                                                                            <  1:460480/139‑1 (MQ=255)
         gACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATg                                                                                                                    >  2:507250/1‑139 (MQ=255)
            cAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGgag                                                                                                                             >  1:226580/1‑127 (MQ=255)
            cAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGgag                                                                                                                             <  2:226580/127‑1 (MQ=255)
                             ctccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGc                                                                                                >  2:206095/1‑139 (MQ=255)
                              tccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCg                                                                                               <  2:358119/139‑1 (MQ=255)
                               ccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCgg                                                                                              >  2:457661/1‑139 (MQ=255)
                               ccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCgg                                                                                              <  1:461257/139‑1 (MQ=255)
                                          aCTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGt                                                                                   >  2:531608/1‑139 (MQ=255)
                                                          gTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATAc                                                                  <  1:216964/139‑1 (MQ=255)
                                                              caacGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGAt                                                                              >  1:89381/1‑124 (MQ=255)
                                                              caacGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGAt                                                                              <  2:89381/124‑1 (MQ=255)
                                                               aacGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAcc                                                              <  1:531608/139‑1 (MQ=255)
                                                                        gCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAggg                                                     <  1:206095/139‑1 (MQ=255)
                                                                            agaaGTACTTATTCATCATTTGTATTACAGAAACTGGGCGCAAGCCCTGTTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCTTTTTCCAGTCAGATAACCTTTAATACGACCCTGCAAGGGGAt                                                  <  1:429581/136‑1 (MQ=255)
                                                                                    tAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTgg                                                                                                                                <  2:484416/52‑1 (MQ=255)
                                                                                    tAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATAcc                                         >  1:201475/1‑139 (MQ=255)
                                                                                      aTACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCat                                       >  1:175600/1‑139 (MQ=255)
                                                                                           tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGgag                                                                                                                             <  1:524954/48‑1 (MQ=255)
                                                                                           tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGgag                                                                                                                             >  2:524954/1‑48 (MQ=255)
                                                                                           tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGa                                                   <  1:474073/122‑1 (MQ=255)
                                                                                           tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGa                                                   >  2:474073/1‑122 (MQ=255)
                                                                                           tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTc                                  >  2:69803/1‑139 (MQ=255)
                                                                                           tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTc                                  >  1:423089/1‑139 (MQ=255)
                                                                                             atTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGgag                                                                                                                             <  1:264870/46‑1 (MQ=255)
                                                                                             atTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGgag                                                                                                                             >  2:264870/1‑46 (MQ=255)
                                                                                                tGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                             >  1:168151/1‑75 (MQ=255)
                                                                                                tGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                             <  2:168151/75‑1 (MQ=255)
                                                                                                  tATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGt                           <  2:175600/139‑1 (MQ=255)
                                                                                                    ttACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAgggg                                                    >  1:249093/1‑112 (MQ=255)
                                                                                                    ttACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAgggg                                                    <  2:249093/112‑1 (MQ=255)
                                                                                                         gAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAgggg                                                    >  2:284885/1‑107 (MQ=255)
                                                                                                         gAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAgggg                                                    <  1:284885/107‑1 (MQ=255)
                                                                                                         gAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAAGGGATGTGTATACCATTATTCCTGAAGTCACTCtt                    <  1:507250/139‑1 (MQ=255)
                                                                                                           aaCAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTAc                  <  2:201475/139‑1 (MQ=255)
                                                                                                               gggCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATa                                                                             <  2:239938/76‑1 (MQ=255)
                                                                                                               gggCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATa                                                                             >  1:239938/1‑76 (MQ=255)
                                                                                                               gggCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAg                                                       >  1:81413/1‑98 (MQ=255)
                                                                                                               gggCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAg                                                       <  2:81413/98‑1 (MQ=255)
                                                                                                                ggCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACt                                                                                                      <  2:205578/50‑1 (MQ=255)
                                                                                                                ggCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACt                                                                                                      >  1:205578/1‑50 (MQ=255)
                                                                                                                    cAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGGGGCACTTTCCAGTTAGATAACCTTTAATACGAACCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTg         >  2:273590/1‑139 (MQ=255)
                                                                                                                          cTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCtgtg   <  2:423089/139‑1 (MQ=255)
                                                                                                                          cTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCtgtg   <  1:457661/139‑1 (MQ=255)
                                                                                                                           tGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGc  <  1:69803/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                
GCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGC  >  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: