Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,850,890	C→T	7.5%	G124G (GGC→GGT)	PP_3398 →	curli fiber surface‑exposed nucleator CsgB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,850,890	0	C	T	7.5%	152.5 / 2.1	53	G124G (GGC→GGT)	PP_3398	curli fiber surface‑exposed nucleator CsgB
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base T (2/2);  total (27/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 = 2.92e-01

TCGTCCAGTCCGGCGGTGCCCAGGAAGCGTACATCTTGCAGGAAGGCAGTGACCTGATGGCCATCATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGC  >  NC_002947/3850745‑3851037                                                                                                                                                  |                                                                                                                                                     tCGTCCAGTCCGGCGGTGCCCAGGAAGCGTACATCTTGCAGGAAGGCAGTGACCTGATGGCCATCATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGAAACAACGCCGCCATCGAGCAGATCGGCAAc                                                                                                                                                   >  5:190897/1‑149 (MQ=255)        gTCCGGCGGTGCCCAGGAAGCGTACATCTTGCAGGAAGGCAGTGACCTGATGGCCATCATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGacaaca                                                                                                                                            <  6:316407/149‑1 (MQ=255)         tCCGGCGGTGCCCAGGAAGCGTACATCTTGCAGGAAGGCAGTGACCTGATGGCCATCATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAg                                                                                                                                           >  5:201443/1‑149 (MQ=255)                gTGCCCAGGAAGCGTACATCTTGCAGGAAGGCAGTGACCTGATGGCCATCATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGc                                                                                                                                   <  4:134320/149‑1 (MQ=255)                   cccAGGAAGCGTACATCTTGCAGGAAGGCAGTGACCTGATGGCCATCATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATc                                                                                                                                <  4:11061/149‑1 (MQ=255)                    ccAGGAAGCGTACATCTTGCAGGAAGGCAGTGACCTGATGGCCATCATCAGCCAAAAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCg                                                                                                                               >  2:325098/1‑149 (MQ=255)                     cAGGAAGCGTACATCTTGCAGGAAGGCAGGGACCTGATGGCCATCATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGt                                                                                                                              <  2:95512/149‑1 (MQ=255)                               aCATCTTGCAGGAAGGCAGTGACCTGATGGCCATCATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGc                                                                                                                    >  2:111037/1‑149 (MQ=255)                                          gAAGGCAGTGACCTGATGGCCATCATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTg                                                                                                                             <  1:321865/129‑1 (MQ=255)                                          gAAGGCAGTGACCTGATGGCCATCATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCa                                                                                                          >  3:161687/1‑148 (MQ=255)                                          gAAGGCAGTGACCTGATGGCCATCATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTAGCCAGAGCGGCAGCAGAAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGACCGGCTCAGGGCTCa                                                                                                          >  7:331387/1‑148 (MQ=255)                                                            gCCATCATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAgg                                                                                        <  4:161687/148‑1 (MQ=255)                                                              catcatCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTg                                                                                     <  1:325098/149‑1 (MQ=255)                                                               atcatcAGGCAACAGGGTAAAGGCAACATTGCGTGGATTCGCCAGTGGGGCATGATCAAAAACGCCGCCATCGAGCATTTCGGCAACGAAAACAGCGCCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGTCTg                                                                                     <  1:84537/149‑1 (MQ=255)                                                                   tcaGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGTAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGAccc                                                                                           <  6:1667/138‑1 (MQ=255)                                                                   tcaGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGTAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGAccc                                                                                           >  5:6787/1‑138 (MQ=255)                                                                   tcaGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGTAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGAccc                                                                                           <  6:6787/138‑1 (MQ=255)                                                                   tcaGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGTAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGAccc                                                                                           >  5:1667/1‑138 (MQ=255)                                                                     aGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATgg                                                                              >  6:246582/1‑149 (MQ=255)                                                                       ccAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGcc                                                                            >  4:116608/1‑149 (MQ=255)                                                                        cAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCa                                                                           <  5:246582/149‑1 (MQ=255)                                                                         aaCAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAg                                                                          >  8:248817/1‑149 (MQ=255)                                                                         aaCAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAg                                                                          >  3:118034/1‑149 (MQ=255)                                                                              ggAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGt                                                                     <  6:201443/149‑1 (MQ=255)                                                                                       aaCAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAc                                                            <  4:118034/149‑1 (MQ=255)                                                                                           gTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCaa                                                        >  5:174798/1‑149 (MQ=255)                                                                                              cGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATAc                                                     <  6:174798/149‑1 (MQ=255)                                                                                               gTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTgg                                                                                    <  3:144099/117‑1 (MQ=255)                                                                                               gTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTgg                                                                                    >  4:144099/1‑117 (MQ=255)                                                                                               gTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATAcc                                                    >  6:243618/1‑149 (MQ=255)                                                                                                  gATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGa                                                  >  3:352017/1‑148 (MQ=255)                                                                                                      cGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGAc                                             >  3:133391/1‑149 (MQ=255)                                                                                                       gCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGAt                                                 <  4:314727/144‑1 (MQ=255)                                                                                                       gCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGAt                                                 >  3:314727/1‑144 (MQ=255)                                                                                                       gCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGAcc                                            >  5:322178/1‑149 (MQ=255)                                                                                                       gCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGAcc                                            >  7:174153/1‑149 (MQ=255)                                                                                                         cAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTg                                          <  5:243618/149‑1 (MQ=255)                                                                                                         cAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTg                                          >  1:353884/1‑149 (MQ=255)                                                                                                          agagCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAc                                                            >  2:68315/1‑130 (MQ=255)                                                                                                          agagCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAc                                                            <  1:68315/130‑1 (MQ=255)                                                                                                             gCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAgg                                      >  2:55355/1‑149 (MQ=255)                                                                                                                           aacGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTc                                                                                                                  <  7:238387/59‑1 (MQ=255)                                                                                                                           aacGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTc                                                                                                                  >  8:238387/1‑59 (MQ=255)                                                                                                                           aacGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTAc                         >  6:54378/1‑148 (MQ=255)                                                                                                                            acGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACaa                       <  4:133391/149‑1 (MQ=255)                                                                                                                            acGCCGCCATCGAGCAGATCGGCAACGACAACACCG‑CCAGCATCGTGCAGTCCCGCTCAGGGCTCTACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACaa                       <  5:367632/149‑1 (MQ=255)                                                                                                                            acGCAGCCATCGAGCAGATCGCCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCGCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACaa                       <  6:190897/149‑1 (MQ=255)                                                                                                                               ccgccATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCAGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAgct                    >  6:176300/1‑149 (MQ=255)                                                                                                                                       gAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGt            <  2:353884/149‑1 (MQ=255)                                                                                                                                       gAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGt            <  1:55355/149‑1 (MQ=255)                                                                                                                                                ggCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCTTGATGTTCAGATCACCCAATACCGATAGAGCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCg   <  1:111037/149‑1 (MQ=255)                                                                                                                                                ggCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCg   <  7:248817/149‑1 (MQ=255)                                                                                                                                                 gCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGc  >  5:300630/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     TCGTCCAGTCCGGCGGTGCCCAGGAAGCGTACATCTTGCAGGAAGGCAGTGACCTGATGGCCATCATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCG‑CCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGC  >  NC_002947/3850745‑3851037

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