Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	887,959	G→A	*179* (TAG→TAA)	rcdA →	transcriptional regulator of csgD and ybiJI, autoregulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	887,959	0	G	A	100.0%	117.2 / NA	37	*179* (TAG→TAA)	rcdA	transcriptional regulator of csgD and ybiJI, autoregulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/21);  total (16/21)

AACAATGGTTTGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAG  >  NC_000913/887826‑888087                                                                                                                                      |                                                                                                                                 aaCAATGGTTTGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGAt                                                                                                                             <  1:1222651/139‑1 (MQ=255)                   gAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACaaaa                                                                                                           >  2:361955/1‑139 (MQ=255)                             cgcgcTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGAt                                                                                                                             <  2:299664/111‑1 (MQ=255)                             cgcgcTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGAt                                                                                                                             >  1:299664/1‑111 (MQ=255)                                     aTGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAAc                                                                                         <  2:717219/139‑1 (MQ=255)                                                   gggATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCg                                                                               <  2:1327187/135‑1 (MQ=255)                                                   gggATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCg                                                                               >  1:1327187/1‑135 (MQ=255)                                                         aCGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGa                                                                                                                              >  2:48470/1‑82 (MQ=255)                                                         aCGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGa                                                                                                                              <  1:48470/82‑1 (MQ=255)                                                           gCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGc                                                                   <  1:371828/139‑1 (MQ=255)                                                             tGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTc                                                                 >  2:463093/1‑139 (MQ=255)                                                              gCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTcc                                                                <  2:659240/139‑1 (MQ=255)                                                                    tGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGAt                                                                                                                             >  2:1070765/1‑72 (MQ=255)                                                                    tGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGAt                                                                                                                             <  1:1070765/72‑1 (MQ=255)                                                                        accgaccgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATggg                                                      <  1:361955/139‑1 (MQ=255)                                                                           gaccgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGAc                                                   >  1:821498/1‑139 (MQ=255)                                                                                tAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATgg                                                                                 >  2:1208543/1‑104 (MQ=255)                                                                                tAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATgg                                                                                 <  1:1208543/104‑1 (MQ=255)                                                                                       cTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAg                                       >  2:376995/1‑139 (MQ=255)                                                                                           cgcgcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGt                                   <  2:144129/139‑1 (MQ=255)                                                                                           cgcgcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGt                                   <  2:821498/139‑1 (MQ=255)                                                                                             cgcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGAt                                                                       >  1:1066398/1‑101 (MQ=255)                                                                                             cgcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGAt                                                                       <  2:1066398/101‑1 (MQ=255)                                                                                              gcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAAcc                                                                                               <  1:63179/76‑1 (MQ=255)                                                                                              gcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAAcc                                                                                               >  2:63179/1‑76 (MQ=255)                                                                                               cgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa                               <  2:1165052/139‑1 (MQ=255)                                                                                               cgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa                               <  2:388872/139‑1 (MQ=255)                                                                                                    agaTTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAAc                                                                                         <  2:259321/76‑1 (MQ=255)                                                                                                    agaTTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAAc                                                                                         >  1:259321/1‑76 (MQ=255)                                                                                                                 ggTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa                               >  1:946514/1‑121 (MQ=255)                                                                                                                 ggTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa                               <  2:946514/121‑1 (MQ=255)                                                                                                                  gTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGccc            <  1:463093/139‑1 (MQ=255)                                                                                                                            ttGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa                               <  2:368057/110‑1 (MQ=255)                                                                                                                            ttGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa                               >  1:368057/1‑110 (MQ=255)                                                                                                                            ttGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAg  >  1:604781/1‑139 (MQ=255)                                                                                                                             tGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAAcc                                                                                               >  2:1039863/1‑45 (MQ=255)                                                                                                                             tGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAAcc                                                                                               <  1:1039863/45‑1 (MQ=255)                                                                                                                                      |                                                                                                                                 AACAATGGTTTGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAG  >  NC_000913/887826‑888087

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