Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,153,230	C→A	intergenic (‑132/+466)	ppc ← / ← argE	phosphoenolpyruvate carboxylase/acetylornithine deacetylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,153,230	0	C	A	100.0%	210.5 / NA	60	intergenic (‑132/+466)	ppc/argE	phosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (31/29);  total (31/29)

TCATATTACCCCAGACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGA  >  NC_000913/4153095‑4153356                                                                                                                                        |                                                                                                                               tcATATTACCCCAGACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGaaa                                                                                                                            >  1:620553/1‑140 (MQ=255) tcATATTACCCCAGACACCCCATCTTATCGTTTAATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGaaa                                                                                                                            >  2:420289/1‑140 (MQ=255)    tatTACCCCAGACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCg                                                                                                                         >  1:58622/1‑140 (MQ=255)      ttACCCCAGACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCt                                                                                                                       >  1:463162/1‑140 (MQ=255)       tACCCCAGACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCtt                                                                                                                      <  1:295707/140‑1 (MQ=255)            cAGACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCa                                                                                                                 <  2:440918/140‑1 (MQ=255)            cAGACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCa                                                                                                                 <  2:488115/140‑1 (MQ=255)              gACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCaaa                                                                                                               >  1:724226/1‑140 (MQ=255)               acacCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAAc                                                                                                              >  1:326543/1‑140 (MQ=255)                  ccccATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGaa                                                                                                           >  1:285905/1‑140 (MQ=255)                   cccATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGaat                                                                                                          >  2:176916/1‑140 (MQ=255)                   cccATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGaat                                                                                                          >  1:409243/1‑140 (MQ=255)                   cccATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGaat                                                                                                          >  2:606679/1‑140 (MQ=255)                   cccATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGaat                                                                                                          >  1:149912/1‑140 (MQ=255)                      aTCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGaataaa                                                                                                       <  2:1987/140‑1 (MQ=255)                      aTCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGaataaa                                                                                                       <  1:116107/140‑1 (MQ=255)                          tATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAGTAGc                                                                                                   <  2:149912/140‑1 (MQ=255)                                  gATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAAttt                                                                                           >  1:19966/1‑140 (MQ=255)                                          tGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCAtt                                                                                   >  2:208793/1‑140 (MQ=255)                                                 ttCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAAt                                                                                             >  1:558435/1‑123 (MQ=255)                                                 ttCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAAt                                                                                             <  2:558435/123‑1 (MQ=255)                                                  tCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCac                                                                           >  2:438629/1‑140 (MQ=255)                                                   cACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCacg                                                                          <  2:409243/140‑1 (MQ=255)                                                   cACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCacg                                                                          >  1:7650/1‑140 (MQ=255)                                                   cACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCacg                                                                          >  1:798927/1‑140 (MQ=255)                                                    aCGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCacga                                                                         <  1:176916/140‑1 (MQ=255)                                                    aCGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCacga                                                                         <  2:476717/140‑1 (MQ=255)                                                    aCGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCacga                                                                         <  2:326543/140‑1 (MQ=255)                                                    aCGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCacga                                                                         <  2:496365/140‑1 (MQ=255)                                                      gTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCacgacg                                                                       >  2:432754/1‑140 (MQ=255)                                                      gTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCacgacg                                                                       >  1:308302/1‑140 (MQ=255)                                                      gTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCacgacg                                                                       <  1:389645/140‑1 (MQ=255)                                                      gTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCacgacg                                                                       >  2:444718/1‑140 (MQ=255)                                                        cGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCt                                                                     >  1:85920/1‑140 (MQ=255)                                                        cGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCt                                                                     >  1:40144/1‑140 (MQ=255)                                                             tGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTaaa                                                                <  2:206608/140‑1 (MQ=255)                                                             tGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTaaa                                                                <  2:316517/140‑1 (MQ=255)                                                                  cgAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGc                                                           >  2:106563/1‑140 (MQ=255)                                                                   gAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCg                                                          <  2:348255/140‑1 (MQ=255)                                                                   gAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCg                                                          <  2:285905/140‑1 (MQ=255)                                                                    aaTATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGt                                                         <  2:305547/140‑1 (MQ=255)                                                                    aaTATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGACAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGt                                                         <  1:420289/140‑1 (MQ=255)                                                                               gctttgctttTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGaa                                              >  2:594948/1‑140 (MQ=255)                                                                               gctttgctttTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGaa                                              >  2:672792/1‑140 (MQ=255)                                                                               gctttgctttTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGaa                                              <  2:235189/140‑1 (MQ=255)                                                                                   tgctttTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGa                                               >  1:423143/1‑135 (MQ=255)                                                                                   tgctttTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGa                                               <  2:423143/135‑1 (MQ=255)                                                                                               tcTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGcc                              <  2:230655/140‑1 (MQ=255)                                                                                               tcTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGcc                              <  1:672792/140‑1 (MQ=255)                                                                                               tcTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGcc                              <  2:798927/140‑1 (MQ=255)                                                                                                cTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGccc                             <  1:438629/140‑1 (MQ=255)                                                                                                     ataAAGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGt                        <  1:444718/140‑1 (MQ=255)                                                                                                         aGCCACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAaca                    >  2:293129/1‑140 (MQ=255)                                                                                                           ccACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACAcc                  >  1:81969/1‑140 (MQ=255)                                                                                                            cACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCa                 <  2:249527/140‑1 (MQ=255)                                                                                                            cACGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCa                 <  2:620553/140‑1 (MQ=255)                                                                                                             aCGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCaa                <  2:85920/140‑1 (MQ=255)                                                                                                              cGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCaaa               >  1:205033/1‑140 (MQ=255)                                                                                                              cGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCaaa               >  1:1526/1‑140 (MQ=255)                                                                                                              cGTAAAAGCGGTGACGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCaaa               >  1:422486/1‑140 (MQ=255)                                                                                                                           aCGTCAAATGCTGAGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGa  >  1:178341/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                               TCATATTACCCCAGACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGA  >  NC_000913/4153095‑4153356

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