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,7540CT95.3% 124.5 / ‑1.5 43R180H (CGT→CAT) allCallantoate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (22/19);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

GTTTGATTCGCCGTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGACGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGCCC  >  NC_000913/544618‑544885
                                                                                                                                        |                                                                                                                                   
gTTTGATTCGCCGTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGAtgc                                                                                                                                   >  2:293736/1‑139 (MQ=255)
  ttGATTCGCCGTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCgg                                                                                                                                 <  1:76993/139‑1 (MQ=255)
               cAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCg                                                                                                                    <  1:293736/139‑1 (MQ=255)
                       ccGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAag                                                                                                            >  1:886423/1‑139 (MQ=255)
                                            cAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcat                                                                                       <  1:573292/139‑1 (MQ=255)
                                                      caccacGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACa                                                                             >  1:900455/1‑139 (MQ=255)
                                                         cacGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGAtgc                                                                                                                                   <  1:219292/82‑1 (MQ=255)
                                                         cacGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGAtgc                                                                                                                                   >  2:219292/1‑82 (MQ=255)
                                                           cGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTggg                                                                                                                       >  1:840445/1‑92 (MQ=255)
                                                           cGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTggg                                                                                                                       <  2:840445/92‑1 (MQ=255)
                                                                        cccATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGACGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCa                                                           >  1:843769/1‑139 (MQ=255)
                                                                                  tCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCtt                                                                >  2:715809/1‑124 (MQ=255)
                                                                                  tCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCtt                                                                <  1:715809/124‑1 (MQ=255)
                                                                                       cacacaGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCgg                                                                                                               <  2:1395593/72‑1 (MQ=255)
                                                                                       cacacaGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCgg                                                                                                               >  1:1395593/1‑72 (MQ=255)
                                                                                         cacaGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAACCTATTTCCTTTGGCATCACAGATATTCCGCAc                                          >  2:387954/1‑139 (MQ=255)
                                                                                            aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCgg                                                                                                                                 <  2:337714/49‑1 (MQ=255)
                                                                                            aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCgg                                                                                                                                 >  1:337714/1‑49 (MQ=255)
                                                                                            aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTaaa                                                                                                      <  2:1103420/76‑1 (MQ=255)
                                                                                            aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTaaa                                                                                                      >  1:1103420/1‑76 (MQ=255)
                                                                                            aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCAcgtc                                       >  2:981553/1‑139 (MQ=255)
                                                                                                 tGTTCAATATGCAGTTCAACAAAGGCTTTGATATCCTGATGCGGAGTTAGTGGGGCGTTc                                                                                                                 >  1:869787/1‑60 (MQ=255)
                                                                                                 tGTTCAATATGCAGCTCAACAAAGGCTTTGATATCCTGATGCGGAGTTAGTGGGGCGTTc                                                                                                                 <  2:869787/60‑1 (MQ=255)
                                                                                                   ttCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTa                                                                       >  1:553108/1‑100 (MQ=255)
                                                                                                   ttCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTa                                                                       <  2:553108/100‑1 (MQ=255)
                                                                                                        tatGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTCAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTGCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGc                           >  2:949314/1‑139 (MQ=255)
                                                                                                        tatGCAGTTCAACAAAGGCTTTAATATCCTGACGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGc                           >  2:826362/1‑139 (MQ=255)
                                                                                                             aGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGcc                      >  1:1086266/1‑139 (MQ=255)
                                                                                                              gTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGccc                     >  1:472534/1‑139 (MQ=255)
                                                                                                                      aaGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAatat             >  1:1152736/1‑139 (MQ=255)
                                                                                                                         gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCAt                                                          <  1:805580/91‑1 (MQ=255)
                                                                                                                         gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCAt                                                          >  2:805580/1‑91 (MQ=255)
                                                                                                                         gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCa                                           <  1:866618/106‑1 (MQ=255)
                                                                                                                         gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCa                                           >  2:866618/1‑106 (MQ=255)
                                                                                                                         gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATAtttt          >  1:30828/1‑139 (MQ=255)
                                                                                                                         gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATAtttt          >  2:662608/1‑139 (MQ=255)
                                                                                                                            ttAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCTAGCCCAAAAATATTTTTAc       <  1:387954/139‑1 (MQ=255)
                                                                                                                            ttAATATCCTGATGCGGAGTCAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTAc       <  1:949314/139‑1 (MQ=255)
                                                                                                                              aaTATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTg     <  2:30828/139‑1 (MQ=255)
                                                                                                                                 atCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGccc  <  2:1086266/139‑1 (MQ=255)
                                                                                                                                   ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCa                                                                                                <  2:1006841/43‑1 (MQ=255)
                                                                                                                                   ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCa                                                                                                >  1:1006841/1‑43 (MQ=255)
                                                                                                                                   ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTc                                                                                         <  2:566992/50‑1 (MQ=255)
                                                                                                                                   ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTc                                                                                         >  1:566992/1‑50 (MQ=255)
                                                                                                                                    cTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGc                                                                                             <  2:536981/45‑1 (MQ=255)
                                                                                                                                    cTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGc                                                                                             >  1:536981/1‑45 (MQ=255)
                                                                                                                                        |                                                                                                                                   
GTTTGATTCGCCGTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGACGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGCCC  >  NC_000913/544618‑544885

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