Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,827,031	T→A	D410V (GAT→GTT)	astB ←	succinylarginine dihydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,827,031	0	T	A	100.0%	187.4 / NA	54	D410V (GAT→GTT)	astB	succinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (34/20);  total (34/20)

GCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCG  >  NC_000913/1826937‑1827124                                                                                               |                                                                                              gcgcTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgca                                                                                          <  1:604150/100‑1 (MQ=255) gcgcTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgca                                                                                          <  1:1326621/100‑1 (MQ=255)     tGGAACGGATAAACCGAACCGATATTCAGTAATTGGCTCAATACATCTGGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgcagc                                                                                        <  2:1290094/98‑1 (MQ=255)     tGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGcagcag                                                                                      <  1:4522/100‑1 (MQ=255)     tGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGcagcag                                                                                      <  1:1092715/100‑1 (MQ=255)      ggAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTAGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGt                                                                                     >  2:1438231/1‑100 (MQ=255)          cGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAgg                                                                                 >  2:759948/1‑100 (MQ=255)                aaccgaaccgAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATcg                                                                           >  2:1426155/1‑100 (MQ=255)                aaccgaaccgAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATcg                                                                           >  1:725067/1‑100 (MQ=255)                 accgaaccgaGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATcgc                                                                          >  1:247962/1‑100 (MQ=255)                  ccgaaccgaGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATcgcg                                                                         >  1:1284016/1‑100 (MQ=255)                         gagaTTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTaa                                                                  >  1:1417030/1‑100 (MQ=255)                           gaTTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACg                                                                <  1:716598/100‑1 (MQ=255)                             ttCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGAt                                                              >  1:164965/1‑100 (MQ=255)                                 gTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCAc                                                          >  2:566251/1‑100 (MQ=255)                                  tAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCAcc                                                         >  1:1083848/1‑100 (MQ=255)                                   aaTTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCAccc                                                        >  2:1337731/1‑100 (MQ=255)                                     ttGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAg                                                      <  1:452628/100‑1 (MQ=255)                                          tCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAACGCGATCGCGGTAGTAACGATCCACCCAGTCAtt                                                 <  2:946813/100‑1 (MQ=255)                                             aTACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAg                                              >  1:597501/1‑100 (MQ=255)                                              tACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAgc                                             >  1:497610/1‑100 (MQ=255)                                                cATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAgcgc                                           >  2:403294/1‑100 (MQ=255)                                                  tCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAt                                         >  1:1349313/1‑100 (MQ=255)                                                  tCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAt                                         >  1:1523641/1‑100 (MQ=255)                                                   ccAGTGCTTCCCGCCCTTCGCGCAGCTATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAtt                                        >  1:954188/1‑100 (MQ=255)                                                   ccAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCGGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAtt                                        >  1:1410826/1‑100 (MQ=255)                                                   ccAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAtt                                        >  1:153924/1‑100 (MQ=255)                                                    cAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTa                                       <  2:876334/100‑1 (MQ=255)                                                    cAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTa                                       <  1:1507208/100‑1 (MQ=255)                                                      gTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTaaa                                     >  1:1520950/1‑100 (MQ=255)                                                        gCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACa                                   >  2:1091522/1‑100 (MQ=255)                                                        gCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACa                                   >  2:676674/1‑100 (MQ=255)                                                        gCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACa                                   >  2:1234746/1‑100 (MQ=255)                                                         cTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAg                                  >  2:464243/1‑100 (MQ=255)                                                                 ccTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGt                          >  2:906337/1‑100 (MQ=255)                                                                 ccTTCGCGCAGCAATTGCGGGTCGGCCAGAACCGCACCAGTAAGGAGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGt                          <  1:670568/100‑1 (MQ=255)                                                                  cTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGtt                         >  2:1380614/1‑100 (MQ=255)                                                                    tCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtca                       <  2:1083841/100‑1 (MQ=255)                                                                    tCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtca                       <  1:566254/100‑1 (MQ=255)                                                                    tCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtca                       <  2:725064/100‑1 (MQ=255)                                                                      gcgcAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatc                     <  2:1204340/100‑1 (MQ=255)                                                                       cgcAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtca                       <  2:121857/97‑1 (MQ=255)                                                                           gcaATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAc                >  1:903771/1‑100 (MQ=255)                                                                           gcaATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAc                <  2:754040/100‑1 (MQ=255)                                                                           gcaATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAc                >  2:221835/1‑100 (MQ=255)                                                                             aaTTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAccg              >  1:395624/1‑100 (MQ=255)                                                                                  cGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGa         >  1:1026267/1‑100 (MQ=255)                                                                                   gggTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGAt        >  1:1054561/1‑100 (MQ=255)                                                                                   gggTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGAt        <  1:172906/100‑1 (MQ=255)                                                                                    ggTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGAtt       <  2:633740/100‑1 (MQ=255)                                                                                    ggTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGAtt       <  1:1069873/100‑1 (MQ=255)                                                                                       cGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCAc    >  1:1084476/1‑100 (MQ=255)                                                                                       cGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCAc    >  1:8343/1‑100 (MQ=255)                                                                                         gCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCg  >  1:77059/1‑100 (MQ=255)                                                                                               |                                                                                              GCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCG  >  NC_000913/1826937‑1827124

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