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%	184.5 / NA	57	Q447Q (CAG→CAA)	ydgA	DUF945 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/31);  total (26/31)

GGATGCGAAACTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAGATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGTTAGAAGATTTTGTTGGCATGTTTGCGATGCCAGCTCTTAA  >  NC_000913/1691057‑1691326                                                                                                                                         |                                                                                                                                       ggATGCGAAACTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATg                                                                                                                                     >  1:1359455/1‑139 (MQ=255) ggATGCGAAACTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATg                                                                                                                                     <  2:1359455/139‑1 (MQ=255)      cGAAACTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCg                                                                                                                                >  1:317114/1‑139 (MQ=255)      cGAAACTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCg                                                                                                                                >  2:131698/1‑139 (MQ=255)              aCCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                                        <  1:131698/139‑1 (MQ=255)                  ttCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCtt                                                                                                                    <  1:649085/139‑1 (MQ=255)                   tCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                                        <  2:479825/134‑1 (MQ=255)                   tCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                                        >  1:479825/1‑134 (MQ=255)                   tCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTg                                                                                                                   <  2:167478/139‑1 (MQ=255)                   tCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTg                                                                                                                   <  2:1283295/139‑1 (MQ=255)                          gATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACa                                                                                                            >  1:1040075/1‑139 (MQ=255)                               ggCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATAcc                                                                                                       >  2:754436/1‑139 (MQ=255)                                 cAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCAt                                                                                                     <  1:754436/139‑1 (MQ=255)                                    cTGAGTTTATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTACATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATcac                                                                                                  >  1:759554/1‑139 (MQ=255)                                         tttATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCAccacca                                                                                             <  2:1415681/139‑1 (MQ=255)                                         tttATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCAccacca                                                                                             >  2:85634/1‑139 (MQ=255)                                           tATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGc                                                                                           <  2:759554/139‑1 (MQ=255)                                              gACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTg                                                                                        <  1:21682/139‑1 (MQ=255)                                              gACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTg                                                                                        <  2:1040075/139‑1 (MQ=255)                                                   aGGTAGCGAAGCTGGAAGGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAat                                                                                    >  1:128454/1‑139 (MQ=255)                                                   aGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAata                                                                                   >  2:1034267/1‑139 (MQ=255)                                                     gTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATg                                                                                 <  1:839773/139‑1 (MQ=255)                                                       aGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAATCTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCt                                                                               <  2:317114/139‑1 (MQ=255)                                                         cGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTaa                                                                             >  1:387727/1‑139 (MQ=255)                                                         cGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTaa                                                                             <  1:85634/139‑1 (MQ=255)                                                          gAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAAc                                                                            <  2:1314389/139‑1 (MQ=255)                                                           aaGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACg                                                                           >  1:104590/1‑139 (MQ=255)                                                               tGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCa                                                                       >  2:857044/1‑139 (MQ=255)                                                                ggAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGcc                                                                                          <  2:1335/119‑1 (MQ=255)                                                                ggAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGcc                                                                                          >  1:1335/1‑119 (MQ=255)                                                                   a‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATa                                                                   <  2:441395/139‑1 (MQ=255)                                                                                 gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGt                                                                                                                                    >  2:338998/1‑61 (MQ=255)                                                                                 gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGt                                                                                                                                    <  1:338998/61‑1 (MQ=255)                                                                                 gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTcc                                                                                                                                 <  2:1307481/64‑1 (MQ=255)                                                                                 gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTcc                                                                                                                                 >  1:1307481/1‑64 (MQ=255)                                                                                 gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTg                                                                                                                   >  2:134161/1‑78 (MQ=255)                                                                                 gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTg                                                                                                                   <  1:134161/78‑1 (MQ=255)                                                                                 gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGc                                                                                       >  1:412345/1‑106 (MQ=255)                                                                                 gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGc                                                                                       <  2:412345/106‑1 (MQ=255)                                                                                 gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGc                                                      >  2:1361875/1‑139 (MQ=255)                                                                                        cGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAAt                                               >  1:401091/1‑139 (MQ=255)                                                                                         gAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATg                                              >  2:1272671/1‑139 (MQ=255)                                                                                            aaaaCTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTg                                                                                                                   <  2:969436/67‑1 (MQ=255)                                                                                            aaaaCTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTg                                                                                                                   >  1:969436/1‑67 (MQ=255)                                                                                            aaaaCTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCg                                           >  1:139689/1‑139 (MQ=255)                                                                                              aaCTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGtt                                         <  2:387727/139‑1 (MQ=255)                                                                                                   gCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGAc                                                                                                             >  2:440073/1‑66 (MQ=255)                                                                                                   gCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGAc                                                                                                             <  1:440073/66‑1 (MQ=255)                                                                                                      aaaCAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGTTAGAAGAtt                                 <  2:128454/139‑1 (MQ=255)                                                                                                             aaGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGTTAGAAGATTTTGTTgg                          <  2:139689/139‑1 (MQ=255)                                                                                                                 tGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGTTAGAAGATTTTGTTGGCATg                      <  1:1272671/139‑1 (MQ=255)                                                                                                                    aGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGa                                                                     >  1:978257/1‑89 (MQ=255)                                                                                                                    aGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGa                                                                     <  2:978257/89‑1 (MQ=255)                                                                                                                      gTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGTTAGAAGATTTTGTTGGCATGTTTGc                 <  2:401091/139‑1 (MQ=255)                                                                                                                          aTCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACgg                                                                          >  1:1181155/1‑78 (MQ=255)                                                                                                                          aTCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACgg                                                                          <  2:1181155/78‑1 (MQ=255)                                                                                                                                     ggCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGTTAGAAGATTTTGTTGGCATGTTTGCGATGCCAGCTCTTaa  <  1:419522/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       GGATGCGAAACTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAA‑GGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAGATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGTTAGAAGATTTTGTTGGCATGTTTGCGATGCCAGCTCTTAA  >  NC_000913/1691057‑1691326

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