Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,156,864	C→T	G5G (GGC→GGT)	argH →	argininosuccinate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,156,864	0	C	T	100.0%	141.6 / NA	44	G5G (GGC→GGT)	argH	argininosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/24);  total (20/24)

TGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTT  >  NC_000913/4156729‑4156998                                                                                                                                        |                                                                                                                                       tGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGgggtgggc                                                                                                                                    >  2:593181/1‑140 (MQ=255)   cGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGt                                                                                                                                  >  2:607373/1‑140 (MQ=255)           cTTCCGGCACTGTTTAAGGGTAGGCCGATGGATACCCGGATTTTAGCTAAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCa                                                                                                                          <  2:378537/140‑1 (MQ=255)           cTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCa                                                                                                                          <  2:994112/140‑1 (MQ=255)             tCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAgg                                                                                                                        >  2:249587/1‑140 (MQ=255)              ccGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGgc                                                                                                                       <  1:962052/140‑1 (MQ=255)                                 tGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCaaa                                                                                                    >  1:585540/1‑140 (MQ=255)                                 tGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCaaa                                                                                                    >  1:576883/1‑140 (MQ=255)                                 tGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATCGCACTTTGGGGTGGGCGTTTTACCCCGGCATCAGATCAACGGTTCaaa                                                                                                    >  2:342942/1‑140 (MQ=255)                                      aTGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAAtt                                                                                               <  2:533851/140‑1 (MQ=255)                                      aTGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAAtt                                                                                               <  2:702194/140‑1 (MQ=255)                                       tGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTc                                                                                              >  2:196562/1‑140 (MQ=255)                                           tACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACg                                                                                          <  2:143229/140‑1 (MQ=255)                                                  aTTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACt                                                                                   <  2:237590/140‑1 (MQ=255)                                                    tttAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTgc                                                                                 <  2:62032/140‑1 (MQ=255)                                                       aGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCt                                                                              <  1:842037/140‑1 (MQ=255)                                                            aaGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGAt                                                                         >  2:20918/1‑140 (MQ=255)                                                                          ggAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGa                                                           <  1:601810/140‑1 (MQ=255)                                                                           gAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAg                                                          <  1:877972/140‑1 (MQ=255)                                                                               cgcAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAgg                                                      <  1:196562/140‑1 (MQ=255)                                                                               cgcAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAgg                                                      <  1:607373/140‑1 (MQ=255)                                                                                 cAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAgg                                                      <  2:198097/138‑1 (MQ=255)                                                                                 cAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAgg                                                      >  1:198097/1‑138 (MQ=255)                                                                                 cAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGat                                                    <  1:670461/140‑1 (MQ=255)                                                                                 cAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGat                                                    >  2:239822/1‑140 (MQ=255)                                                                                    cTTTCCGGCATTGAATTTCAAATTAAGGAAACAGAGTTACGGCACTTTTGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATAtt                                                 <  2:308292/140‑1 (MQ=255)                                                                                       tCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATAttgtt                                              >  1:751797/1‑140 (MQ=255)                                                                                        ccGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATAttgttg                                             <  2:751797/140‑1 (MQ=255)                                                                                         cGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGGGCTTTGATTACTGTTTGGCGGGGCAGGATATTGTTgg                                            >  1:291036/1‑140 (MQ=255)                                                                                         cGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTgg                                            <  1:988912/140‑1 (MQ=255)                                                                                         cGGCATTGAATTTCAAAATAAGCAAACAGAGGTATCGCATTTTGGGGTCGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTgg                                            <  2:414292/140‑1 (MQ=255)                                                                                                   tttCAAAATAAGGAAACAGAGTTATGGCGCTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGcc                                  <  2:950206/140‑1 (MQ=255)                                                                                                    ttCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCt                                 <  1:249587/140‑1 (MQ=255)                                                                                                        aaaTAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTc                             >  2:229668/1‑140 (MQ=255)                                                                                                          aTAAGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCa                           <  1:593181/140‑1 (MQ=255)                                                                                                             aGGAAACAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAg                        >  2:534334/1‑140 (MQ=255)                                                                                                                   cAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTTTGGCGGAGCAGGATATTGTTGGCTCCGTGGCCTGGTCCAAAGCCCTgg                  >  2:346499/1‑140 (MQ=255)                                                                                                                   cAGAGTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTATGGCGGAGTAGGGTATTGTTGGCTCGGTGGGCTGGTCCAAAGCCCTgt                  >  2:275125/1‑139 (MQ=255)                                                                                                                      agTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCa               >  2:551484/1‑140 (MQ=255)                                                                                                                      agTTATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCa               >  2:627030/1‑140 (MQ=255)                                                                                                                        ttATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGCGGTCTGGTCTAAAGCCCTGGTCACg             >  1:426671/1‑140 (MQ=255)                                                                                                                         tATGGCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTCGCCTGGTCCAAAGCCCTGGTCACgg            <  2:691745/140‑1 (MQ=255)                                                                                                                             gCACTTTGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAgg        >  1:172901/1‑140 (MQ=255)                                                                                                                                  ttGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCCGTGGCCTAGTCGAAAGCCCTGATCACGGTAGGCgtgt   >  1:309831/1‑140 (MQ=255)                                                                                                                                   tGGGGTGGGCGTTTTTCCCAGGCAGTAGATCAACGGTGCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGtt  <  1:484007/140‑1 (MQ=255)                                                                                                                                   tGGGGTGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGtt  <  1:229668/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                       TGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTT  >  NC_000913/4156729‑4156998

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