Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 870,468 A→G Q678Q (CAA→CAG pdeI → EAL domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273870,4680AG100.0% 74.6 / NA 21Q678Q (CAA→CAGpdeIEAL domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/12);  total (9/12)

CATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGTTGAGATCACTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAAGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGACAGCATAAAGGTGGACTCGCAATTCACCGGCGATATTTTCGGTACTGCGGGAAA  >  NZ_CP009273/870337‑870596
                                                                                                                                   |                                                                                                                                
cATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGTTGAGATCACTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTa                                                                                                                            <  1:498501/138‑1 (MQ=255)
      cGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGTTGAGATCACTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGatat                                                                                                                     <  2:727531/139‑1 (MQ=255)
                                                 gCCTCTGTGTTGAGATCACTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATtgat                                                                                                             <  1:139231/104‑1 (MQ=255)
                                                 gCCTCTGTGTTGAGATCACTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATtgat                                                                                                             >  2:139231/1‑104 (MQ=255)
                                                          ttGAGATCACTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGaaa                                                                  <  2:345419/138‑1 (MQ=255)
                                                                  aCTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATc                                                            >  1:214607/1‑136 (MQ=255)
                                                                  aCTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATc                                                            >  1:826957/1‑136 (MQ=255)
                                                                  aCTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCa                                                          <  1:853332/138‑1 (MQ=255)
                                                                   cTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGa                                                        <  1:611723/139‑1 (MQ=255)
                                                                   cTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGa                                                        <  2:214607/139‑1 (MQ=255)
                                                                      aaaCCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGCCTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGAc                                                       <  2:369511/137‑1 (MQ=255)
                                                                       aaCCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGACAg                                                     >  1:312409/1‑138 (MQ=255)
                                                                                  aGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGACAGCATAAAGGTGGa                                         >  1:484114/1‑139 (MQ=255)
                                                                                    cGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGACAGCATAAAGGTGGACt                                       <  2:644401/139‑1 (MQ=255)
                                                                                                      tATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGACAGCATAAAGGTGGACTCGCAATTCACCGGCGat                      <  2:76941/138‑1 (MQ=255)
                                                                                                                        cTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGa                                                        <  1:551920/86‑1 (MQ=255)
                                                                                                                        cTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGa                                                        >  2:551920/1‑86 (MQ=255)
                                                                                                                        cTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGACAGCATAAAGGTGGACTCGCAATTCACCGGCGATATTTTCGGTACTGCGGGa    >  1:564091/1‑138 (MQ=255)
                                                                                                                          gCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGACAGCATAAAGGTGGACTCGCAATTCACCGGCGATATTTTCGGTACTGCGGGaaa  >  2:274953/1‑138 (MQ=255)
                                                                                                                           cGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGa                                                        >  1:131543/1‑83 (MQ=255)
                                                                                                                           cGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGa                                                        <  2:131543/83‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                
CATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGTTGAGATCACTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAAGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGACAGCATAAAGGTGGACTCGCAATTCACCGGCGATATTTTCGGTACTGCGGGAAA  >  NZ_CP009273/870337‑870596

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