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_000913544,7540CT97.6% 130.6 / NA 42R180H (CGT→CAT) allCallantoate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (24/17);  minor base A (1/0);  total (25/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TTGATTCGCCGTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGACGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGCCCCA  >  NC_000913/544620‑544887
                                                                                                                                      |                                                                                                                                     
ttGATTCGCCGTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCgg                                                                                                                                   <  1:917933/139‑1 (MQ=255)
ttGATTCGCCGCTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCgg                                                                                                                                   <  1:851684/139‑1 (MQ=255)
        ccGTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTg                                                                                                                           <  1:498475/139‑1 (MQ=255)
                  ttACCGTATAACGACGCTGCCCGACAATTGCAATCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCgg                                                                                                                                   <  2:79791/121‑1 (MQ=255)
                  ttACCGTATAACGACGCTGCCCGACAATTGCAATCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCgg                                                                                                                                   >  1:79791/1‑121 (MQ=255)
                                             ttGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcatcgc                                                                                      >  2:406567/1‑139 (MQ=255)
                                                       cacGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGtt                                                                                                                    <  1:1111049/99‑1 (MQ=255)
                                                       cacGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGtt                                                                                                                    >  2:1111049/1‑99 (MQ=255)
                                                       cacGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACaaaa                                                                            >  2:1102350/1‑139 (MQ=255)
                                                          gCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGAc                                                                                >  1:1080020/1‑132 (MQ=255)
                                                          gCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGAc                                                                                <  2:1080020/132‑1 (MQ=255)
                                                              attgattgCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcatcgc                                                                                      >  1:571120/1‑122 (MQ=255)
                                                              attgattgCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcatcgc                                                                                      <  2:571120/122‑1 (MQ=255)
                                                                                tCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTggg                                                                                                                         <  2:1262142/69‑1 (MQ=255)
                                                                                tCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTggg                                                                                                                         >  1:1262142/1‑69 (MQ=255)
                                                                                          aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCa                                                                                                  >  2:162653/1‑82 (MQ=255)
                                                                                          aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCa                                                                                                  <  1:162653/82‑1 (MQ=255)
                                                                                            cccTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCAcgtcgt                                       >  2:416515/1‑139 (MQ=255)
                                                                                                 ttCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCa                                             <  1:198854/128‑1 (MQ=255)
                                                                                                 ttCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCa                                             >  2:198854/1‑128 (MQ=255)
                                                                                                  tCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGAt                                 >  2:275269/1‑139 (MQ=255)
                                                                                                    aaTATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTc                               >  1:814607/1‑139 (MQ=255)
                                                                                                      tatGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGc                             >  1:996330/1‑139 (MQ=255)
                                                                                                           aggtCAACAAAGGCTTTAATAGCCAGAAGCGGAGCTAGTGTGGCGGTCGGGAGAGTAGATCCGCTAGCCTTCATGGCATCGAGAAATGAATTTCCTGTGGCATCACAGCTATTCCGCATGTCGTCAGAATTCGCTAGcc                        >  2:692622/4‑139 (MQ=255)
                                                                                                           aGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCa                                             >  1:501028/1‑118 (MQ=255)
                                                                                                           aGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCa                                             <  2:501028/118‑1 (MQ=255)
                                                                                                           aGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGcc                        >  2:48096/1‑139 (MQ=255)
                                                                                                           aGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGcc                        >  2:97058/1‑139 (MQ=255)
                                                                                                            gTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTc                                                                                                                   >  2:449856/1‑47 (MQ=255)
                                                                                                            gTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTc                                                                                                                   <  1:449856/47‑1 (MQ=255)
                                                                                                            gTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGccc                       >  2:77823/1‑139 (MQ=255)
                                                                                                             ttCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCa                      >  1:902514/1‑139 (MQ=255)
                                                                                                                aacaaAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCaaaa                   >  2:1091705/1‑139 (MQ=255)
                                                                                                                  agacagCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCCGAAGAGTGAATCCGCAACCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTGAGGATTCGCCAGCCCaaaact                 <  2:789478/134‑3 (MQ=255)
                                                                                                                       gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCaa                                                                                                 >  2:337679/1‑54 (MQ=255)
                                                                                                                       gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCaa                                                                                                 <  1:337679/54‑1 (MQ=255)
                                                                                                                       gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGc                                                              >  2:273829/1‑89 (MQ=255)
                                                                                                                       gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGc                                                              <  1:273829/89‑1 (MQ=255)
                                                                                                                          ttAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTAc         <  2:902514/139‑1 (MQ=255)
                                                                                                                            aaTATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTg       <  2:149077/139‑1 (MQ=255)
                                                                                                                                 ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGCCCCa  >  2:1038027/1‑139 (MQ=255)
                                                                                                                                 ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGCCCCa  >  2:777128/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                     
TTGATTCGCCGTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGACGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGCCCCA  >  NC_000913/544620‑544887

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