Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	718,705	G→T	R586S (CGC→AGC)	speF ←	ornithine decarboxylase, degradative

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	718,705	0	G	T	100.0%	155.8 / NA	45	R586S (CGC→AGC)	speF	ornithine decarboxylase, degradative
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (18/27);  total (18/27)

TTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGA  >  NC_000913/718574‑718839                                                                                                                                    |                                                                                                                                       ttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGTAAAAAAAATGAATTTTTTTACCAAAATATCCCAGGGGAAAAACGGGGTAGCTCCCTTCAt                                                                                                                                >  2:292052/1‑140 (MQ=255) ttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCAt                                                                                                                                <  1:1091894/140‑1 (MQ=255)  tGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATg                                                                                                                               <  1:1176489/140‑1 (MQ=255)   gCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGc                                                                                                                              <  2:871018/140‑1 (MQ=255)   gCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGc                                                                                                                              <  2:369406/140‑1 (MQ=255)      gATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGt                                                                                                                           <  2:539932/140‑1 (MQ=255)         tCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTg                                                                                                                        <  1:1076747/140‑1 (MQ=255)          cATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGAGCCAAAAAATGCATTTTCTTTCACAACCCATGCCGGGGGTAAAACGCGGAGGTCTCTTCTTGCTGttttt                                                                                                                       >  2:103753/1‑138 (MQ=255)                gACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATgg                                                                                                                 <  2:746499/140‑1 (MQ=255)                       gggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAg                                                                                                          <  1:967885/140‑1 (MQ=255)                       gggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAg                                                                                                          <  1:814547/140‑1 (MQ=255)                        ggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGc                                                                                                         <  1:292052/140‑1 (MQ=255)                         gAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCa                                                                                                        <  2:565240/140‑1 (MQ=255)                           aGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACt                                                                                                      <  2:748719/140‑1 (MQ=255)                             tGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAACATGCATTTTCTGACACAACAGACGCAGGGTATAAACGGCGTAGCTCTCTTCCTTCTGTTTTTTGATGGGAGGTAGCACCTc                                                                                                    >  2:459211/1‑140 (MQ=255)                               cTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTACACAACTGACGCAGGGGATAAACGGCGGAGGTCTCTTCCTTCTGTTTTTTGATGGAAGGGAGGACTTCtt                                                                                                  >  2:814547/1‑139 (MQ=255)                                 ccTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCt                                                                                                <  1:45754/140‑1 (MQ=255)                                      cGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATgg                                                                                           <  1:1044348/140‑1 (MQ=255)                                       ggAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACATGCATTTTCTTGCACAACCGAACCAGGGGATAAACCGCGGCGCTCTCTTCTTTCTTTTTTTTGATGGAAGGGAGCACCTCTGCTAAttgg                                                                                          >  2:714353/1‑137 (MQ=255)                                       ggAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATTCATTTTCTTACACAAATGACGCAGGGGATAAACGGCGTAGCTCTCTTTATTCTGTTTTTTGATGGAAGGGAGGACTTCTTCTAAttgg                                                                                          >  2:1244881/1‑137 (MQ=255)                                             tctcTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGTTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCg                                                                                    <  2:665874/140‑1 (MQ=255)                                                          ttgttTCACGTTGTGGCGGGCATAAAAAACATGCATTTTCTGACACAAACGACGCAGGGGGTAAACCGCGGAGCTCTCTTCCTTCTGTTTTTTGATGGGAGGGAGGACTTCTGCTAATGGCGCGTCGGGCTCCAGGAgat                                                                       >  2:724365/1‑138 (MQ=255)                                                            gttTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGtttt                                                                     <  2:455708/140‑1 (MQ=255)                                                                cACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGa                                                                 <  1:724365/140‑1 (MQ=255)                                                                cACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGa                                                                 <  1:1125418/140‑1 (MQ=255)                                                                                    aaTCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGcaaca                                             <  2:924618/140‑1 (MQ=255)                                                                                        aTGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGc                                                 >  1:810001/1‑132 (MQ=255)                                                                                        aTGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGc                                                 <  2:810001/132‑1 (MQ=255)                                                                                        aTGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAACGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGtt                                         >  1:1179290/1‑140 (MQ=255)                                                                                         tGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGcaaca                                             <  2:1097734/135‑1 (MQ=255)                                                                                         tGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGcaaca                                             >  1:1097734/1‑135 (MQ=255)                                                                                         tGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTg                                        >  1:491571/1‑140 (MQ=255)                                                                                               tCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAg                                  <  1:278558/140‑1 (MQ=255)                                                                                                   gACACAACTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTg                              <  1:1232324/140‑1 (MQ=255)                                                                                                        aaCTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCat                         >  2:461973/1‑140 (MQ=255)                                                                                                          cTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCatat                       >  2:544250/1‑140 (MQ=255)                                                                                                          cTGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCatat                       >  2:1048206/1‑140 (MQ=255)                                                                                                           tGACGCAGGGTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATAtc                      <  2:198056/140‑1 (MQ=255)                                                                                                                   ggaataACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCg              >  1:364933/4‑140 (MQ=255)                                                                                                                   ggTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCg              <  1:1271378/140‑1 (MQ=255)                                                                                                                    gTATAACCGGCGTAGCTCTCTTCATGCTGTTTGTCGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTCCCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCgg             >  2:305570/1‑140 (MQ=255)                                                                                                                     tataACCGGCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGa            >  2:429539/1‑140 (MQ=255)                                                                                                                             gCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTcagcag    >  2:749772/1‑140 (MQ=255)                                                                                                                             gCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTcagcag    >  2:1183826/1‑140 (MQ=255)                                                                                                                             gCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTcagcag    >  1:1250192/1‑140 (MQ=255)                                                                                                                             gCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTcagcag    >  1:225582/1‑140 (MQ=255)                                                                                                                             gCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTc         <  2:1063764/135‑1 (MQ=255)                                                                                                                             gCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTc         >  1:1063764/1‑135 (MQ=255)                                                                                                                             gCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTCCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGt          >  2:1058766/1‑134 (MQ=255)                                                                                                                             gCGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTCCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGt          <  1:1058766/134‑1 (MQ=255)                                                                                                                              cGTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAgg   >  1:446289/1‑140 (MQ=255)                                                                                                                               gTAGCTCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGa  >  2:1224726/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                       TTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGA  >  NC_000913/718574‑718839

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