Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 495,361 A→G 17.4% T23T (ACT→ACC)  aes ← acetyl esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273495,3610AG17.4% 50.5 / 5.5 23T23T (ACT→ACC) aesacetyl esterase
Reads supporting (aligned to +/- strand):  ref base A (9/10);  new base G (2/2);  total (11/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

ATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTATCTGGCTGCGGACAAAAGAGACGTGTTTCCACCTGCCCATATTTTGTTGGAACCATGTAAGCTCTGGTTGCCATTTCTGGAGCGCCCGCATTCCAGAATCGGCGCTCAAGCGTGTAATACTGTCGTTGCTCAGCAATCGTTCCCGTTGCGGGCCAGGGCGGTAAATCCGGCTGAAGAGTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACTTC  >  CP009273/495118‑495619
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                  
aTCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTATCTGGCTGCGGACAAAAGAGACGTGTTTCCACCTGCCCATATTTTGTTGGAACCATGTAAGCTCTGGTTGCCATTTCTGGAGCGCCCGCATTCCAGAATCGGCGCTCAAGCGTGTAATACTGTCGTTGCTCAGCAATCGTTCCCGTTGCGGGCCAGGGCGGTAAATCCGGCTGAAGAGTATTCACAACGGTCTTCATTTcagcag                                                                                                                                                                                                                                        <  1:175510/272‑1 (MQ=255)
                                  catCCATGCAAATAAAATAGCGGTGCTTGGCTATCTGGCTGCGGACAAAAGAGACGTGTTTCCACCTGCCCATATTTTGTTGGAACCATGTAAGCGCTGGTTGCCATTTCTGGAGCGCCCGCATTCCAGAATCGGCGCTCAAGCGTGTAATACTGTCGTTGCTCAGCAATCGTTCCCGTTGCGGGCCAGGGCGGTAAATCCGGCTGAAGAGTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCtt                                                                                                                                                                                                   <  1:263721/273‑1 (MQ=255)
                                             ataaaataGCGTCGCTGGGCTATCTGGCTGCGGACAAAAGAGACGTGTTTCCACCTGCCCATATTTTGTTGGAACCATGTAAGCTCTGGTTGCCATTTCTGGAGCGCCCGCATTCCAGAATCGGCGCTCAAGCGGGTAATACTGTCGTTGCTCAGCAATCGTTCCCGTTGCGGGCCAGGGAGGTAAATCCGGCTGAAGAGAATTAACAACGTACTTCATTTCAGCAGAAATAAGGTCCCGACCAGGTAGTTTGTTTTCTGGATNNTTCCaaacag                                                                                                                                                                                        >  2:60786/1‑272 (MQ=255)
                                                                ttATCTGGCTGCGGACAAAAGAGACGTGTTTCCACCTGCCCTTATTTTGTTGGAACCATGTAAGCTCTGGTTGCCATTTCTGGAGCGCCCGCATTACAGAATCGGCGCTCAAGCGTGTAATACTGTCGTTGCTCAGCAATCGTGCCCGTTGCGGGCCAGGGCGGTAAATCCGGCTGAAGAGTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCAt                                                                                                                                                                     <  1:60786/274‑1 (MQ=255)
                                                                                       aCGTGTTTCCACCTGCCCATATTTTGTTGGAACCATGTAAGCTCTGGTTGCCATTTCTGGAGCGCCCGCATTCCAGAATCGGCGCTCAAGCGTGTAATACTGTCGTTGCTCAGCAATCGTTCCCGTTGCGGGCCAGGGCGGTAAATCCGGCTGAAGGGTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATAAAAAAACACCTTTCAAATTACGTAATTTTAAGGAACCACTGCCATTaaa                                                                                                                                              >  2:196306/1‑275 (MQ=255)
                                                                                           gtTTCCACCTGCCCATATTTTGTTGGAACCATGTAAGCTCTGGTTGCCATTTCTGGAGCGCCCGCATTCCAGAATCGGCGCTCAAGCGTGTAATACTGTCGTTGCTCAGCAATCGTTCCCGTTGCGGGCCAGGGCGGTAAATCCGGCTGAAGGGTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTg                                                                                                                                                      <  1:196308/263‑1 (MQ=255)
                                                                                                                     aaCCATGTAAGCTCTGGTTGCCATTTCTGGAGCGCCCGCATTCCAGAATCGGCGCTCAAGCGTGTAATACTGTCGTTGCTCAGCAATCGTTCCCGTTGCGGGCCAGGGCGGTAAATCCGGCTGAAGAGTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAg                                                                                                                                                               >  2:221864/1‑228 (MQ=255)
                                                                                                                         aTGTAAGCTCTGGTTGCCATTTCTGGAGCGCCCGCATTCCAGAATCGGCGCTCAAGCGTGTAATACTGTCGTTGCTCAGCAATCGTTCCCGTTGCGGGCCAGGGCGGTAAATCCGGCTGAAGAGTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAAc                                                                                                                                                           <  1:221866/228‑1 (MQ=255)
                                                                                                                          tGTAAGCTCTGGTTGCCATTTCTGGAGCGCCCGCATTCCAGAATCGGCGCTCAAGCGTGTAATACTGTCGTTGCACAGCAATCGTTCCCGTTGCGTGCCTGGGCGGTAAATCCGGCAGAAGAGTATTCACAACGGTCTTCATTTCAGCAGAACAAAGGTCCAGAACAGGTAGATTGTTTTCCg                                                                                                                                                                                                       >  2:188621/1‑183 (MQ=255)
                                                                                                                              aGCTCTGGTTGCCATTTCTGGAGCGCCCGCATTCCAGAATCGGCGGTCAAGCGTGTAATACTGTCGTTGCTCAGCAATCGTTCCCGTTGCGGGCCAGGGCGGTAAATCCGGCTGAAGAGTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCtt                                                                                                                                                                                                   <  1:188622/183‑1 (MQ=255)
                                                                                                                                   tGGTTGCCATTTCTGGAGCGCCCGCATTCCAGAATCGGCGCTCAAGCGTGTAATACTGTCGTTGCTCAGCAATCGTTCCCGTTGCGGGCCAGGGCGGTAAATCCGGCTGAaga                                                                                                                                                                                                                                                                    <  2:510931/113‑1 (MQ=255)
                                                                                                                                           aTTTCTGGAGCGCCCGCATTCCAGAATCGGCGCTCAAGCGTGTAATACTGTCGTTGCTCAGCAATCGTTCCCGTTGCGGGCCAGGGCGGTAAATCCGGCTGAAGAGTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAAt                                                                                                                                                                             >  1:202299/1‑192 (MQ=255)
                                                                                                                                               cTGGAGCGCCCGCATTCCAGAATCGGCGCTCAAGCGTGTAATACTGTCGTTGCTCAGCAATCGTTCCCGTTGCGGGCCAGGGCGGTAAATCCGGCTGTAGAGTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGa                                                                                      <  1:303114/275‑1 (MQ=255)
                                                                                                                                                                                     tAATACTGTCGTTGCTCAGCAATCGTTCCCGTTGCGGGCCAGGGCGGTAAATCCGGCTGAAGAGTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTAGCGAGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTaaa                                                >  1:96879/1‑275 (MQ=255)
                                                                                                                                                                                             tCGTTGCTCAGCAATCGTTCCCGTTGCGGGCCAGGGCGGTAAATCCGGCTGAAGAGTATTc                                                                                                                                                                                                                                                              >  1:159691/1‑61 (MQ=255)
                                                                                                                                                                                                 tGCTCAGCAATCGTTCCCGTTGCGGGCCAGGGCGGTAAATCCGGCTGAAGAGTATTCACaa                                                                                                                                                                                                                                                          <  2:159690/61‑1 (MQ=255)
                                                                                                                                                                                                  gCTCAGCAATCGTTCCCGTTGCGGGCCAGGGCGGTAAATCCGGCTGAAGAGTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGtt                                                             >  2:486855/1‑249 (MQ=255)
                                                                                                                                                                                                      aGCAATCGTTCCCGTTGCGGGCCAGGGCGGTAAATCCGGCTGAAGAGTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAacc                                                         <  1:486862/249‑1 (MQ=255)
                                                                                                                                                                                                                         ggCCAGGGCGGTAAATCCGGCTGAAGGGTATTCACAATGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGa                                                                                                                                                             >  2:26580/1‑130 (MQ=255)
                                                                                                                                                                                                                             aGGGCGGTAAATCCGGCTGAAGGGTATTCACAATGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCa                                                                                                                                                         <  1:26580/130‑1 (MQ=255)
                                                                                                                                                                                                                                    tAAATCCGGCTGAAGAGTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTTTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCACGGTAAACGTAAAACCAAAAATTAcaac  >  2:166458/1‑271 (MQ=255)
                                                                                                                                                                                                                                        tCCGGCTGAAGAGTATTCACCACGGTCTTAATTTCAGCAGAAATa                                                                                                                                                                                                                                   <  1:393991/45‑1 (MQ=255)
                                                                                                                                                                                                                                         ccGGCTGAAGAGTAGTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGt                                                                                          >  1:203214/1‑181 (MQ=255)
                                                                                                                                                                                                                                             cggAAGAGTATTCACCACGGTCTTCATTTCCACAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCTCTGCCACTATATTGCGATCCCGCCTGCTGATATTGAAACTCGGTGCGTCTCGCGCCCTCCCGTCAGa                                                                                      <  2:203212/179‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                  
ATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCAAATAAAATAGCGTCGCTGGGCTATCTGGCTGCGGACAAAAGAGACGTGTTTCCACCTGCCCATATTTTGTTGGAACCATGTAAGCTCTGGTTGCCATTTCTGGAGCGCCCGCATTCCAGAATCGGCGCTCAAGCGTGTAATACTGTCGTTGCTCAGCAATCGTTCCCGTTGCGGGCCAGGGCGGTAAATCCGGCTGAAGAGTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACTTC  >  CP009273/495118‑495619

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