Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,103,590	A→G	intergenic (‑394/‑361)	csgD ← / → csgB	DNA‑binding transcriptional dual regulator CsgD/curlin, minor subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,103,590	0	A	G	100.0%	186.2 / NA	49	intergenic (‑394/‑361)	csgD/csgB	DNA‑binding transcriptional dual regulator CsgD/curlin, minor subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/15);  total (34/15)

GTACATTTGGTTTTTATTGCACAATTTTAAAAAATCATACAAATGGTGATAACTTACTAATAATGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATT  >  NC_000913/1103463‑1103721                                                                                                                                |                                                                                                                                    gTACATTTGGTTTTTATTGCACAATTTTAAAAAATCATACAAATGGTGATAACTTACTAATAATGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGa                                                                                                                           <  1:1157152/138‑1 (MQ=255) gTACATTTGGTTTTTATTGCACAATTTTAAAAAATCATACAAATGGTGATAACTTACTAATAATGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGa                                                                                                                           >  2:1157152/1‑138 (MQ=255)      tttGGTTTTTATTGCACAATTTTAAAAAATCATACAAATGGTGATAACTTACTAATAATGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAAtttt                                                                                                                    >  1:744284/1‑140 (MQ=255)         ggTTTTTATTGCACAATTTTAAAAAATCATACAAATGGTGATAACTTACTAATAATGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTa                                                                                                                 >  1:801235/1‑140 (MQ=255)         ggTTTTTATTGCACAATTTTAAAAAATCATACAAATGGTGATAACTTACTAATAATGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTa                                                                                                                 <  1:446771/140‑1 (MQ=255)               tATTGCACAATTTTAAAAAATCATACAAATGGAGATAACTTACTAATAATGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGATCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCaa                                                                                                           >  2:144996/1‑140 (MQ=255)                      caATTTTAAAAAATCATACAAATGGTGATAACTTACTAATAATGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACg                                                                                                    >  1:556321/1‑140 (MQ=255)                         ttttAAAAAATCATACAAATGGTGATAACTTACTAATAATGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAgg                                                                                                 >  1:397426/1‑140 (MQ=255)                                   tCATACAAATGGTGATAACTTACTAATAATGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACgg                                                                                       >  2:1053405/1‑140 (MQ=255)                                   tCATACAAATGGTGATAACTTACTAATAATGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACgg                                                                                       >  1:784072/1‑140 (MQ=255)                                      tACAAATGGTGATAACTTACTAATAATGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCtt                                                                                    <  1:432782/140‑1 (MQ=255)                                           aTGGTGATAACTTACTAATAATGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTgcg                                                                               <  1:120390/140‑1 (MQ=255)                                             ggTGATAACTTACTAATAATGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCa                                                                             >  1:724999/1‑140 (MQ=255)                                                     cTTACTAATAATGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACatat                                                                     >  1:1349055/1‑140 (MQ=255)                                                         cTAATAATGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCgca                                                                 <  2:309814/140‑1 (MQ=255)                                                            ataatGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCgcagca                                                              >  1:212759/1‑140 (MQ=255)                                                            ataatGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCgcagca                                                              >  1:378198/1‑140 (MQ=255)                                                                tGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAGCGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCa                                                          >  2:546161/1‑140 (MQ=255)                                                                tGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCa                                                          >  1:943993/1‑140 (MQ=255)                                                                tGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCa                                                          >  1:810053/1‑140 (MQ=255)                                                                tGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCa                                                          >  1:692793/1‑140 (MQ=255)                                                                  cATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGc                                                        >  2:1215414/1‑140 (MQ=255)                                                                   atatAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCg                                                       >  1:371620/1‑140 (MQ=255)                                                                    tataAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGa                                                      >  1:1363074/1‑140 (MQ=255)                                                                    tataAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGa                                                      <  1:182127/140‑1 (MQ=255)                                                                        aaaaaTATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACggg                                                  >  2:790303/1‑140 (MQ=255)                                                                         aaaaTATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGc                                                 <  1:1053405/140‑1 (MQ=255)                                                                            atatTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCaag                                              >  2:191541/1‑140 (MQ=255)                                                                               tttCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCaagaag                                           >  2:181834/1‑140 (MQ=255)                                                                                ttCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAagaaga                                          <  1:1372790/140‑1 (MQ=255)                                                                                 tCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAgaagaa                                         >  2:1431651/1‑140 (MQ=255)                                                                                    gtgtAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGa                                      <  2:692793/140‑1 (MQ=255)                                                                                    gtgtAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGa                                      <  2:397426/140‑1 (MQ=255)                                                                                          tCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCt                                <  1:1141856/140‑1 (MQ=255)                                                                                            cTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTgg                              >  2:724915/1‑140 (MQ=255)                                                                                               tCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGc                           >  1:70496/1‑140 (MQ=255)                                                                                               tCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGc                           >  2:653716/1‑140 (MQ=255)                                                                                               tCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGc                           >  2:619014/1‑140 (MQ=255)                                                                                               tCGTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGc                           >  2:1249830/1‑140 (MQ=255)                                                                                                 gTCATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCtt                         >  2:905163/1‑140 (MQ=255)                                                                                                   cATGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGAATGTCTGGTGCtttt                       >  2:1135349/1‑140 (MQ=255)                                                                                                    aTGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCttttt                      >  2:1184533/1‑140 (MQ=255)                                                                                                    aTGTAAAACGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCttttt                      >  1:1020912/1‑140 (MQ=255)                                                                                                        aaaaCGTTCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATa                  <  2:1349055/140‑1 (MQ=255)                                                                                                               tCTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGaaa           >  2:354457/1‑140 (MQ=255)                                                                                                                cTTGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGaaaa          >  2:1326270/1‑140 (MQ=255)                                                                                                                 ttGTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAAc         <  1:619014/140‑1 (MQ=255)                                                                                                                  tctgttttCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACg        <  1:144996/136‑1 (MQ=255)                                                                                                                   gTTTTTTCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACgg       <  2:744284/140‑1 (MQ=255)                                                                                                                        ttCTCCACGCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGAtt  <  2:966141/140‑1 (MQ=255)                                                                                                                                |                                                                                                                                    GTACATTTGGTTTTTATTGCACAATTTTAAAAAATCATACAAATGGTGATAACTTACTAATAATGCATATAAAAAATATTTCGGTGTAGTCCTTTCGTCATGTAAAACGTTCTTGTTTTTTCTCCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTTGTCACGGCTTGTGCGCAAGACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACGGAGATT  >  NC_000913/1103463‑1103721

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