Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,085,824	G→A	Q30* (CAG→TAG)	speA ←	biosynthetic arginine decarboxylase, PLP‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,085,824	0	G	A	100.0%	220.5 / NA	68	Q30* (CAG→TAG)	speA	biosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/32);  total (36/32)

GCGCGAGATCGACGCGAGCTTCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTAT  >  NC_000913/3085692‑3085958                                                                                                                                     |                                                                                                                                       gcgcGAGATCGACGCGAGCTTCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCt                                                                                                                                  >  1:125613/1‑139 (MQ=255)          cGACGCGAGCTTCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAAc                                                                                                                         >  2:356699/1‑139 (MQ=255)               cgAGCTTCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAAcctcct                                                                                                                    >  1:1497821/1‑139 (MQ=255)                         ggACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCg                                                                                                          <  1:457604/139‑1 (MQ=255)                          gACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCgt                                                                                                         <  2:125613/139‑1 (MQ=255)                               ccgggtccgggCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTac                                                                                                    >  1:829178/1‑139 (MQ=255)                                cgggtccgggCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTaca                                                                                                   <  2:1497821/139‑1 (MQ=255)                                cgggtccgggCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTaca                                                                                                   >  2:1388115/1‑139 (MQ=255)                                        ggCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTg                                                                                                                             <  2:254305/105‑1 (MQ=255)                                        ggCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTg                                                                                                                             >  1:254305/1‑105 (MQ=255)                                         gCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCg                                                                                          <  2:829178/139‑1 (MQ=255)                                                       ggCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAgg                                                                            >  2:702107/1‑139 (MQ=255)                                                         cccAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCa                                                                          <  2:1570815/139‑1 (MQ=255)                                                                           tagtagTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATgtcgt                                                        <  1:356699/139‑1 (MQ=255)                                                                              tagtTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAg                                                     <  1:702107/139‑1 (MQ=255)                                                                                   aTTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGc                                                                                                                                   >  1:1509666/1‑56 (MQ=255)                                                                                   aTTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGc                                                                                                                                   <  2:1509666/56‑1 (MQ=255)                                                                                     tGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCa                                                                          >  2:1307169/1‑111 (MQ=255)                                                                                     tGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCa                                                                          <  1:1307169/111‑1 (MQ=255)                                                                                      gCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCa                                                                          >  1:928859/1‑110 (MQ=255)                                                                                      gCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCa                                                                          <  2:928859/110‑1 (MQ=255)                                                                                      gCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCg                                             >  2:1373751/1‑139 (MQ=255)                                                                                        cccACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGaa                                           >  2:681176/1‑139 (MQ=255)                                                                                          caccaGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCgt                                                                                                         >  1:267871/1‑75 (MQ=255)                                                                                          caccaGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCgt                                                                                                         <  2:267871/75‑1 (MQ=255)                                                                                          caccaGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGc                                              >  2:85323/1‑134 (MQ=255)                                                                                          caccaGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGc                                              <  1:85323/134‑1 (MQ=255)                                                                                          caccaGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAAcg                                         >  1:1291285/1‑138 (MQ=255)                                                                                           accaGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACAt                                                             >  1:1646378/1‑118 (MQ=255)                                                                                           accaGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACAt                                                             <  2:1646378/118‑1 (MQ=255)                                                                                           accaGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCt                                        >  1:1155981/1‑139 (MQ=255)                                                                                               ggCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACGTCaaa                                    <  2:1291285/139‑1 (MQ=255)                                                                                                 cAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAatt                                  <  1:1336926/139‑1 (MQ=255)                                                                                                         tAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCaaa                                    <  2:33163/129‑1 (MQ=255)                                                                                                         tAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCaaa                                    >  1:33163/1‑129 (MQ=255)                                                                                                         tAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAAttattttatt                          >  1:379875/1‑139 (MQ=255)                                                                                                         tAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAAttattttatt                          <  1:681176/139‑1 (MQ=255)                                                                                                         tAAGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAAttattttatt                          <  1:1373751/139‑1 (MQ=255)                                                                                                          aaGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCgtag                                                                                                       <  2:763804/61‑1 (MQ=255)                                                                                                          aaGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCgtag                                                                                                       >  1:763804/1‑61 (MQ=255)                                                                                                           aGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGc                                                                                                                            >  1:687390/1‑39 (MQ=255)                                                                                                           aGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGc                                                                                                                            <  2:687390/39‑1 (MQ=255)                                                                                                           aGTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTaa                        >  1:638438/1‑139 (MQ=255)                                                                                                            gTACGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAg                       >  2:1461740/1‑139 (MQ=255)                                                                                                              aCGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATg                                                                                                               <  1:321298/49‑1 (MQ=255)                                                                                                              aCGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCa                                                                          >  1:1468554/1‑86 (MQ=255)                                                                                                              aCGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCa                                                                          <  2:1468554/86‑1 (MQ=255)                                                                                                              aCGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATgtcg                                                         <  2:505274/103‑1 (MQ=255)                                                                                                              aCGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATgtcg                                                         >  1:505274/1‑103 (MQ=255)                                                                                                              aCGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACa                                                  >  2:176509/1‑110 (MQ=255)                                                                                                              aCGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACa                                                  <  1:176509/110‑1 (MQ=255)                                                                                                              aCGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAgtg                     >  1:83589/1‑139 (MQ=255)                                                                                                              aCGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAgtg                     >  2:428924/1‑139 (MQ=255)                                                                                                              aCGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAg                       >  1:447749/1‑137 (MQ=255)                                                                                                              aCGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAg                       <  2:447749/137‑1 (MQ=255)                                                                                                              aCGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCAGAATg                                                                                                               >  2:321298/1‑49 (MQ=255)                                                                                                               cGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATg                                                            <  2:643670/99‑1 (MQ=255)                                                                                                               cGCAGCATCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATg                                                            >  1:643670/1‑99 (MQ=255)                                                                                                               cGCAGCATCTTGCTGGCTTCCTAGCAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAgtgt                    >  1:1012188/1‑139 (MQ=255)                                                                                                                   gcaTCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTaaaa                <  2:1155981/139‑1 (MQ=255)                                                                                                                    caTCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAAc               <  1:428924/139‑1 (MQ=255)                                                                                                                     aTCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCa                                                                          >  1:780735/1‑79 (MQ=255)                                                                                                                     aTCTTGCTGGCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCa                                                                          <  2:780735/79‑1 (MQ=255)                                                                                                                              gCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAAttattttatt                          <  1:1340258/118‑1 (MQ=255)                                                                                                                              gCTTCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAAttattttatt                          >  2:1340258/1‑118 (MQ=255)                                                                                                                                 tCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGATGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTAt  >  1:983131/1‑139 (MQ=255)                                                                                                                                 tCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATgtcgt                                                        <  2:692775/85‑1 (MQ=255)                                                                                                                                 tCCTAGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATgtcgt                                                        >  1:692775/1‑85 (MQ=255)                                                                                                                                     |                                                                                                                                       GCGCGAGATCGACGCGAGCTTCCGGGACGTCCGGGTCCGGGCACACGCTAATGTGGCCCAGCTCGTTAACGTCATAGTAGTTATTGCCCCACCAGGCAATATTGTAAGTACGCAGCATCTTGCTGGCTTCCTGGGAGCTCATTGCAACCTCCTGCATGGAGCGTAGTACACCGTGTTCGCCCGCTGACGAAGGCAAACCCATAGACATGTCGTCAGACATAGCGAACCTCAAATTATTTTATTAAGTGTAAAACAGTTAACGACTAT  >  NC_000913/3085692‑3085958

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