Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	544,754	C→T	R180H (CGT→CAT)	allC ←	allantoate amidohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	544,754	0	C	T	100.0%	185.0 / NA	58	R180H (CGT→CAT)	allC	allantoate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/27);  total (31/27)

CCGTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGACGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGCCCCAG  >  NC_000913/544628‑544888                                                                                                                               |                                                                                                                                       ccGTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTg                                                                                                                            <  1:997643/139‑1 (MQ=255)        gCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGtt                                                                                                                     <  1:38745/139‑1 (MQ=255)                         cgCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGc                                                                                                    >  2:1073666/1‑139 (MQ=255)                           cTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCaa                                                                                                  >  2:1061426/1‑139 (MQ=255)                             gCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGc                                                                                                <  1:1073666/139‑1 (MQ=255)                                      ttGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcatcgc                                                                                       <  2:1019160/139‑1 (MQ=255)                                              accacGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCgg                                                                                                                                    <  2:466422/86‑1 (MQ=255)                                              accacGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCgg                                                                                                                                    >  1:466422/1‑86 (MQ=255)                                              accacGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTg                                                                                                                            >  2:397774/1‑94 (MQ=255)                                              accacGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTg                                                                                                                            <  1:397774/94‑1 (MQ=255)                                              accacGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcat                                                                                          >  2:109360/1‑128 (MQ=255)                                              accacGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcat                                                                                          <  1:109360/128‑1 (MQ=255)                                                             tgCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTggg                                                                                                                          >  2:324054/1‑81 (MQ=255)                                                             tgCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTggg                                                                                                                          <  1:324054/81‑1 (MQ=255)                                                             tgCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACaa                                                                               >  1:607143/1‑124 (MQ=255)                                                             tgCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACaa                                                                               <  2:607143/124‑1 (MQ=255)                                                             tgCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTgg                                                                >  1:193944/1‑139 (MQ=255)                                                                ccATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCAt                                                             <  1:1061426/139‑1 (MQ=255)                                                                         tCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTgg                                                                >  2:1245422/1‑127 (MQ=255)                                                                         tCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTgg                                                                <  1:1245422/127‑1 (MQ=255)                                                                             gCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCgg                                                                                                                                    <  1:441886/55‑1 (MQ=255)                                                                             gCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCgg                                                                                                                                    >  2:441886/1‑55 (MQ=255)                                                                              cacacaGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTggg                                                                                                                          >  1:434479/1‑64 (MQ=255)                                                                              cacacaGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTggg                                                                                                                          <  2:434479/64‑1 (MQ=255)                                                                                cacaGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCAc                                             >  1:1334188/1‑139 (MQ=255)                                                                                   aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAg                                                        >  2:1327270/1‑125 (MQ=255)                                                                                   aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAg                                                        <  1:1327270/125‑1 (MQ=255)                                                                                   aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCAcgtc                                          >  1:453942/1‑139 (MQ=255)                                                                                   aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGATAGTGGGGCGTTCGGAAGAGTGAAGCCGCAAGCCTTAATCGCATCGACAAAACTATTACCTTTGGCATAACAGATATTCCGCAcgtc                                          >  2:305875/1‑139 (MQ=255)                                                                                     cccTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTggg                                                                                                                          >  2:1258932/1‑57 (MQ=255)                                                                                     cccTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTggg                                                                                                                          <  1:1258932/57‑1 (MQ=255)                                                                                          ttCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCg                                                                                                     >  2:517945/1‑73 (MQ=255)                                                                                          ttCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCg                                                                                                     <  1:517945/73‑1 (MQ=255)                                                                                               tatGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGc                              >  2:820429/1‑139 (MQ=255)                                                                                               tatGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGc                              >  1:365758/1‑139 (MQ=255)                                                                                                    aGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcatcgc                                                                                       >  1:84039/1‑77 (MQ=255)                                                                                                    aGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcatcgc                                                                                       <  2:84039/77‑1 (MQ=255)                                                                                                      ttCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACa                                                                                <  2:418989/82‑1 (MQ=255)                                                                                                      ttCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACa                                                                                >  1:418989/1‑82 (MQ=255)                                                                                                       tCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCaa                      >  2:387993/1‑139 (MQ=255)                                                                                                           caaAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCtt                                                                   <  1:363432/90‑1 (MQ=255)                                                                                                           caaAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCtt                                                                   >  2:363432/1‑90 (MQ=255)                                                                                                              aGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGATTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGc                                                               >  2:950684/1‑91 (MQ=255)                                                                                                              aGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGc                                                               <  1:950684/91‑1 (MQ=255)                                                                                                                gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTc                                                                                                                    >  2:1178155/1‑36 (MQ=255)                                                                                                                gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTc                                                                                                                    <  1:1178155/36‑1 (MQ=255)                                                                                                                gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTaa                                                                                                          >  1:790250/1‑46 (MQ=255)                                                                                                                gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCtt                                                                   <  2:201706/85‑1 (MQ=255)                                                                                                                gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCtt                                                                   >  1:201706/1‑85 (MQ=255)                                                                                                                gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCaaa                     >  1:251297/1‑131 (MQ=255)                                                                                                                gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCaaa                     <  2:251297/131‑1 (MQ=255)                                                                                                                gCTTTAATATCCTGATGCGGAGTTAATGGGGCGTTCGGAAGAGTaa                                                                                                          <  2:790250/46‑1 (MQ=255)                                                                                                                          ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcatcgca                                                                                      <  2:1004195/56‑1 (MQ=255)                                                                                                                          ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcatcgca                                                                                      >  1:1004195/1‑56 (MQ=255)                                                                                                                          ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACa                                                                                >  1:1094552/1‑62 (MQ=255)                                                                                                                          ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACa                                                                                <  2:1094552/62‑1 (MQ=255)                                                                                                                          ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGCCCCa   >  1:367392/1‑139 (MQ=255)                                                                                                                           cTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGCCCCag  <  2:453942/139‑1 (MQ=255)                                                                                                                               |                                                                                                                                       CCGTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGACGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGCCCCAG  >  NC_000913/544628‑544888

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