Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 539,434 G→A L96L (CTG→CTA allB → allantoinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913539,4340GA100.0% 72.8 / NA 21L96L (CTG→CTAallBallantoinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/11);  total (10/11)

GAAGGTTATGAAACCGGTACTCGCGCAGCGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGG  >  NC_000913/539357‑539507
                                                                             |                                                                         
gAAGGTTATGAAACCGGTACTCGCGCAGCGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGctacct                                                                        <  1:108073/81‑1 (MQ=255)
    gTTATGAAACCGGTACTCGCGCAGCGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTACCTGCaa                                                                    <  2:108074/81‑1 (MQ=255)
       aTGAAACCGGTACTCGCGCAGCGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTACCTGCAACgg                                                                 >  2:108107/1‑81 (MQ=255)
             ccGGTACTCGCGCAGCGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTACCTGCAACGGTTGAcc                                                           <  2:108075/81‑1 (MQ=255)
                 tACTCGCGCAGCGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTACCTGCAACGGTTGACCgcgc                                                       <  1:108076/81‑1 (MQ=255)
                 tACTCGCGCAGCGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTACCTGCAACGGTTGACCgcgc                                                       >  1:108098/1‑81 (MQ=255)
                          aGCGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTACCTGCAACGGTTGACCGCGCTTCAATTg                                               >  2:108111/1‑80 (MQ=255)
                          aGCGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTACCTGCAACGGTTGACCGCGCTTCAATTg                                               >  1:108103/1‑80 (MQ=255)
                            cGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTACCTGCAACGGTTGACCGCGCTTCAATTGAAc                                            >  2:108109/1‑81 (MQ=255)
                                  aaGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTACCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGt                                      <  2:108077/81‑1 (MQ=255)
                                   aGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTACCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGtt                                     <  1:108078/81‑1 (MQ=255)
                                   aGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTACCTGCAACGGTTGACCGCGCTTCAATTGAACTAAAGtt                                     <  2:108079/81‑1 (MQ=255)
                                       ggtATCACCACCATGATCGAAATGCCGCTCAACCAGCTACCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGAt                                 >  1:108110/1‑81 (MQ=255)
                                         tATCACCACCATGATCGAAATGCCGCTCAACCAGCTACCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGc                               <  1:108080/81‑1 (MQ=255)
                                             accaccATGATCGAAATGCCGCTCAACCAGCTACCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGc                            >  2:108105/1‑80 (MQ=255)
                                                  caTGATCGAAATGCCGCTCAACCAGCTACCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAgg                      >  2:108106/1‑81 (MQ=255)
                                                        cGAAATGCCGCTCAACCAGCTACCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCt                >  2:108124/1‑81 (MQ=255)
                                                         gAAATGCCGCTCAACCAGCTACCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTg               <  2:108081/81‑1 (MQ=255)
                                                            aTGCCGCTCAACCAGCTACCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACt            <  2:108082/81‑1 (MQ=255)
                                                             tGCCGCTCAACCAGCTACCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTa           <  1:108083/81‑1 (MQ=255)
                                                                  cTCAACCAGCTACCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGAt      <  2:108084/81‑1 (MQ=255)
                                                                      aCCAGCTACCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCgg  >  2:108115/1‑81 (MQ=255)
                                                                             |                                                                         
GAAGGTTATGAAACCGGTACTCGCGCAGCGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGG  >  NC_000913/539357‑539507

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