Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 444,779 C→T 61.6% W440* (TGG→TAG)  cyoB ← cytochrome o ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273444,7790CT61.6% 20.7 / 39.1 44W440* (TGG→TAG) cyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base C (6/11);  new base T (13/14);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 5.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.24e-01

CTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAG  >  CP009273/444643‑444908
                                                                                                                                        |                                                                                                                                 
ctgGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCa                                                                                                                       >  2:144679/1‑149 (MQ=255)
 tgGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCTCATACCGTGGCATCAAGGCAACGAAGAAGCCGAGGCTCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCAt                                                                                                                      <  2:248380/149‑1 (MQ=255)
 tgGCTCAAACGACGGGTCATGCCCATGAAGCCAAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGAGCCAGAACCAGAACGCGCGTTTGCCCAAGGTTTCGTTCAGTTTGAAAACGAACGCTTTAGGCTACCGCTAGGTCac                                                                                                                      >  2:406476/1‑148 (MQ=255)
        aaCGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCg                                                                                                               >  1:272102/1‑149 (MQ=255)
              gggTCATGCCCATGAAGCCCAGCGCATACTGTGGCCTAAATGCACCGACGAAGCCAATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTGGATACAGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAagcag                                                                                                         <  2:120380/149‑1 (MQ=255)
                gTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGcc                                                                                                       >  1:356304/1‑149 (MQ=255)
                    tGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAGCCGAACGCTTTAGGCTACCGGTAGGTCATCCCTGCGAAGCAGCCGAAg                                                                                                   <  1:249602/149‑1 (MQ=255)
                     gCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGa                                                                                                  >  1:408198/1‑149 (MQ=255)
                       ccATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGAcc                                                                                                >  2:10688/1‑149 (MQ=255)
                            aaGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCAcgcc                                                                                           >  1:130630/1‑149 (MQ=255)
                                    cgcAGACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCgatga                                                                                   <  2:469013/149‑1 (MQ=255)
                                           aGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGAACACgat                                                                            >  2:386661/1‑147 (MQ=255)
                                            gTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTa                                                                           <  2:323651/149‑1 (MQ=255)
                                                 aTAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGGCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGATATGGaa                                                                      >  2:331548/1‑149 (MQ=255)
                                                 aTAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGaa                                                                      <  1:413287/149‑1 (MQ=255)
                                                        cAACGAAGAAGCCGATGATCCAGAACAACGAAAAGCGGTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAGGACCACGCCGCCGATGATCACGGTATGGGCGTGCGCa                                                               >  1:113277/1‑149 (MQ=255)
                                                          aCGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAAt                                                             <  2:130630/149‑1 (MQ=255)
                                                             aagaagCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAg                                                          >  1:244331/1‑149 (MQ=255)
                                                              agaaaCCGATGATCCAGAACCAGAACGCTCGGTTACCCCAGATTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGAAGACCACGGTAGGGCAGGACGCaaccggg                                                         >  2:131770/1‑144 (MQ=255)
                                                                aagCCTATGATCCAGACCCGGAACGCGCGGTTACCCCAGGTTTCGTTCAGTGTGAAACCGAACGCTTTAGGCCACCAGTAGGGCATCCCTACGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGaa                                                       <  2:266697/149‑1 (MQ=255)
                                                                aagCCGATGGTCCAGAAGCCGACCGCGCGCTTGCTCCTGTTTCCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGACGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAggga                                                       <  2:266748/149‑3 (MQ=255)
                                                                aagCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGaa                                                       <  1:70712/149‑1 (MQ=255)
                                                                   ccGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCCGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCGTGCGAAGCAGCCGAAGACCACGCCGCCGATGAACAGGTAATGGAAGTGCGCaaccgcgcacga                                                    >  2:225630/1‑139 (MQ=255)
                                                                     gatgatCCAGAACCCCAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCCGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCagatcggaagagca                                                  <  1:273963/149‑14 (MQ=255)
                                                                      atgatCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCt                                                 <  2:301119/149‑1 (MQ=255)
                                                                         attcagaaccagaacGCGCGTTTACCCCCGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGtt                                              <  2:164130/146‑1 (MQ=255)
                                                                          tCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTa                                             >  2:60084/1‑149 (MQ=255)
                                                                           ccagaaccagaaCGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTAt                                            <  2:356304/149‑1 (MQ=255)
                                                                            ccgcaccagaccGCGCGTTTACCCCAGGCTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATg                                           <  1:6754/145‑1 (MQ=255)
                                                                            cagaaccagaacGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATg                                           <  1:324060/149‑1 (MQ=255)
                                                                                 ccagaacGCGCGTTTACCCCAGGTTTCGTTTAGTTTCAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATTACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGaa                                      <  1:1881/149‑1 (MQ=255)
                                                                                 ccagaacGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGaa                                      <  1:10688/149‑1 (MQ=255)
                                                                                      acGCGCGTTTACCCCAGGTTTCGTTCAGTTCGAAACCAAACGCCTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGACCACGTTATGGAAGTGCGCAATCAGGAACAGGCCGTTATGCAGAACGAAg                                 >  2:430242/1‑149 (MQ=255)
                                                                                      acGCGCGTTTACCCCAGGTTTCGGTCAGTTCGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTAATGCAGAACGAAg                                 >  2:430246/1‑149 (MQ=255)
                                                                                            gTTTTCCCCAGGTTTCGTTCAGTTTGAAACCGACCGCTGTAGGCTACCCGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCcgcg                           <  1:47844/149‑1 (MQ=255)
                                                                                              ttACCCCAGGTTTCGTTCACTTTGAACCCGAACGCTTTATGTTACCAGTAGGTCAGCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATCGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGcc                         <  2:59494/149‑1 (MQ=255)
                                                                                               tatccctgtcttCCTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGccc                        <  2:47975/139‑1 (MQ=255)
                                                                                                aCCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCg                       <  2:408198/149‑1 (MQ=255)
                                                                                                          tCGTTCAGTTTGAAACCGAACGCTTTGGGCTACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGACCACGCTATGGAAGTGCGGAATCAGGAACAGGCGGTTATGCAGAACGACGTCCGCGGCCggacccgcccg             >  2:320721/1‑140 (MQ=255)
                                                                                                           cGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagc            <  1:60084/149‑1 (MQ=255)
                                                                                                                 gTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACg      >  1:297041/1‑149 (MQ=255)
                                                                                                                    tGAAACCGAACGTCTGAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAACACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCa   <  1:19106/149‑1 (MQ=255)
                                                                                                                    ggAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCCGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCa   <  2:294658/148‑1 (MQ=255)
                                                                                                                     gAAACCGAACGCTTTAGGCTACCAGTAGGTGATCCCTGCGAAGCAGCCGAAGACCACGCAGCCGATGATCGAGTTATGGAAGTGCGCAATCAGGAGCGAGCTGTTATGAAGAGAGAAGTCCGCGCCCGGTACGGCCAGCAGCCCGCCac  >  2:75155/1‑148 (MQ=255)
                                                                                                                                        |                                                                                                                                 
CTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAG  >  CP009273/444643‑444908

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑