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

CAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCC‑GCGAAGAGGCCCTTCA‑GTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAAT  >  ACB122/2085354‑2085608
                                                                                                                                |                                                                                                                                
cAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTAGCGGGCGTGccc                                                                                                                            >  5:105577/1‑135 (MQ=255)
    gACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCGGGCGTAACGAAAAAATCCCCGGCTTAGGGCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTGGCCC‑GCGa                                                                                                                       <  4:382735/135‑1 (MQ=255)
    gACACTCGCGGGGCGGGCAAGCCGCTGAAGCTGCGTCGCTACCGCGCGTATTGATCCTGGCGAAACCAAAAAATCCCCGGCTTAGGGCGGGTATTTTTTTTACCTGACCCGCCCCTTCCGCGGGCTTCCCC‑GCGa                                                                                                                       <  6:380997/135‑1 (MQ=255)
     acacTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCc‑gcgc                                                                                                                       >  8:24029/1‑133 (MQ=255)
     acacTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCC‑GCGaa                                                                                                                      >  6:123110/1‑135 (MQ=255)
           gcgTGGCTGGCAAGCAGCTGAAGATGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTAACGGCCATTTAGCCGGCTTTCCC‑GCGAAGCGGcc                                                                                                                >  6:60893/1‑135 (MQ=255)
            cgTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGGGTGCCC‑GCGAAGAGGccc                                                                                                               >  7:163922/1‑135 (MQ=255)
            cgTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGGCC‑GCGCAGAGGgcc                                                                                                               >  1:308195/1‑135 (MQ=255)
                  gggcAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACCAAAAAATCCCCGGCTTAGGGCGGGTATTTTTTTTACCTGTACCGGCCCTTGCGCGGGCTTGCCC‑GCGAAGAGGCCCTTCA‑Gt                                                                                                        <  1:296081/134‑1 (MQ=255)
                    gcAAGCAGCTGAAGCTGGGTCGCTACCGCGAGTATTGAGCGTGGCTTAACTAAAAAATCCCCGGGTTAGGTCGGGTTTTTTTTTTACCTGTACCGGCCCTTTCGCGGGGTTGCCC‑GGGAAGAGGCCCTTCA‑GTCg                                                                                                      <  8:62443/135‑1 (MQ=255)
                    gcAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTAGCGGGGGTGGCC‑GCGCAGGGGCCCGTGA‑Gtg                                                                                                       >  6:360386/1‑133 (MQ=255)
                          agcTGAAGCTGCGGCGCTACCGCGAGTATTTATCCGGGCGTAACTAAAAAATCCCGGGCTTCGGGCGGGTATTTTTTTTACCTGTACCGCCCCTTTCGCGGGCTGGCCC‑GCGAAGAGGCCCTTCA‑GTCGAACCGa                                                                                                <  5:123110/135‑1 (MQ=255)
                                aGCTGCGTGGCTTCCGCGAGTATTGAGCCTGGCTTAACGAAAAAATCCCCTGCTTAGGGCGGGTTTTTTTATTACCTGACCCGGCCCTTGCCGGGGGGCCCCG‑GGGAAGAGGCCCTTCA‑GTCGAACCGATAGAt                                                                                           <  6:214984/134‑1 (MQ=255)
                                aGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTACCGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTCCCGGCCCTTTCGCGGGCTTGCCC‑GCGAAGAGGCCCTTCA‑GTCGAACCGATAGAt                                                                                           <  7:267291/134‑1 (MQ=255)
                                         gCTACCGCGACGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCC‑GCGAACAGGCCCATCA‑GTCGAACCGAAAGATGTCGATcccc                                                                                 >  8:346742/1‑135 (MQ=39)
                                              cgcgAGGATTGACCCGGGCGTACCGAAAAAATCCCCGGCTTAGGGGGGGTTTTTTTTTTACCTGAACGGCCCCTTGCGCGGGCTTCCCG‑GCGAAGAGGCCCTTCA‑GTCGAACCGATAGATGTCCATCCCCAgcgc                                                                            <  3:188672/135‑1 (MQ=255)
                                               gcgAGTATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCC‑GCGAAGAGGCCCTTCA‑GTCGAACCGATAGATGTCCATCCCCAgcgcg                                                                           <  7:386759/135‑1 (MQ=255)
                                                     tttGCTCCTGGCGAACCCAAAAAAGCCCGGGCTTAGGGCGGTTATTTTTTTAACGTGAACCGGCCGTTGCGGGGGCTGCCCC‑GCGAAGAGGCCCTTCA‑GTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAgg                                                                     <  6:352195/134‑1 (MQ=255)
                                                            cTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCC‑GCGAAGCGGCCCGTCC‑GTCGAACCGATAGATGTCCCTCCCCAGCCCGCCCCGGGTAAAcc                                                              >  6:133974/1‑135 (MQ=39)
                                                             gggCGTACCGAAAAAAGCCCGGGCTTAGGGGGGGGATTTTTTTTACCTGTCCCGGCCCTTTCGCGGGCTTGCCC‑GGGAAGAGGCCCTTCA‑GTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAAccc                                                             <  8:45164/134‑1 (MQ=255)
                                                                 gTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTGGCGGGCTTGGCC‑CCGGACAGGCCCCTCATGTGCAACGGATTG‑TGTCCAAGCCCCGCGGGCCCACGGGAAAacccg                                                           >  3:21411/1‑132 (MQ=255)
                                                                 gTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCC‑GCGAAGAGGCCCTTCA‑GTCGAACCGATAGATGTACATCCCCCGCGCGCCCAGGGTAAACCCCTGa                                                         >  4:232194/1‑135 (MQ=255)
                                                                  tAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGGGTGGCC‑GGGGAGAGGCCCGTCA‑GTCGCGCCGCTTGCTGTGCCTGCCCCGCGCGCCCCGGGGTAACCCCGGAt                                                        >  1:244252/1‑135 (MQ=255)
                                                                  tAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTAGCGGGGGTTCCC‑CCGTATCGGCCCGTCA‑GTGGGACCCAGACAATTCCCGACCCAGCGCGCCCCGGGGATAccccggct                                                        >  8:20327/1‑131 (MQ=255)
                                                                  tAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTAGCGGGGGTGGCC‑GCGAAGCGGCCCGTCA‑GTCGACGCGATTGGTGTCCATCCCCCGCGCGCCCACGGTAAACCCCTGag                                                        >  2:259446/1‑134 (MQ=255)
                                                                    ccaaaaaaaCCCCCGGCTTAGGTGGGGTTTTTTTTTTACCTCTACCGCCCGTTTCGCGGGGTTGCCC‑GCGAAGAGGCCCTTCA‑GTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATg                                                       <  5:69667/131‑1 (MQ=255)
                                                                                           cGGGTATTTTTTTTACCTGTACGGGCCATTTCGGGGGCTGGCGC‑GCGAAGAGGCCCTTCA‑GTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGc                               <  5:226146/135‑1 (MQ=255)
                                                                                           cGGGTATTTTTTTTACCGGTACGGCCCGTTTGGGGGGGTTCCCC‑GCGAAGAGGCCCTTCA‑GTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGc                               <  1:61117/135‑1 (MQ=255)
                                                                                           cGGGGATTTTTTTTACCTGAACCGGCCCTGGCGCGGGCTTCCCC‑GCGAAGAGGCCCTTCA‑GTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGc                               <  7:282256/135‑1 (MQ=255)
                                                                                             gggTTTTTTTTTACCTGTACCGCCCCTTCCGGGGGCTTCCCC‑GCGAAGAGGCCCTTCA‑GTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGccc                             <  6:68914/135‑1 (MQ=255)
                                                                                                            tGTACCGGCCCTTGCGCGGGCGTGCCC‑GCGAAGAGGCCCTTCA‑GTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACa              <  6:454280/135‑1 (MQ=255)
                                                                                                                   gCCCTTTCGCGGGCTTGCCC‑GCGAAGAGGCCCTTCA‑GTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGgcag       <  4:153962/135‑1 (MQ=255)
                                                                                                                   gCCCTTCCCCGGGCGGGCCG‑GCCAAGAGGCCCTCCA‑GACGAACCGATAGATGTCCAACCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGgcag       <  1:259446/135‑1 (MQ=255)
                                                                                                                          gctgggCCTGCCCGGCGAAGAGGCCCTTCA‑GTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAAt  <  5:60893/131‑1 (MQ=255)
                                                                                                                          cgcgGGCTGGCCC‑GCTAAGAGGCCCTTCA‑GTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAAt  <  3:167090/133‑1 (MQ=255)
                                                                                                                                |                                                                                                                                
CAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCC‑GCGAAGAGGCCCTTCA‑GTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAAT  >  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: ‑