Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 444,519 G→A R527C (CGT→TGT)  cyoB ← cytochrome o ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273444,5190GA95.9% 168.9 / ‑1.8 50R527C (CGT→TGT) cyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (26/21);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-01

GGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCGC  >  NZ_CP009273/444439‑444599
                                                                                |                                                                                
ggCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTggtc                                                                         <  2:617131/90‑1 (MQ=255)
   aaaGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTggtcgcg                                                                      >  2:885210/1‑90 (MQ=255)
      gTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTggtcgcggtc                                                                   >  1:199468/1‑90 (MQ=255)
                ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                         >  2:357615/1‑90 (MQ=255)
                ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                         >  2:679397/1‑90 (MQ=255)
                ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                         >  1:456454/1‑90 (MQ=255)
                ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                         >  2:618743/1‑90 (MQ=255)
                ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                         >  1:307317/1‑90 (MQ=255)
                ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                         >  2:521752/1‑90 (MQ=255)
                ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                         >  1:191801/1‑90 (MQ=255)
                ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGGCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                         >  2:823471/1‑90 (MQ=255)
                    gCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGt                                                     <  1:149035/90‑1 (MQ=255)
                      ggggAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTAc                                                   <  2:668524/90‑1 (MQ=255)
                         gAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATc                                                <  2:274119/90‑1 (MQ=255)
                              ggTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAAt                                           <  1:871390/90‑1 (MQ=255)
                                 tGCCCACTTCAGCGTACGTCCACCCCACGGGTCGCCAGGCAGGTCCCCGTTCTGGTCGCGGTCGCGAAAAGAAAAGCACCTGTGAATAAc                                        >  1:737271/1‑90 (MQ=255)
                                 tGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAAc                                        >  1:126915/1‑90 (MQ=255)
                                 tGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAAc                                        >  1:836787/1‑90 (MQ=255)
                                 tGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAAc                                        >  1:345501/1‑90 (MQ=255)
                                       cTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTggtcgc                                                                       <  1:716576/53‑1 (MQ=255)
                                       cTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTggtcgc                                                                       >  2:716576/1‑53 (MQ=255)
                                       cTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa                                  <  1:618743/90‑1 (MQ=255)
                                       cTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa                                  <  2:359436/90‑1 (MQ=255)
                                       cTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa                                  <  1:129226/90‑1 (MQ=255)
                                       cTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa                                  <  2:419142/90‑1 (MQ=255)
                                        tCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCag                                 >  1:93519/1‑90 (MQ=255)
                                            gCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGaa                             >  2:194571/1‑90 (MQ=255)
                                                 cGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATAccc                        >  1:792324/1‑90 (MQ=255)
                                                 cGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATAccc                        >  1:999865/1‑90 (MQ=255)
                                                 cGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa                                  <  1:861449/80‑1 (MQ=255)
                                                 cGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa                                  >  2:861449/1‑80 (MQ=255)
                                                 cGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATAccc                        >  1:362268/1‑90 (MQ=255)
                                                 cGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATAccc                        >  2:299236/1‑90 (MQ=255)
                                                 cGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATAccc                        >  1:490400/1‑90 (MQ=255)
                                                  ggCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCa                       <  2:126915/90‑1 (MQ=255)
                                                   gCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAg                      <  2:199468/90‑1 (MQ=255)
                                                   gCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAg                      >  2:990772/1‑90 (MQ=255)
                                                     cACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAgcg                    <  2:456454/90‑1 (MQ=255)
                                                     cACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAgcg                    <  2:191801/90‑1 (MQ=255)
                                                      aCCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAgcgc                   >  2:127279/1‑90 (MQ=255)
                                                          cACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATc               <  2:196597/90‑1 (MQ=255)
                                                          cACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATc               <  1:127279/90‑1 (MQ=255)
                                                                 cGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa                                  <  1:9107/64‑1 (MQ=255)
                                                                 cGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa                                  >  2:9107/1‑64 (MQ=255)
                                                                 cGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTg        >  2:504399/1‑90 (MQ=255)
                                                                 cGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTg        >  1:997453/1‑90 (MQ=255)
                                                                 cGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTg        >  2:773404/1‑90 (MQ=255)
                                                                     aGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCAcc    <  1:194571/90‑1 (MQ=255)
                                                                       tCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCgc  <  1:101459/90‑1 (MQ=255)
                                                                       tCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCgc  <  1:885210/90‑1 (MQ=255)
                                                                                |                                                                                
GGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCGC  >  NZ_CP009273/444439‑444599

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