Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,102,221	T→C	E107E (GAA→GAG)	csgE ←	curlin secretion specificity factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,102,221	0	T	C	100.0%	205.8 / NA	58	E107E (GAA→GAG)	csgE	curlin secretion specificity factor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/27);  total (31/27)

GAAATAAGCATGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGC  >  NC_000913/1102087‑1102350                                                                                                                                       |                                                                                                                                  gAAATAAGCATGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGcctctt                                                                                                                               <  1:521556/139‑1 (MQ=255)  aaaTAAGCATGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCctcttc                                                                                                                              >  2:428918/1‑139 (MQ=255)           tGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                     <  1:428918/139‑1 (MQ=255)             aGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                     <  2:764638/137‑1 (MQ=255)             aGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                     >  1:764638/1‑137 (MQ=255)             aGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCa                                                                                                                   >  2:1118091/1‑139 (MQ=255)              gTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATc                                                                                                                    >  2:1147797/1‑137 (MQ=255)              gTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATc                                                                                                                    <  1:1147797/137‑1 (MQ=255)                      actGCATGTTTGACACGCATGGCAGCCTCCGAACAATTTTTTATTTAGAATTCCTCATTCGCCCAGTCGCCTGTACTTAATAACGCCTGATTTATCTGGCGACTATTTAGTGTCTCCGCAATTTGAATCAGTAGAAAgt                                                                                                          >  1:1176272/1‑138 (MQ=255)                       ctGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                     >  1:1167487/1‑127 (MQ=255)                       ctGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                     <  2:1167487/127‑1 (MQ=255)                       ctGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAgac                                                                                                         >  1:1306706/1‑139 (MQ=255)                       ctGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAgac                                                                                                         >  1:1318340/1‑139 (MQ=255)                            tGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAg                                                                                                    >  1:1554262/1‑139 (MQ=255)                                               ccTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCtttt                                                                                 >  2:846105/1‑139 (MQ=255)                                                cTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAgacga                                                                                                       <  2:1374335/116‑1 (MQ=255)                                                cTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAgacga                                                                                                       >  1:1374335/1‑116 (MQ=255)                                                 tCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                     >  2:1264886/1‑101 (MQ=255)                                                 tCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                     <  1:1264886/101‑1 (MQ=255)                                                    gAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGaa                                                                                                                      <  2:1325801/97‑1 (MQ=255)                                                    gAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGaa                                                                                                                      >  1:1325801/1‑97 (MQ=255)                                                                          catcatGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAAt                                                      <  1:1118091/139‑1 (MQ=255)                                                                                gcgcCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTc                                                >  2:1082262/1‑139 (MQ=255)                                                                                    cAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGa                                            <  2:502415/139‑1 (MQ=255)                                                                                     aaaTCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAt                                           >  1:1329342/1‑139 (MQ=255)                                                                                      aaTCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                                          <  2:547109/139‑1 (MQ=255)                                                                                      aaTCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                                          <  1:1082262/139‑1 (MQ=255)                                                                                         cGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGAc                                       <  2:1306706/139‑1 (MQ=255)                                                                                          gCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGcctctt                                                                                                                               >  1:323050/1‑50 (MQ=255)                                                                                          gCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGcctctt                                                                                                                               <  2:1345010/50‑1 (MQ=255)                                                                                          gCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGcctctt                                                                                                                               <  2:323050/50‑1 (MQ=255)                                                                                          gCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGcctctt                                                                                                                               >  1:1345010/1‑50 (MQ=255)                                                                                           cccGTACTTAATAACGCCTGATTTATCTGGCGCCGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCg                                     >  1:1406870/1‑139 (MQ=255)                                                                                            ccGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGt                                    >  1:1074377/1‑139 (MQ=255)                                                                                                aCTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGaa                                                                                                                      <  2:137019/53‑1 (MQ=255)                                                                                                aCTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGaa                                                                                                                      >  1:137019/1‑53 (MQ=255)                                                                                                     ataaCGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAAc                                                                             >  2:62846/1‑89 (MQ=255)                                                                                                     ataaCGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAAc                                                                             <  1:62846/89‑1 (MQ=255)                                                                                                      taaCGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAgacgac                                                                                                      <  2:534709/63‑1 (MQ=255)                                                                                                      taaCGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAgacgac                                                                                                      >  1:534709/1‑63 (MQ=255)                                                                                                         cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGt                                    >  1:978949/1‑126 (MQ=255)                                                                                                         cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAg                       >  1:814700/1‑139 (MQ=255)                                                                                                         cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAg                       >  2:886124/1‑139 (MQ=255)                                                                                                         cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGAAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                                          >  2:1087152/1‑120 (MQ=255)                                                                                                         cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGAAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                                          <  1:1087152/120‑1 (MQ=255)                                                                                                         cGCCTGATTTATCTGGCGACGATTTAATGCCTCTTCAGTTTGAATCAGTGCAACGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACATCCTGATTGACCGt                                    <  2:978949/126‑1 (MQ=255)                                                                                                          gCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGc                      >  2:1234917/1‑139 (MQ=255)                                                                                                           ccTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCt                     >  2:783972/1‑139 (MQ=255)                                                                                                            cTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCtt                    >  1:377507/1‑139 (MQ=255)                                                                                                              gATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTcc                  <  1:1007654/139‑1 (MQ=255)                                                                                                              gATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTcc                  <  1:1341099/139‑1 (MQ=255)                                                                                                              gATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCCTTATAGTGATCCAGCTTcc                  <  2:1074377/139‑1 (MQ=255)                                                                                                                   aTCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATc             <  2:1554262/139‑1 (MQ=255)                                                                                                                    tCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACgg                                                                           <  2:722773/76‑1 (MQ=255)                                                                                                                    tCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACgg                                                                           >  1:722773/1‑76 (MQ=255)                                                                                                                    tCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGACTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCg            >  1:1054702/1‑139 (MQ=255)                                                                                                                        gCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATcc                         <  2:216281/122‑1 (MQ=255)                                                                                                                        gCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATcc                         >  1:216281/1‑122 (MQ=255)                                                                                                                         cgacgaTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGtt                                                                                                  >  2:1348946/1‑48 (MQ=255)                                                                                                                         cgacgaTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGtt                                                                                                  <  1:1348946/48‑1 (MQ=255)                                                                                                                          gacgaTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGt                                                                                        >  2:1571255/1‑57 (MQ=255)                                                                                                                          gacgaTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGt                                                                                        <  1:1571255/57‑1 (MQ=255)                                                                                                                              aTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGc  <  1:846105/139‑1 (MQ=255)                                                                                                                              aTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGc  <  2:814700/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                  GAAATAAGCATGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGC  >  NC_000913/1102087‑1102350

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