Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,108,382	G→A	R258R (CGC→CGT)	speC ←	ornithine decarboxylase, biosynthetic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,108,382	0	G	A	100.0%	102.2 / NA	32	R258R (CGC→CGT)	speC	ornithine decarboxylase, biosynthetic
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/10);  total (22/10)

CACACAGATGCCCAACGGTATCGATCACCTGACGGGCGTTATAGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGGCGCAGATACTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTCCAGATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATG  >  NC_000913/3108246‑3108510                                                                                                                                         |                                                                                                                                  cacacaGATGCCCAACGGTATCGATCACCTGACGGGCGTTATAGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGacgc                                                                                                                                <  2:396487/140‑1 (MQ=255)   cacaGATGCCCAACGGTATCGATCACCTGACGGGCGTTATAGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGGGCAACGTCGCGAATTTGCTGACGCAg                                                                                                                              >  2:673938/1‑140 (MQ=255)   cacaGATGCCCAACGGTATCGATCACCTGACGGGCGTTATAGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAg                                                                                                                              >  1:521215/1‑140 (MQ=255)    acaGATGCCCAACGGTATCGATCACCTGACGGGCGTTATAGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGa                                                                                                                             >  1:519156/1‑140 (MQ=255)    acaGATGCCCAACGGGATCGATCACGTGAGGGGCGTTTTAGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGa                                                                                                                             <  1:259785/140‑1 (MQ=255)        aTGCCCAACGGTATCGATCACCTGACGGGCGTTATAGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATAct                                                                                                                         >  1:711277/1‑140 (MQ=255)         tGCCCAACGGTATCGATTACCTCACGGGCGTTATACACAGTGCCGTCAAAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATActc                                                                                                                        <  2:318263/140‑1 (MQ=255)         tGCCCAACGGTATCGATCACCTGACGGGCGTTATAGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATActc                                                                                                                        <  2:380319/140‑1 (MQ=255)         tGCCCAACGGTATCGATCACCTGACGGGCGTTATAGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATActc                                                                                                                        <  2:431433/140‑1 (MQ=255)         tGCCCAACGGTATCGATCACCTGACGGGCGTTATAGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATActc                                                                                                                        <  1:563086/140‑1 (MQ=255)           cccAACGGTATCGATCACCTGACGGGCGTTATAGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCtt                                                                                                                      <  2:488558/140‑1 (MQ=255)              aaCGGTATCGATCACCTGACGGGCGTTATAGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGTGCAACGTAGCGAATTTGCTGACGCAGATACTCTTCAt                                                                                                                   >  1:297052/1‑140 (MQ=255)                           aCCTGACGGGCGTTATAGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGTgcgc                                                                                                      >  1:477823/1‑140 (MQ=255)                            ccTGACGGGCGTTATAGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGTGCGCa                                                                                                     <  2:671016/140‑1 (MQ=255)                             cTGACGGGCGTTATAGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGTGCGCAt                                                                                                    >  1:171444/1‑140 (MQ=255)                                    gCGTTATAGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGTGCGCATCAATA‑cc                                                                                            <  2:477823/140‑1 (MQ=255)                                     cGTTATAGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGTGCGCATCAATA‑CCg                                                                                           >  1:562535/1‑140 (MQ=255)                                          taGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCTTTAGAATAGTGCGGATCAATA‑CCGCCAAt                                                                                      >  1:248537/1‑140 (MQ=255)                                            gACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGTTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATCCTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGa                                                                                    <  2:519156/140‑1 (MQ=255)                                                gTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGAAAcc                                                                                >  1:362576/1‑140 (MQ=255)                                                 tGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGAAACCa                                                                               >  2:735102/1‑140 (MQ=255)                                                 tGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGAAGATACTCGTCAATTAAACAGTGCGGATCAAAA‑CCGCCAATGAAACCa                                                                               >  2:351387/1‑140 (MQ=255)                                                      tCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGGCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGACCGCATCAATA‑CCAGCAATGAAACCAAACgg                                                                          >  2:337114/1‑140 (MQ=255)                                                           ggTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGAAACCAAACGGGTTGc                                                                     >  2:677078/1‑140 (MQ=255)                                                               cccAGCGGAAAAATCGCCGGGCGATGGGGGCGCGGCAGGTCGGCTTTTTGTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCATAAAAACAGTGCGCATCAATA‑CCGCCAATGAAACCAAACGGGTTGCGTGa                                                                 <  2:551062/140‑1 (MQ=255)                                                                ccAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGAAACCAAACGGGTTGCGTGaa                                                                >  1:551062/1‑140 (MQ=255)                                                                  aGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCATTAAGACAGTGCGCATCGATA‑GCGCCAATGAAACCAAACGGGTTGCGTGAAGc                                                              >  2:291069/1‑140 (MQ=255)                                                                   gCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCt                                                             >  1:394203/1‑140 (MQ=255)                                                                     tGAAAAATCGCCAGGCGATCCGGGCGGGGCAGGCCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATCCTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTc                                                           <  2:89743/140‑1 (MQ=255)                                                                      gAATAATCGCCAGGCGATACGGTCGGGGCCGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATCCTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTcc                                                          <  1:295479/140‑1 (MQ=255)                                                                            tCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTCCAGATAg                                                    >  1:145005/1‑140 (MQ=255)                                                                            tCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTCCAGATAg                                                    >  1:678784/1‑140 (MQ=255)                                                                              gCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTc                                                           <  2:675282/131‑1 (MQ=255)                                                                              gCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTc                                                           >  1:675282/1‑131 (MQ=255)                                                                              gCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTCCAGATAGAc                                                  >  1:89743/1‑140 (MQ=255)                                                                                        aggggCGGGGCAGGTCGGTTTTCTCGGGCGCCACGCCGCAAATTTGCGACCGCAGATCCTTTTCATTAAAACAGTGCGCATAAATA‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTCCAGATAGACCGGCGTCGcc                                        <  2:248537/138‑1 (MQ=255)                                                                                         ggggCGGGGCAGGTCGGCTTTTTCGGGCGCAACGTCGCAAATTTGCTGACGAAGATCCTCTTCATAAAACCGGTGCGCACAAATC‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTCCAGATAGACCGGCGTCGccc                                       <  2:655637/139‑1 (MQ=255)                                                                                         ggggCGCGGCAGGTCGGCTTTTTCGGGCGCAACGTCGCAAATTTGCTGACGCAGATCCTCTTCATAAAAACAGGGCGCAACAATA‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTCCAGATAGACCGGCGTCGccc                                       <  2:16548/139‑1 (MQ=255)                                                                                             cgcgGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGTTGACGCAGATACTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGAAACCAAACGGGCTGCGTGAAGCTACCTGATAGATCGGCGTCGCCCCCGc                                   >  1:338480/1‑140 (MQ=255)                                                                                                       gggttttttCGGGCCCATCGTCTGGAATTTGCTGACGAATATACTCTTCATAAAAACAGTGCCCACAAATA‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTCCAGATAGACCGGCGTCGCCCCCGCCTGAATCAgc                         <  2:297052/136‑1 (MQ=255)                                                                                                          ttttttCGGGCGCAAGCCCGCAAATTTTTTTGCGAGGATCCTCTTAATCAAAACAGTGCGCA‑AAATACCCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTCCAGATAGACCGGCGTCGCCCCCGCCTGAATCAgcgcg                      <  2:338480/139‑1 (MQ=255)                                                                                                               tCTGGCGCAACGTCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGAAAACAAACGGGTTGCGTGAAGCTTCCAAATAAACCCGCGTCGCCCCCGCCTGAATCAGCGCGCCGtg                 >  2:240713/1‑140 (MQ=255)                                                                                                                    cgcAACGTCGCAAATGTGCTGCTGCAGATCCTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTCCAGATAGACCGGCGTCG‑CCCCGCCTGAATCAGCGCGCCGTGATGAtt           <  1:291069/140‑1 (MQ=255)                                                                                                                     gcAACGTCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGTGCGCATAAATA‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTCCAGATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGAtt           <  2:584771/140‑1 (MQ=255)                                                                                                                          ggcgcaAATTTGCGGACGCAGATCCTCTTCATAAAAACAGTGCGCACAAAGC‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTCCAGATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACt      <  2:477239/138‑1 (MQ=255)                                                                                                                          gccgcaAATTGGCTGACACAGATACCTTTTATAAAAAAAGGGCGCATAAATC‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTCCAGATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACt      <  2:70194/138‑1 (MQ=255)                                                                                                                          gTCGCGATTTTGCTGACGCAGATCCTTTTCATAAAAACAGTGCGCATCAATC‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTCCAGATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACt      <  1:673938/140‑1 (MQ=255)                                                                                                                           tCGCGAATTTGCTGACGCAGATACTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTCCAGATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACtt     <  2:394203/140‑1 (MQ=255)                                                                                                                              caattttGCTGACGCAGATCCTCTTCATAAAAACAGTCCCCATCAATC‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTCCAGATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATg  <  2:93419/136‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  CACACAGATGCCCAACGGTATCGATCACCTGACGGGCGTTATAGACAGTGCCGTCATAGGTTCCCAGCTGAATAATCGCCAGGCGATACGGGCGCGGCAGGTCGGCTTTTTCTGGCGCAACGTCGCGAATTTGCTGGCGCAGATACTCTTCATTAAAACAGTGCGCATCAATA‑CCGCCAATGAAACCAAACGGGTTGCGTGAAGCTTCCAGATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATG  >  NC_000913/3108246‑3108510

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