Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,229,765	A→G	A126A (GCT→GCC)	pepE ←	peptidase E, alpha‑aspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,229,765	0	A	G	100.0%	225.9 / NA	60	A126A (GCT→GCC)	pepE	peptidase E, alpha‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/33);  total (27/33)

ACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCA  >  NC_000913/4229630‑4229895                                                                                                                                        |                                                                                                                                   aCGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGggcaggca                                                                                                                                 >  2:15102/1‑139 (MQ=255)               tATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTgc                                                                                                                   >  1:674567/1‑139 (MQ=255)                 tGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTgcgc                                                                                                                 >  2:387919/1‑139 (MQ=255)                  ggCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCt                                                                                                                <  1:818459/139‑1 (MQ=255)                    ccTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTcc                                                                                                              <  1:15102/139‑1 (MQ=255)                    ccTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTcc                                                                                                              >  2:448582/1‑139 (MQ=255)                      ttCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAg                                                                                                            <  2:460963/139‑1 (MQ=255)                      ttCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAg                                                                                                            <  1:251692/139‑1 (MQ=255)                       tCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGc                                                                                                           >  1:324911/1‑139 (MQ=255)                               cgcgTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATAc                                                                                                   >  2:1036745/1‑139 (MQ=255)                                      gTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAgcgc                                                                                            <  2:1234538/139‑1 (MQ=255)                                       tGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGcc                                                                                           <  2:385939/139‑1 (MQ=255)                                                ggTTGATTTGCAGCGGGAACAGATTTAGTGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAAc                                                                                  <  1:263380/139‑1 (MQ=255)                                                ggTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAgcgc                                                                                            >  2:769435/1‑129 (MQ=255)                                                ggTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAgcgc                                                                                            <  1:769435/129‑1 (MQ=255)                                                ggTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAAc                                                                                  <  2:802087/139‑1 (MQ=255)                                                ggTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAAc                                                                                  <  1:679487/139‑1 (MQ=255)                                                ggTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAAc                                                                                  <  2:324911/139‑1 (MQ=255)                                                               ggAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAAcc                                                                                 >  2:710985/1‑125 (MQ=255)                                                               ggAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAAcc                                                                                 <  1:710985/125‑1 (MQ=255)                                                               ggAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAAc                                                                                  >  1:379867/1‑124 (MQ=255)                                                               ggAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAAc                                                                                  <  2:379867/124‑1 (MQ=255)                                                                gAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAg                                                                                                            <  2:371197/97‑1 (MQ=255)                                                                gAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAg                                                                                                            >  1:371197/1‑97 (MQ=255)                                                                         tAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCagcagc                                                         >  2:1426488/1‑139 (MQ=255)                                                                                tCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTcc                                                                                                              >  2:10013/1‑79 (MQ=255)                                                                                tCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTcc                                                                                                              <  1:10013/79‑1 (MQ=255)                                                                                tCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAgcgc                                                                                            >  1:749246/1‑97 (MQ=255)                                                                                tCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAgcgc                                                                                            <  2:749246/97‑1 (MQ=255)                                                                                      ccTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGc                                            <  1:1036745/139‑1 (MQ=255)                                                                                         tGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTcc                                                                                                              <  1:1491301/70‑1 (MQ=255)                                                                                         tGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTcc                                                                                                              >  2:1491301/1‑70 (MQ=255)                                                                                         tGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACt                                         >  1:784112/1‑139 (MQ=255)                                                                                              aTCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTcc                                                                                                              >  1:154655/1‑65 (MQ=255)                                                                                              aTCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTcc                                                                                                              <  2:154655/65‑1 (MQ=255)                                                                                                cGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAg                                  >  1:958462/1‑139 (MQ=255)                                                                                                    aaTCGGCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAAc                              <  1:542120/139‑1 (MQ=255)                                                                                                         gCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTcc                                                                                                              <  1:1106994/54‑1 (MQ=255)                                                                                                         gCATATCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTcc                                                                                                              >  2:1106994/1‑54 (MQ=255)                                                                                                             atCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACt                             >  1:670643/1‑131 (MQ=255)                                                                                                             atCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACt                             <  2:670643/131‑1 (MQ=255)                                                                                                              tCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAAc                                                                                  <  2:199978/77‑1 (MQ=255)                                                                                                              tCGTTGGTGGTACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAAc                                                                                  >  1:199978/1‑77 (MQ=255)                                                                                                                        tACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTcc                                                                                                              <  2:789293/39‑1 (MQ=255)                                                                                                                        tACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTcc                                                                                                              <  1:683363/39‑1 (MQ=255)                                                                                                                        tACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTcc                                                                                                              >  2:683363/1‑39 (MQ=255)                                                                                                                        tACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTcc                                                                                                              >  1:789293/1‑39 (MQ=255)                                                                                                                        tACGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCccgccgcc          >  1:808855/1‑139 (MQ=255)                                                                                                                         aCGAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCccgccgccg         >  1:663915/1‑139 (MQ=255)                                                                                                                           gAATAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGAc       >  2:190530/1‑139 (MQ=255)                                                                                                                            aaTAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCTAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACa      <  1:967914/139‑1 (MQ=255)                                                                                                                            aaTAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACa      <  2:1436205/139‑1 (MQ=255)                                                                                                                            aaTAGTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACa      <  2:766140/139‑1 (MQ=255)                                                                                                                                gTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCa  <  1:190530/139‑1 (MQ=255)                                                                                                                                gTTGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCa  <  1:1107792/139‑1 (MQ=255)                                                                                                                                 ttGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAgcg                                                                                             <  1:307131/47‑1 (MQ=255)                                                                                                                                 ttGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAgcg                                                                                             >  2:307131/1‑47 (MQ=255)                                                                                                                                  tGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACa                                                                                                  <  2:671288/41‑1 (MQ=255)                                                                                                                                  tGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACa                                                                                                  >  1:671288/1‑41 (MQ=255)                                                                                                                                  tGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCa  <  2:1226767/137‑1 (MQ=255)                                                                                                                                  tGGGCAGGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCa  >  1:1226767/1‑137 (MQ=255)                                                                                                                                        |                                                                                                                                   ACGGGTTTCACCTTTATGGCCTTCCGGCAGCGCGTTGGTGAAGTGCGGGTTGATTTGCAGCGGGAACAGATTTAGCGCATCGAAACCTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCA  >  NC_000913/4229630‑4229895

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