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%	160.9 / NA	47	E107E (GAA→GAG)	csgE	curlin secretion specificity factor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/25);  total (22/25)

GAAATAAGCATGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCCT  >  NC_000913/1102087‑1102352                                                                                                                                       |                                                                                                                                    gAAATAAGCATGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGcctctt                                                                                                                                 <  2:273997/139‑1 (MQ=255)         cATGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGa                                                                                                                         >  2:1148879/1‑139 (MQ=255)           tGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                       <  1:422679/139‑1 (MQ=255)           tGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                       <  2:1035598/139‑1 (MQ=255)           tGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                       <  2:1106415/139‑1 (MQ=255)            gAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATc                                                                                                                      <  1:1031241/139‑1 (MQ=255)             aGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCa                                                                                                                     >  1:404879/1‑139 (MQ=255)                  aaCTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCa                                                                                                                >  2:716363/1‑139 (MQ=255)                       ctGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                       >  2:301360/1‑127 (MQ=255)                       ctGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                       <  1:301360/127‑1 (MQ=255)                       ctGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAgac                                                                                                           >  1:376307/1‑139 (MQ=255)                            tGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAg                                                                                                      >  1:498216/1‑139 (MQ=255)                                  acacGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTtct                                                                                                >  1:464125/1‑139 (MQ=255)                                        aTTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGt                                                                                          <  1:165252/139‑1 (MQ=255)                                                  ccGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCaa                                                                                <  1:636058/139‑1 (MQ=255)                                                       cAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGaa                                                                           <  2:376307/139‑1 (MQ=255)                                                            ttttATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATaaa                                                                      <  1:716363/139‑1 (MQ=255)                                                            ttttATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATaaa                                                                      <  1:661831/139‑1 (MQ=255)                                                             tttatttaGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATaaaa                                                                     <  2:498216/139‑1 (MQ=255)                                                             tttatttaGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATaaaa                                                                     <  2:464125/139‑1 (MQ=255)                                                               tatttaGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATaaaaaa                                                                   >  1:519501/1‑139 (MQ=255)                                                                          catcatGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAAt                                                        <  2:519501/139‑1 (MQ=255)                                                                                gcgcCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTc                                                  >  2:976412/1‑139 (MQ=255)                                                                                gcgcCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTc                                                  >  1:1218495/1‑139 (MQ=255)                                                                                gcgcCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTc                                                  <  1:186610/139‑1 (MQ=255)                                                                                     aaaTCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                       >  1:235505/1‑65 (MQ=255)                                                                                     aaaTCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                       <  2:235505/65‑1 (MQ=255)                                                                                      aaTCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                                            <  1:976412/139‑1 (MQ=255)                                                                                      aaTCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                                            <  2:1218495/139‑1 (MQ=255)                                                                                           cccGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTATCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCg                                       >  2:824000/1‑139 (MQ=255)                                                                                                         cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAg                                                                                                             >  2:976345/1‑55 (MQ=255)                                                                                                         cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAg                                                                                                             <  1:976345/55‑1 (MQ=255)                                                                                                         cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATaaaa                                                                     <  2:626600/95‑1 (MQ=255)                                                                                                         cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATaaaa                                                                     >  1:626600/1‑95 (MQ=255)                                                                                                         cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATcc                           <  2:548708/137‑1 (MQ=255)                                                                                                         cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATcc                           >  1:548708/1‑137 (MQ=255)                                                                                                          gCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTAATCCAGc                        >  1:1127795/1‑139 (MQ=255)                                                                                                           ccTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCt                       >  2:104024/1‑139 (MQ=255)                                                                                                            cTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGt                                                                 >  2:150359/1‑96 (MQ=255)                                                                                                            cTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGt                                                                 <  1:150359/96‑1 (MQ=255)                                                                                                            cTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAAt                                                        >  1:587651/1‑105 (MQ=255)                                                                                                            cTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAAt                                                        <  2:587651/105‑1 (MQ=255)                                                                                                              gATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTcc                    >  1:24388/1‑139 (MQ=255)                                                                                                                        gCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCa          <  1:104024/139‑1 (MQ=255)                                                                                                                           acgaTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGt                                                                 <  2:40302/81‑1 (MQ=255)                                                                                                                           acgaTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGt                                                                 >  1:40302/1‑81 (MQ=255)                                                                                                                           acgaTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTg       >  1:1205196/1‑139 (MQ=255)                                                                                                                                ttAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCCt  <  1:824000/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                    GAAATAAGCATGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCCT  >  NC_000913/1102087‑1102352

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