Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 2,085,482 C→G 24.4% intergenic (+71/+20) deaD → / ← ygiD ATP‑dependent RNA helicase DeaD/4,5‑DOPA dioxygenase extradiol

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1222,085,4820CG24.4% 34.1 / 3.1 35intergenic (+71/+20)deaD/ygiDATP‑dependent RNA helicase DeaD/4,5‑DOPA dioxygenase extradiol
Reads supporting (aligned to +/- strand):  ref base C (9/17);  new base G (5/4);  total (14/21)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAAT  >  ACB122/2085354‑2085608
                                                                                                                                |                                                                                                                                
cAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTAGCGGGCGTGccc                                                                                                                            >  5:105577/1‑135 (MQ=255)
    gACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCGGGCGTAACGAAAAAATCCCCGGCTTAGGGCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTGGCCCGCGa                                                                                                                       <  4:382735/135‑1 (MQ=255)
    gACACTCGCGGGGCGGGCAAGCCGCTGAAGCTGCGTCGCTACCGCGCGTATTGATCCTGGCGAAACCAAAAAATCCCCGGCTTAGGGCGGGTATTTTTTTTACCTGACCCGCCCCTTCCGCGGGCTTCCCCGCGa                                                                                                                       <  6:380997/135‑1 (MQ=255)
     acacTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCcgcgc                                                                                                                       >  8:24029/1‑133 (MQ=255)
     acacTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGaa                                                                                                                      >  6:123110/1‑135 (MQ=255)
           gcgTGGCTGGCAAGCAGCTGAAGATGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTAACGGCCATTTAGCCGGCTTTCCCGCGAAGCGGcc                                                                                                                >  6:60893/1‑135 (MQ=255)
            cgTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGGGTGCCCGCGAAGAGGccc                                                                                                               >  7:163922/1‑135 (MQ=255)
            cgTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGGCCGCGCAGAGGgcc                                                                                                               >  1:308195/1‑135 (MQ=255)
                  gggcAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACCAAAAAATCCCCGGCTTAGGGCGGGTATTTTTTTTACCTGTACCGGCCCTTGCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGt                                                                                                        <  1:296081/134‑1 (MQ=255)
                    gcAAGCAGCTGAAGCTGGGTCGCTACCGCGAGTATTGAGCGTGGCTTAACTAAAAAATCCCCGGGTTAGGTCGGGTTTTTTTTTTACCTGTACCGGCCCTTTCGCGGGGTTGCCCGGGAAGAGGCCCTTCAGTCg                                                                                                      <  8:62443/135‑1 (MQ=255)
                    gcAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTAGCGGGGGTGGCCGCGCAGGGGCCCGTGAGtg                                                                                                       >  6:360386/1‑133 (MQ=255)
                          agcTGAAGCTGCGGCGCTACCGCGAGTATTTATCCGGGCGTAACTAAAAAATCCCGGGCTTCGGGCGGGTATTTTTTTTACCTGTACCGCCCCTTTCGCGGGCTGGCCCGCGAAGAGGCCCTTCAGTCGAACCGa                                                                                                <  5:123110/135‑1 (MQ=255)
                                aGCTGCGTGGCTTCCGCGAGTATTGAGCCTGGCTTAACGAAAAAATCCCCTGCTTAGGGCGGGTTTTTTTATTACCTGACCCGGCCCTTGCCGGGGGGCCCCGGGGAAGAGGCCCTTCAGTCGAACCGATAGAt                                                                                           <  6:214984/134‑1 (MQ=255)
                                aGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTACCGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTCCCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGAt                                                                                           <  7:267291/134‑1 (MQ=255)
                                         gCTACCGCGACGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAACAGGCCCATCAGTCGAACCGAAAGATGTCGATcccc                                                                                 >  8:346742/1‑135 (MQ=39)
                                              cgcgAGGATTGACCCGGGCGTACCGAAAAAATCCCCGGCTTAGGGGGGGTTTTTTTTTTACCTGAACGGCCCCTTGCGCGGGCTTCCCGGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAgcgc                                                                            <  3:188672/135‑1 (MQ=255)
                                               gcgAGTATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAgcgcg                                                                           <  7:386759/135‑1 (MQ=255)
                                                     tttGCTCCTGGCGAACCCAAAAAAGCCCGGGCTTAGGGCGGTTATTTTTTTAACGTGAACCGGCCGTTGCGGGGGCTGCCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAgg                                                                     <  6:352195/134‑1 (MQ=255)
                                                            cTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGCGGCCCGTCCGTCGAACCGATAGATGTCCCTCCCCAGCCCGCCCCGGGTAAAcc                                                              >  6:133974/1‑135 (MQ=39)
                                                             gggCGTACCGAAAAAAGCCCGGGCTTAGGGGGGGGATTTTTTTTACCTGTCCCGGCCCTTTCGCGGGCTTGCCCGGGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAAccc                                                             <  8:45164/134‑1 (MQ=255)
                                                                 gTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTGGCGGGCTTGGCCCCGGACAGGCCCCTCATGTGCAACGGATTGTGTCCAAGCCCCGCGGGCCCACGGGAAAacccg                                                           >  3:21411/1‑132 (MQ=255)
                                                                 gTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTACATCCCCCGCGCGCCCAGGGTAAACCCCTGa                                                         >  4:232194/1‑135 (MQ=255)
                                                                  tAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGGGTGGCCGGGGAGAGGCCCGTCAGTCGCGCCGCTTGCTGTGCCTGCCCCGCGCGCCCCGGGGTAACCCCGGAt                                                        >  1:244252/1‑135 (MQ=255)
                                                                  tAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTAGCGGGGGTTCCCCCGTATCGGCCCGTCAGTGGGACCCAGACAATTCCCGACCCAGCGCGCCCCGGGGATAccccggct                                                        >  8:20327/1‑131 (MQ=255)
                                                                  tAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTAGCGGGGGTGGCCGCGAAGCGGCCCGTCAGTCGACGCGATTGGTGTCCATCCCCCGCGCGCCCACGGTAAACCCCTGag                                                        >  2:259446/1‑134 (MQ=255)
                                                                    ccaaaaaaaCCCCCGGCTTAGGTGGGGTTTTTTTTTTACCTCTACCGCCCGTTTCGCGGGGTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATg                                                       <  5:69667/131‑1 (MQ=255)
                                                                                           cGGGTATTTTTTTTACCTGTACGGGCCATTTCGGGGGCTGGCGCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGc                               <  5:226146/135‑1 (MQ=255)
                                                                                           cGGGTATTTTTTTTACCGGTACGGCCCGTTTGGGGGGGTTCCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGc                               <  1:61117/135‑1 (MQ=255)
                                                                                           cGGGGATTTTTTTTACCTGAACCGGCCCTGGCGCGGGCTTCCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGc                               <  7:282256/135‑1 (MQ=255)
                                                                                             gggTTTTTTTTTACCTGTACCGCCCCTTCCGGGGGCTTCCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGccc                             <  6:68914/135‑1 (MQ=255)
                                                                                                            tGTACCGGCCCTTGCGCGGGCGTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACa              <  6:454280/135‑1 (MQ=255)
                                                                                                                   gCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGgcag       <  4:153962/135‑1 (MQ=255)
                                                                                                                   gCCCTTCCCCGGGCGGGCCGGCCAAGAGGCCCTCCAGACGAACCGATAGATGTCCAACCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGgcag       <  1:259446/135‑1 (MQ=255)
                                                                                                                          gctgggCCTGCCCGGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAAt  <  5:60893/131‑1 (MQ=255)
                                                                                                                          cgcgGGCTGGCCCGCTAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAAt  <  3:167090/133‑1 (MQ=255)
                                                                                                                                |                                                                                                                                
CAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAAT  >  ACB122/2085354‑2085608

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: