Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,574,062	A→T	C273S (TGC→AGC)	asd ←	aspartate‑semialdehyde dehydrogenase, NAD(P)‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,574,062	0	A	T	100.0%	193.7 / NA	55	C273S (TGC→AGC)	asd	aspartate‑semialdehyde dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (27/28);  total (27/28)

CGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCACTC  >  NC_000913/3573966‑3574155                                                                                                 |                                                                                              cGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTgcg                                                                                            <  2:1780400/100‑1 (MQ=255)    cTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTgcgc                                                                                           <  2:987127/98‑1 (MQ=255)      ttCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATg                                                                                       <  1:1772689/100‑1 (MQ=255)          cccACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGcccc                                                                                   >  2:1352282/1‑100 (MQ=255)          cccACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAATGAATACCTGGCTATGGCTGCGCAATGcccc                                                                                   >  1:210185/1‑100 (MQ=255)           ccACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCg                                                                                  <  2:1850166/100‑1 (MQ=255)           ccACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCg                                                                                  <  1:2149851/100‑1 (MQ=255)             aCGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGac                                                                                >  2:1830179/1‑100 (MQ=255)                   tgtgCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGcac                                                                          <  2:485438/100‑1 (MQ=255)                    gtgCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGcaca                                                                         <  2:104046/100‑1 (MQ=255)                    gtgCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGcaca                                                                         >  1:2225590/1‑100 (MQ=255)                    gtgCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGcaca                                                                         >  1:817253/1‑100 (MQ=255)                    gtgCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGcaca                                                                         <  2:542349/100‑1 (MQ=255)                     tgCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGcacac                                                                        <  2:327633/100‑1 (MQ=255)                     tgCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGcacac                                                                        <  2:669680/100‑1 (MQ=255)                           gccagcAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAAc                                                                  >  2:2102109/1‑100 (MQ=255)                             cagcagTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCa                                                                <  2:2106696/100‑1 (MQ=255)                              agcagTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCAt                                                               >  2:1319468/1‑100 (MQ=255)                                cagTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATc                                                              <  2:2129514/99‑1 (MQ=255)                                  gatcttcCACGGTAGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTAc                                                           >  1:1581749/3‑100 (MQ=255)                                   ttcttcCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTAcc                                                          >  2:1698028/1‑100 (MQ=255)                                    tcttcCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCg                                                         <  1:1420139/100‑1 (MQ=255)                                     cttcCACGGTCGGAATAGACACATCTTTTTTTTTTTTTTTTTTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCg                                                         >  2:552200/1‑99 (MQ=255)                                     cttcCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCgg                                                        >  1:2135681/1‑100 (MQ=255)                                       tcCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGaa                                                      >  2:1273639/1‑100 (MQ=255)                                         cACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAAtt                                                    >  1:492504/1‑100 (MQ=255)                                             gTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACgg                                                <  2:846948/100‑1 (MQ=255)                                             gTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATATACCGGAATTAAgg                                                >  1:2104877/1‑100 (MQ=255)                                              tCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGa                                               <  1:1460968/100‑1 (MQ=255)                                                ggAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAg                                             <  2:2225582/100‑1 (MQ=255)                                                ggAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAg                                             <  2:860752/100‑1 (MQ=255)                                                  aaTAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGAt                                           <  2:1581744/100‑1 (MQ=255)                                                   aTAGACACATCTTTTTTTTTTTTATTTGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATGACGGAAGAtg                                          >  2:1529382/1‑100 (MQ=255)                                                   aTAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGAtg                                          >  1:1242904/1‑100 (MQ=255)                                                   aTAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGAtg                                          >  1:1242905/1‑100 (MQ=255)                                                     aGACACATCTTTTTTTTTTTTATTTTTGAATGCCTGGCTGTGGCTACGCAATGCCCCGACACGCACACATAAATCATCTACCGGAATTACGGAAGAtgt                                         >  2:2026329/1‑99 (MQ=255)                                                       acacATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGtt                                      <  1:552201/100‑1 (MQ=255)                                                        cacaTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTg                                     <  2:439685/100‑1 (MQ=255)                                                            tCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTTCCGTAATTCCGGTAGAtggtttatgtg                                 <  1:2079540/100‑10 (MQ=25)                                                             cTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGAt                                >  1:1375011/1‑100 (MQ=255)                                                             cTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGAt                                >  2:629825/1‑100 (MQ=255)                                                               tttttCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCt                              <  2:451399/100‑1 (MQ=255)                                                                ttttCAATTTAATAGTGAATGCCGGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAGCCACCTACCGGAATGACGGAAGATGTGTTGAGGATCt                              <  2:1983029/99‑1 (MQ=255)                                                                 tttCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttg                            <  1:1830184/100‑1 (MQ=255)                                                                   tCAATTTAATAGTGAATGTCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttgtt                          <  1:2026336/100‑1 (MQ=255)                                                                   tCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttgtt                          >  2:747558/1‑100 (MQ=255)                                                                     aaTTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTgg                        <  1:2102116/100‑1 (MQ=255)                                                                      aTTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATG‑CCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGtt                      <  2:803671/100‑1 (MQ=255)                                                                          aaTAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTcc                   >  2:1849954/1‑100 (MQ=255)                                                                           aTAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCg                  <  2:566905/100‑1 (MQ=255)                                                                                gAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGc             >  2:2088673/1‑100 (MQ=255)                                                                                     ccTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCttt        >  1:1072342/1‑100 (MQ=255)                                                                                     ccTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTCCCGCCTGCCCttt        >  2:251671/1‑100 (MQ=255)                                                                                       tggctgtggctgCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTcc      <  1:1849959/100‑1 (MQ=255)                                                                                       tggctgtggctgCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTcc      <  1:1529386/100‑1 (MQ=255)                                                                                       tggctgtggctgCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTcc      <  1:2191383/100‑1 (MQ=255)                                                                                           tgtgGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCActc  >  1:928230/1‑100 (MQ=255)                                                                                           tgtgGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCActc  >  1:911313/1‑100 (MQ=255)                                                                                                 |                                                                                              CGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCACTC  >  NC_000913/3573966‑3574155

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