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%	149.3 / NA	46	R180H (CGT→CAT)	allC	allantoate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/24);  total (22/24)

TTGATTCGCCGTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGACGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGCCCCA  >  NC_000913/544620‑544887                                                                                                                                       |                                                                                                                                      ttGATTCGCCGTTCAGCGTTATCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTACCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCgg                                                                                                                                   <  1:796417/139‑1 (MQ=255) ttGATTCGCCGTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCgg                                                                                                                                   <  2:469641/139‑1 (MQ=255)         ccGTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTg                                                                                                                           <  2:1157066/139‑1 (MQ=255)                       cGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAaga                                                                                                             >  1:333427/1‑139 (MQ=255)                                acgCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCg                                                                                                    >  1:1030981/1‑139 (MQ=255)                                acgCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCg                                                                                                    <  1:1144312/139‑1 (MQ=255)                                 cgCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCACATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGc                                                                                                   <  2:1191394/139‑1 (MQ=255)                                   cTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTggg                                                                                                                         <  2:423801/115‑1 (MQ=255)                                   cTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTggg                                                                                                                         >  1:423801/1‑115 (MQ=255)                                      cccGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGcc                                                                                              >  1:1118748/1‑139 (MQ=255)                                              ttGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcatcgc                                                                                      >  2:172994/1‑139 (MQ=255)                                                gCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcatcgcat                                                                                    <  2:534953/139‑1 (MQ=255)                                                    tCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGAc                                                                                <  2:1030981/139‑1 (MQ=255)                                                           gCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTa                                                                         <  2:333427/139‑1 (MQ=255)                                                                ttgattgCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACa                                                                               <  2:624203/128‑1 (MQ=255)                                                                ttgattgCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACa                                                                               >  1:624203/1‑128 (MQ=255)                                                                ttgattgCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTa                                                                         >  1:323586/1‑134 (MQ=255)                                                                ttgattgCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTa                                                                         <  2:323586/134‑1 (MQ=255)                                                                ttgattgCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTcc                                                                    >  1:1161818/1‑139 (MQ=255)                                                                          aTTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATca                                                          <  1:1357251/139‑1 (MQ=255)                                                                                      cacacaGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGc                                              >  2:1244258/1‑139 (MQ=255)                                                                                       acacaGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCa                                             <  2:1161818/139‑1 (MQ=255)                                                                                         acaGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcatcgcat                                                                                    >  1:1224478/1‑98 (MQ=255)                                                                                         acaGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcatcgcat                                                                                    <  2:1224478/98‑1 (MQ=255)                                                                                           aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcatcgc                                                                                      <  2:1326786/94‑1 (MQ=255)                                                                                           aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcatcgc                                                                                      >  1:1326786/1‑94 (MQ=255)                                                                                                       tatGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTgg                                                               >  1:385308/1‑105 (MQ=255)                                                                                                       tatGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTgg                                                               <  2:385308/105‑1 (MQ=255)                                                                                                       tatGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGc                             >  2:862024/1‑139 (MQ=255)                                                                                                       tatGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGc                             >  2:658980/1‑139 (MQ=255)                                                                                                       tatGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGc                             >  1:342285/1‑139 (MQ=255)                                                                                                          gCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAg                          <  2:1118748/139‑1 (MQ=255)                                                                                                            aGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGa                                                                                                                <  2:787758/51‑1 (MQ=255)                                                                                                            aGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGa                                                                                                                >  1:787758/1‑51 (MQ=255)                                                                                                            aGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCAcgtcgt                                       >  2:1326298/1‑124 (MQ=255)                                                                                                            aGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCAcgtcgt                                       <  1:1326298/124‑1 (MQ=255)                                                                                                                caacaaAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCaaa                    <  2:342285/139‑1 (MQ=255)                                                                                                                 aacaaAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGtt                                                                                                                    <  1:539688/42‑1 (MQ=255)                                                                                                                 aacaaAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGtt                                                                                                                    >  2:539688/1‑42 (MQ=255)                                                                                                                 aacaaAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGtt                                                                                                                    >  1:1080588/1‑42 (MQ=255)                                                                                                                 aacaaAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGtt                                                                                                                    <  2:1080588/42‑1 (MQ=255)                                                                                                                      aGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATAtt              >  1:326405/1‑139 (MQ=255)                                                                                                                                 tCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGcccc   <  1:862024/139‑1 (MQ=255)                                                                                                                                  ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCaa                                                                                                 >  2:1027742/1‑44 (MQ=255)                                                                                                                                  ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCaa                                                                                                 <  1:1027742/44‑1 (MQ=255)                                                                                                                                  ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGCCCCa  >  1:50552/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                      TTGATTCGCCGTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGACGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGCCCCA  >  NC_000913/544620‑544887

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