Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,691,192	G→A	Q447Q (CAG→CAA)	ydgA →	DUF945 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,691,192	0	G	A	100.0%	132.0 / NA	42	Q447Q (CAG→CAA)	ydgA	DUF945 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/23);  total (19/23)

GAAACTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAGATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGTTAGAAGATTTTGTTGGCATGTTTGCGATGCCAGCTCTTA  >  NC_000913/1691063‑1691325                                                                                                                                  |                                                                                                                                      gAAACTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGt                                                                                                                              <  2:983779/139‑1 (MQ=255)  aaaCTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTc                                                                                                                             >  1:603992/1‑139 (MQ=255)    aCTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTg                                                                                                                           >  1:1157166/1‑139 (MQ=255)     cTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGa                                                                                                                          <  1:139807/139‑1 (MQ=255)     cTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGa                                                                                                                          <  1:184183/139‑1 (MQ=255)        aCCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                                       <  2:725489/139‑1 (MQ=255)                 gTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAgg                                                                                                              <  2:1157166/139‑1 (MQ=255)                       atGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATa                                                                                                        >  2:1108477/1‑139 (MQ=255)                         ggCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAAt                                                                                                         <  2:51924/136‑1 (MQ=255)                         ggCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAAt                                                                                                         >  1:51924/1‑136 (MQ=255)                              cTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                                       <  2:244214/117‑1 (MQ=255)                              cTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                                       <  1:363539/117‑1 (MQ=255)                              cTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                                       >  1:244214/1‑117 (MQ=255)                              cTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                                       >  2:363539/1‑117 (MQ=255)                                   tttATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGTAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATcacca                                                                                               <  2:239323/136‑1 (MQ=255)                                   tttATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGTAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATcacca                                                                                               >  1:239323/1‑136 (MQ=255)                                   tttATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAgg                                                                                                              >  2:834238/1‑121 (MQ=255)                                   tttATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGACAAATGTTCCGTCTGACCACCTTGCAgg                                                                                                              <  1:834238/121‑1 (MQ=255)                                             aGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTg                                                                                                                  >  1:1284349/1‑107 (MQ=255)                                             aGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTg                                                                                                                  <  2:1284349/107‑1 (MQ=255)                                                    gAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAAc                                                                           <  2:1176507/139‑1 (MQ=255)                                                         tGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCa                                                                      >  1:1058823/1‑139 (MQ=255)                                                             aGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATa                                                                  <  2:1154786/139‑1 (MQ=255)                                                                          gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                                       <  1:637402/73‑1 (MQ=255)                                                                          gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                                       >  1:180499/1‑73 (MQ=255)                                                                          gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                                       <  2:180499/73‑1 (MQ=255)                                                                          gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                                       >  2:637402/1‑73 (MQ=255)                                                                          gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATa                                                                                                        >  2:334477/1‑88 (MQ=255)                                                                          gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATa                                                                                                        <  1:334477/88‑1 (MQ=255)                                                                          gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAg                                                                                           <  2:1116890/101‑1 (MQ=255)                                                                          gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAg                                                                                           >  1:1116890/1‑101 (MQ=255)                                                                               aGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGa                                                                                                             >  1:598168/1‑78 (MQ=255)                                                                               aGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGa                                                                                                             <  2:598168/78‑1 (MQ=255)                                                                               aGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGc                                                                                      <  2:1338185/101‑1 (MQ=255)                                                                               aGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGc                                                                                      >  1:1338185/1‑101 (MQ=255)                                                                                 cGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCTTTACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAAt                                              <  2:603992/139‑1 (MQ=255)                                                                                                                   aTCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTg                                                                                                                  >  2:407109/1‑37 (MQ=255)                                                                                                                   aTCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTg                                                                                                                  <  1:407109/37‑1 (MQ=255)                                                                                                                        cAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGTTAGAAGATTTTGTTgg                         >  2:262339/1‑121 (MQ=255)                                                                                                                        cAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGTTAGAAGATTTTGTTgg                         <  1:262339/121‑1 (MQ=255)                                                                                                                             gggCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGTTAGAAGATTTTGTTGGCATGTTTGCGATGCCAGCTCTTa  >  1:942822/1‑139 (MQ=255)                                                                                                                             gggCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGTTAGAAGATTTTGTTGGCATGTTTGCGATGCCAGCTCTTa  <  2:942822/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                      GAAACTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAGATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGTTAGAAGATTTTGTTGGCATGTTTGCGATGCCAGCTCTTA  >  NC_000913/1691063‑1691325

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: ‑