Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,588,006	C→A	K36N (AAG→AAT)	ugpQ ←	glycerophosphodiester phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,588,006	0	C	A	100.0%	164.2 / NA	47	K36N (AAG→AAT)	ugpQ	glycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (21/26);  total (21/26)

CCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACC  >  NC_000913/3587870‑3588133                                                                                                                                         |                                                                                                                                ccGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTa                                                                                                                              <  2:63793/140‑1 (MQ=255)  cGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTa                                                                                                                              <  1:805568/139‑1 (MQ=255)  cGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTa                                                                                                                              >  2:805568/1‑139 (MQ=255)      tCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGAcc                                                                                                                         <  1:29974/140‑1 (MQ=255)      tCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGAcc                                                                                                                         <  2:216141/140‑1 (MQ=255)               aGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGc                                                                                                                >  1:632201/1‑140 (MQ=255)                         gCCAGTTCAGTTCACCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCg                                                                                                      >  2:877412/1‑140 (MQ=255)                           cagttcagttCACCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGt                                                                                                        <  1:846207/136‑1 (MQ=255)                           cagttcagttCACCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGt                                                                                                        >  2:846207/1‑136 (MQ=255)                           cagttcagttCACCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGAt                                                                                                    >  1:718601/1‑140 (MQ=255)                              ttcagttcaCCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTTGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGa                                                                                                 >  2:859356/1‑140 (MQ=255)                                  gttcaCCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTTGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGcc                                                                                             <  1:859356/140‑1 (MQ=255)                                     caCCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAgg                                                                                          <  1:713071/140‑1 (MQ=255)                                     caCCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAgg                                                                                          <  1:575171/140‑1 (MQ=255)                                     caCCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAgg                                                                                          <  2:330826/140‑1 (MQ=255)                                         cgcgACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGgtgt                                                                                      <  1:625378/140‑1 (MQ=255)                                             aCGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTc                                                                                  >  2:248551/1‑140 (MQ=255)                                                   cAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGc                                                                            <  2:718601/140‑1 (MQ=255)                                                   cAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTGGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGc                                                                            <  2:304045/140‑1 (MQ=255)                                                          tGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTa                                                                     >  1:330557/1‑140 (MQ=255)                                                           gCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTAc                                                                    >  1:702682/1‑140 (MQ=255)                                                           gCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTAc                                                                    <  2:632201/140‑1 (MQ=255)                                                           gCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTAc                                                                    >  2:435330/1‑140 (MQ=255)                                                                    ttCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGcc                                                                           >  2:125886/1‑124 (MQ=255)                                                                    ttCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGcc                                                                           <  1:125886/124‑1 (MQ=255)                                                                     tCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACg                                                          <  2:702682/140‑1 (MQ=255)                                                                      cGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGa                                                         >  1:447713/1‑140 (MQ=255)                                                                       gagaTTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGAt                                                        <  2:447713/140‑1 (MQ=255)                                                                                    tGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCg                                           <  1:877412/140‑1 (MQ=255)                                                                                       aGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCgggg                                        >  2:638836/1‑140 (MQ=255)                                                                                        gCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGa                                       >  2:474471/1‑140 (MQ=255)                                                                                                aTCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCa                               <  2:550844/140‑1 (MQ=255)                                                                                                 tctcGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAg                              >  1:357914/1‑140 (MQ=255)                                                                                                      ccATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCg                                           <  2:90304/122‑1 (MQ=255)                                                                                                      ccATCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCg                                           >  1:90304/1‑122 (MQ=255)                                                                                                         tCTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCAt                      <  2:330557/140‑1 (MQ=255)                                                                                                          cTTTCGATAACTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATa                     >  2:22886/1‑140 (MQ=255)                                                                                                                   aCTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTc            <  1:248551/140‑1 (MQ=255)                                                                                                                    cTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGt                             >  1:746200/1‑122 (MQ=255)                                                                                                                    cTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGt                             <  2:746200/122‑1 (MQ=255)                                                                                                                    cTTCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTcc           >  2:563538/1‑140 (MQ=255)                                                                                                                     ttCGCGTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCt          >  1:777866/1‑140 (MQ=255)                                                                                                                       cgcgTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCg                 <  2:300877/131‑1 (MQ=255)                                                                                                                       cgcgTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCg                 >  1:300877/1‑131 (MQ=255)                                                                                                                         cgTCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGttt      <  1:22886/140‑1 (MQ=255)                                                                                                                           tCAAATTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCa    <  2:793382/140‑1 (MQ=255)                                                                                                                             aaaTTCGATCATATTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCAcc  <  1:474471/140‑1 (MQ=255)                                                                                                                                    atcataTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCa    <  2:467745/131‑1 (MQ=255)                                                                                                                                    atcataTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCa    >  1:467745/1‑131 (MQ=255)                                                                                                                                         |                                                                                                                                CCGCATCCACGCGCAGTAAATCCTGCCAGTTCAGTTCACCCGCGACGCCCCAGCCGTTGCTGGTACGTTCGAGATTGTCGTCATGGAGCAGGAAGATCTCGCCATCTTTCGATAACTTCGCGTCAAATTCGATCATCTTATGACCGTATTTTGCCCCGACGTCGATTGACGCCAGGGTGTTTTCCGGGGCCAGCTTACCGCCGCCACGATGAGCGACGATGCGGGGATAAGGCCAGTTACTCATACTCGTTGTCCTGTTTCACC  >  NC_000913/3587870‑3588133

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