Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 5,096,415 C→G 25.2% intergenic (+28/‑162) CJ019_04792 → / → CJ019_04793 Putative pterin‑4‑alpha‑carbinolaminedehydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,096,4150CG25.2% 81.6 / 7.9 61intergenic (+28/‑162)CJ019_04792/CJ019_04793Putative pterin‑4‑alpha‑carbinolaminedehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/22);  new base G (10/4);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGAC  >  ACB143/5096282‑5096546
                                                                                                                                     |                                                                                                                                   
cACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGccccgt                                                                                                                                   >  2:383215/1‑134 (MQ=255)
   gttGACCTGGTGGAGGCACTCTATCAAAGGCATGAACCGCATCTACTTCTTCATGTGTGCACGCAATTACTTGGTGGCGGAAACGGCTGAAGTCCGGAAGTAAACACCTGAAGGATGGGGCCAGCGGCCCCttttt                                                                                                                                <  6:140997/134‑1 (MQ=255)
   cGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCttttt                                                                                                                                <  2:145011/136‑1 (MQ=255)
   cGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCttttt                                                                                                                                <  2:222485/136‑1 (MQ=255)
   cGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCttttt                                                                                                                                <  4:383989/136‑1 (MQ=255)
   cGTGACCTGGGGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGAACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAAAACCAGAAGGAAGGGGCCAGAGGCCCCttttt                                                                                                                                <  4:233616/136‑1 (MQ=255)
     tGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTc                                                                                                                              <  2:9105/136‑1 (MQ=255)
     tGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCAGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCtttttt                                                                                                                               >  3:387933/1‑135 (MQ=255)
        ccTGGGGGAGCCACTCTATCAACGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACACGGGGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATg                                                                                                                           <  4:10475/136‑1 (MQ=255)
                aGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAATGCCGGAAGTAAACACCAGAAGGAAGGGGTCAGAGGCCCCTTTTTTAATGGGGCatat                                                                                                                   >  6:182005/1‑136 (MQ=255)
                aGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTAATGGGGCatat                                                                                                                   >  2:224992/1‑136 (MQ=255)
                aGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCATTTTTCATGGGGCatat                                                                                                                   >  7:241003/1‑136 (MQ=255)
                    aCTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTAATGGGGCatattttc                                                                                                               >  7:245306/1‑132 (MQ=255)
                       cGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCt                                                                                                            >  6:383413/1‑136 (MQ=255)
                         aTCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATTGGGCATATGTGCGCTGt                                                                                                          >  2:210175/1‑136 (MQ=255)
                         aTCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGAGCTGt                                                                                                          >  2:362163/1‑136 (MQ=255)
                                     cACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGTAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCaa                                                                                                <  6:213057/134‑1 (MQ=255)
                                        cGCAACTAATTAATCATGTGTGCACGCACTTACATGGTGGCGGAAAGGGCTTAAGGCCGTAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGAAAAACCg                                                                                           <  4:193595/136‑1 (MQ=255)
                                             cGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGTTTTTGATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACaaa                                                                                      >  3:32993/1‑136 (MQ=255)
                                             cGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCGTTTTTGATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGCCaa                                                                                       >  8:89454/1‑135 (MQ=255)
                                             cGACTTCATCATGTGTGCACGCACTGACAAGGTGGAGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGACAGAGGCCCGTTTTTTGATGGGGCATATGTGCGGTGTTGGACCGCTAAACCGTCaat                                                                                      >  6:62585/1‑135 (MQ=255)
                                                            tgCACGCACTGAAAAGGTGGCGTAAACGGCTTAAGGCCGTAAGTAAACACCATATGGAAGGGGCATGAGGCCCGTTTTTTCATGGGGATTATGTGCGCTGTTCGCCCGAAAAACCGACAAACGACCGGGaaaaaaa                                                                       <  5:314643/136‑1 (MQ=255)
                                                              cacgcacTTACATGGTGGGTTAAACGGCTTATGGCGGTAAGTAAACACCACACGGACGGGGCCAGAGGCCCGTTTTTTCATGGGGCATATGTGGGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAAcc                                                                     <  5:288219/136‑1 (MQ=255)
                                                              cacgcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGTTTTTTCATGGGGCATATGTTCGCTGTTCCCCCGCAACACCGACAAACGACCGGGAAAAAAAc                                                                      >  2:384415/1‑135 (MQ=255)
                                                              cacgcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGTTTTTTAATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAAc                                                                      >  3:340437/1‑135 (MQ=255)
                                                               acgcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTAATTGGGCATATGTGCGCTGTTCGCCCGAAAAACCGAAAAACGACCGGGaaaaaaaac                                                                     >  1:114553/1‑133 (MQ=255)
                                                                 gcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAACCCGACAAACGACCGGGAAAAAAACCaa                                                                   >  7:270811/1‑135 (MQ=255)
                                                                      gACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGt                                                              >  8:12506/1‑135 (MQ=255)
                                                                      gACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACAGACAAACGACCGGGAAAAAAACAAAACTGt                                                              >  2:79586/1‑135 (MQ=255)
                                                                      gACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTAATGGGGCAAATTTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAAACAAACTGt                                                              >  6:69399/1‑135 (MQ=255)
                                                                       aCAAGGTGGCGTAAAGGGCTGAAGGCCGGAAGTAAACACCAGAAGTAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTc                                                             <  7:121491/135‑1 (MQ=255)
                                                                        cAAGGTGGCGTAACCGGCTGAAGGCCGGAAGTAAACACCAGAAGTAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCAt                                                           <  3:390847/136‑1 (MQ=255)
                                                                        cAAGGTGCCTGAAACGGTTTAAGGCCGTAATTAAACACCAGAAGGAAGGGGCCAGAGGCCCGTTTTTTCATGGGGCATATGTGCTCTGTTGGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCAt                                                           <  3:59446/136‑1 (MQ=255)
                                                                        cAAGGGGTCGTAAACGGCTGAAGGCGGTATTTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATTTGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCAt                                                           <  7:218275/136‑1 (MQ=255)
                                                                        aaaGGTGGCGTAACCGGCTCAAGCCCGTAAGAAAACACCAGACGGACGGGGCCAGCGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGAAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCAt                                                           <  7:338725/135‑1 (MQ=255)
                                                                         aaGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCTAGAGGCCCGTTTTTTAATGGGGCATATGTGCGCTGTTCGCCCGCCAAACCGCCAAACGACCCGGAAAAAAACCAAACTGTCAta                                                          >  7:338045/1‑135 (MQ=255)
                                                                          cgttGGCGTAAACGGATGAAGGCCGTATGAAAACACATCAAGGACGGGCCCAGATCCCCTTTTTTTAAGGGGTCATATGTGCGCTGTTCGCCCGAAAAACCCACAAACGACCGGGAAAAAAACCAAACTGTCATcc                                                         <  4:81288/133‑1 (MQ=255)
                                                                          cggTGGCGTAAACGGCTGATGGCCGGAAGTAAACACCAGACGGCCGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTGCGCCCGAAAACCCGACAAACGACCGGGAAAAAAACCAAACTGTCATcc                                                         <  3:215793/135‑1 (MQ=255)
                                                                          aGGTGGCGGAACCGGGTGATGGCGGTAAGTAAACACCACACGCCGGGGGCCAGAGGCCCTTTTTTTCATGGGGCATATGTGCGCTGTGCGCCCGCAAACCCCACAAACGACCGGGAAAAAAACCAAACTGTCATcc                                                         <  1:132681/136‑1 (MQ=255)
                                                                          aGGTGGCGGAAACGGCTGAAGGCCGGAATTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATcc                                                         <  1:49399/136‑1 (MQ=255)
                                                                          aGGGGGCGGAAACGGCGGAAGGCCGGAAGTAACCACCAGAAGTAAGGGGCCAGAGGCGCCTTTTTTCATGTGGCATATGTGCGCTGGTCGCCGGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATcc                                                         <  6:267621/136‑1 (MQ=255)
                                                                           gggggCGTAAACGGCTGACGGCGGTAAGAAAACACCACACGCATGGGGCCCGCGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCCAAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCa                                                        <  8:164239/133‑1 (MQ=255)
                                                                              ggCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGAGCGCTGTTCGCCCGCAAAACCGACAAACGCCAGGGAAAAAAACCAAACTGTCAACCacac                                                     >  8:160353/1‑136 (MQ=255)
                                                                                  gAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTtgt                                                 >  8:139473/1‑136 (MQ=255)
                                                                                  gAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCttttt                                                                                                                                >  5:372592/1‑57 (MQ=255)
                                                                                  gAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCttttt                                                                                                                                <  6:372592/57‑1 (MQ=255)
                                                                                   aaaCGGCTGACGGCGGGATGTAAACACCATATGGAAGGGGCCAGCGGCCCCTTTTTTCATGGGGAATATGTGCGCGTTTCGCCCGAAAACCCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTtgtg                                                <  1:267283/136‑1 (MQ=255)
                                                                                    aaCGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACACACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTgtgt                                               >  6:171199/1‑136 (MQ=255)
                                                                                               gCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGTTTTTGATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAACAACCAAACTGTCATCCAGACTTGTGTATCCTGccc                                      >  6:319157/1‑134 (MQ=255)
                                                                                                     aGTAACCCCCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCGGTGCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCtttt                              <  7:79527/136‑1 (MQ=255)
                                                                                                            acCAGAAGGAAGGGGCCAGAGGCCCGTTTTTTAAAGGGGCATATTTGCGCTGTTGGCCCCCCAAACCGACAAACGACCGGGAAAAAAACAAAACAGTCAAACAGACTTGTTTATACTGCCCGCGCGTTTTTaacg                        >  6:20392/1‑133 (MQ=255)
                                                                                                                 aaggaaggGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCCAACTGTCATCCAGACTTGTGTATCCGGCCCGCGCTTTTTGAAGCTACaa                   >  7:179319/1‑135 (MQ=255)
                                                                                                                   ggacggGGCCCGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAAcg                <  3:304047/136‑1 (MQ=255)
                                                                                                                   ggacggGGCCCGAGCCCCCTTTTTTCTTGGGGCATATGTGCGCGGTGCCCCCGCAAACCCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAAcg                <  4:318610/136‑1 (MQ=255)
                                                                                                                    cacgggcCCACAGCCCCATTTTTTATTGGGGCATATGGGCGCTGTCCGCCCGCAAACCCCAAAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAAcgc               <  5:111594/133‑1 (MQ=255)
                                                                                                                         ggCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGCCCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGccgg          >  2:39088/1‑134 (MQ=255)
                                                                                                                         ggCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATcc                                                         <  7:243698/89‑1 (MQ=255)
                                                                                                                         ggCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATcc                                                         >  8:243698/1‑89 (MQ=255)
                                                                                                                         ggCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGCGCTTTTTGAAGCGACAAACGCGCCTg          >  2:65699/1‑136 (MQ=255)
                                                                                                                           ccAGAGGCCCGGTTTTTGATGGGGCATATTTGCGGGGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCGGgc        >  4:317096/1‑136 (MQ=255)
                                                                                                                           ccAGAGGCCCCTTTTTTTATGGGGCATATGTGCGCAGTTCGCCCGCCAAAACGACAAACGACCGGGAAAAAAACCAAACTGGCATCCAGAAGTGTGGATACTGCCCGAGCGTTTTGATGAGTCAAACGCGCCAGgc        >  2:323226/1‑136 (MQ=255)
                                                                                                                                 gCCCCGTTTTTCATGGGGCATATGTGCGCGGTTCGCCCGCAGAACCGCCAAACGCCCCGGAAAAAAACCAAACTGTCATCCCGACTTGTTATTCCTGCCCGACCGTTTTGAAGCTACAAAATCGCCGGGCGCAgcc  >  8:308857/1‑134 (MQ=255)
                                                                                                                                     |                                                                                                                                   
CACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGAC  >  ACB143/5096282‑5096546

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