Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,243,428	T→C	E438E (GAA→GAG)	uxaC ←	uronate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,243,428	0	T	C	100.0%	175.8 / NA	50	E438E (GAA→GAG)	uxaC	uronate isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/15);  total (35/15)

GTTATCCAGCGCATGGATCTTGATGTATTGCATATCAACCCCAGACCTTAGTTCAGTTCAATGGCGAAGTAGTCACGCGCATTGTTAAAGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCTTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGT  >  NC_000913/3243282‑3243569                                                                                                                                                   |                                                                                                                                              gTTATCCAGCGCATGGATCTTGATGTATTGCATATCAACCCCAGACCTTAGTTCAGTTCAATGGCGAAGTAGTCACGCGCATTGTTAAAGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTcc                                                                                                                                            >  2:353830/1‑150 (MQ=255)  ttATCCAGCGCATGGATCTTGATGTATTGCATATCAACCCCAGACCTTAGTTCAGTTCAATGGCGAAGTAGTCACGCGCATTGTTAAAGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCAc                                                                                                                                          >  1:153977/1‑151 (MQ=255)              tGGATCTTGTTGTATTGCATACCACCCCCAGACCTTAGTTCAGTTCAATGGCGAAGTAGTCACGCGCTTTGTTAAAGCAAATATTTTTCCCCACCCCGCCCAGCGGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGa                                                                                                                               <  1:274299/150‑1 (MQ=255)               gggtCTTGATATTTTGCATTTCACCCCCAGACCTTATTTCAGTTCAATGGCGAATTAGTCACGCGCATTGTTAAAGCAAATATTTTTCACCATCCCCCCCAGCATGTTCGTCTCCGCCGGTTCTTCGCCCGCCTCCTCCCAGCGACCGatc                                                                                                                             <  2:514014/148‑1 (MQ=255)                           aTTGCATATCAACCCCAGACCTTAGTTCAGTTCAATGGCGAAGTAGTCACGCGCATTGTTAAAGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGcac                                                                                                                    >  1:298789/1‑148 (MQ=255)                                atatCAACCCCAGACCTTAGTTCAGTTCAATGGCGAAGTAGTCACGCGCATTGTTAAAGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATgcg                                                                                                            >  1:305104/1‑151 (MQ=255)                                     acccccAGACCTTAGTTCAGTTCAAGGGCGAAGTAGTCACGCGCATTGTTAAAGAAAATATTTTCCCCCACCCCGCCCAGCAGGTTGATGCCCGCCGGTGCTCCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGcggcgg                                                                                                        <  1:70570/148‑1 (MQ=255)                                       ccccAGACCTTAGTTCAGTACAAGGGCGAAGTAGTCACGCGCATTGTTAAAGCAAATATTTTCCCCCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATgcgg                                                                                                           <  1:82646/145‑1 (MQ=255)                                         ccaacctttAGTTCAGTACAACGGCAAAATAGCCCGGCGCATTTTTAAAGGAAATATTTTACACCGTCTCGCCCAGCAGGTTGATGGCCACGGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAATTa                                                                                                   <  1:830867/144‑1 (MQ=255)                                            gACCTTAGTTCAGTTCAATGGCGAAGTAGTCACGCGCATTGTTAAAGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTc                                                                                                >  2:824355/1‑151 (MQ=255)                                              ccTTAGTTCAGTTCAATGGCGAAGTAGTCACGCGCATTGTTAAAGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGt                                                                                              >  1:385343/1‑151 (MQ=255)                                                     tcagttcaATGGCGAAGTAGTCACGCGCATTGTTAAAGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGgt                                                                                       >  1:52702/1‑151 (MQ=255)                                                          tcaATGGCGAAGTAGTCACGCGCATTGTTAAAGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATg                                                                                  >  2:301153/1‑151 (MQ=255)                                                            aaTGGCGAAGTAGTCACGCGCATTGTTAAAGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGAc                                                                                >  1:555452/1‑151 (MQ=255)                                                            aaTGGCGAAGTAGTCACGCGCATTGTAAAAGAAATTTTTTTTCCCCATCGCGCCCAGCAGGTTGTTGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGAc                                                                                <  2:66965/151‑1 (MQ=255)                                                             aTGGCGAAGTAGTCACGCGCATTGTTAAAGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGCAGTATTCGTGACGGGAGTATGa                                                                                 >  1:809033/1‑149 (MQ=255)                                                                   aaGTAGTCACGCGCATTGTTAAAGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGc                                                                         >  1:610778/1‑151 (MQ=255)                                                                   aaGTAGTCACGCGCATTGTTAAAGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGc                                                                         >  2:599274/1‑151 (MQ=255)                                                                   aaGTAGTCACGCGAATTGTTAAAGAAAATATTTTTCCCCACCTCGCCCAGCAGGTTGATGCCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGc                                                                         <  1:794499/151‑1 (MQ=255)                                                                      tagtCACGCGCATTGTTAAAGCAAATATTTTTCACCATCTCGCCCAGCGGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTAc                                                                      >  2:746800/1‑151 (MQ=255)                                                                              cgcATTGTTAAAGCAAATATTTTTCCCCATCTCCCCCAGCGGGTTGATGCCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTgtc                                                               <  2:744379/150‑1 (MQ=255)                                                                               gcATTGTAAAAGCAAATATTTTCCCCCCTCCCCCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACgg                                                                    <  1:402693/144‑1 (MQ=255)                                                                                cATTGTTAAAGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATTATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTgtcagt                                                            >  2:816351/1‑151 (MQ=255)                                                                                cATTGTTAAAGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTgtcagt                                                            >  1:20932/1‑151 (MQ=255)                                                                                   tGTTAAAGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTgtcagt                                                            >  1:737708/1‑148 (MQ=255)                                                                                     ttAAAGAAAATATTTTTCCCCACCCCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTgtcagtca                                                          <  1:475772/148‑1 (MQ=255)                                                                                       aaaGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCAt                                                      >  1:594473/1‑150 (MQ=255)                                                                                        aaGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATAc                                                    >  2:205928/1‑151 (MQ=255)                                                                                          gCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCg                                                  >  2:31380/1‑151 (MQ=255)                                                                                            aaaTATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGAc                                                >  1:337024/1‑151 (MQ=255)                                                                                             aaTATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCCTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGAc                                                >  2:92814/1‑150 (MQ=255)                                                                                               tatTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACaa                                              >  1:21418/1‑150 (MQ=255)                                                                                                  ttttCACCATCCCGCCCAGCAGGTTGATGCCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGc                                           <  2:156636/150‑1 (MQ=255)                                                                                                   tttCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGCgg                                         >  2:386825/1‑151 (MQ=255)                                                                                                       aCCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGCGGCTc                                      >  1:216368/1‑150 (MQ=255)                                                                                                       aCCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGCGGCTc                                      >  1:380873/1‑150 (MQ=255)                                                                                                        ccATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCGGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGCGGCTc                                      >  2:43927/1‑149 (MQ=255)                                                                                                        ccACCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGCGGCTcag                                    <  2:790474/151‑1 (MQ=255)                                                                                                         cATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGCGGCTcag                                    >  1:53743/1‑150 (MQ=255)                                                                                                                cccAGCAGGTTGATGCCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGCGGCTCAGCAGACCg                             <  1:659776/150‑1 (MQ=255)                                                                                                                  cagcagGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGCGGCTCAGCAGACCGAGc                          >  1:94148/1‑151 (MQ=255)                                                                                                                     cagGTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGCGGCTCAGCAGACCGAGCTgc                       >  1:550146/1‑151 (MQ=255)                                                                                                                       ggTTGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGCGGCTCAGCAGACCGAGCTgcgc                     >  1:670261/1‑151 (MQ=255)                                                                                                                         ttGATGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGcc                    >  1:475506/1‑150 (MQ=255)                                                                                                                             tGTCCGCCGGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTg               >  1:29891/1‑151 (MQ=255)                                                                                                                                ccgccgGTGCTTCGCCCGCCTCCCCCCAGCGACCGATCATCTGGCACAGAATGCGGCGGGAGTATTCGTTACGGGGGTTTGACAGGAAGCTACGGCAGTCAGTCAGCATACCGCCAAAGCGGCCCAGCAGGCCGAGCTGCGCCAGTTGGGt            >  1:818615/1‑151 (MQ=255)                                                                                                                                    cgGTGCTTCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATc        >  1:650813/1‑151 (MQ=255)                                                                                                                                          ttCGCCCGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGt  >  1:106427/1‑151 (MQ=255)                                                                                                                                            cccccGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGa     <  2:625636/144‑1 (MQ=255)                                                                                                                                              cccGCCTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACg   <  1:329284/146‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                              GTTATCCAGCGCATGGATCTTGATGTATTGCATATCAACCCCAGACCTTAGTTCAGTTCAATGGCGAAGTAGTCACGCGCATTGTTAAAGCAAATATTTTTCACCATCTCGCCCAGCAGGTTGATGTCCGCCGGTGCTTCGCCCGCTTCCACCCAGCGACCGATCATCTGGCACAGAATGCGGCGGAAGTATTCGTGACGGGTGTATGACAGGAAGCTACGGCTGTCAGTCAGCATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGT  >  NC_000913/3243282‑3243569

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑