Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 96,480 T→G 41.2% intergenic (+121/‑72) cynT → / → CJ019_00104 Carbonic anhydrase 1/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT244096,4800TG41.2% 15.4 / 5.6 34intergenic (+121/‑72)cynT/CJ019_00104Carbonic anhydrase 1/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/16);  new base G (8/6);  total (12/22)
Fisher's exact test for biased strand distribution p-value = 2.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTACTAAGCGACCCGTCCGCCGCTAGAACCTTGATAGCCGGCTGCACCCTTCAGGGGTGTGGCGCGGCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCTGCGCGCTTGTCAGGGGCAAACGTGCCCACGGCCAGGGGAATGGCCGTGCGACTAGAGAAAGGAGAAACACGGTGAACATTACACAGTTGAAAGCGGCCCTGCCGCGTGAGTTGCTGGCGTCGGTGGTG  >  KT2440/96353‑96608
                                                                                                                               |                                                                                                                                
ctactaAGCGACCCGTCCGCCGCTAGAACCTTGATAGCCGGCTGCACCCTTCAGGGGTGTGGCGCGGCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCGgcgcgcg                                                                                                                           >  8:227510/1‑134 (MQ=255)
ctactaAGCGACCCGTCCGCCGCTAGAACCTTGATAGCCGGCTGCACCCTTCAGGGGTGTGGCGCGGCCGTCGGCGGCCGTGGATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCCCCGgcccgcg                                                                                                                           >  5:132775/1‑134 (MQ=255)
 tactaAGCGACCCGTCCGCCGCTAGAACCTTGATAGCCGGCTGCACCCTTCAGGGGTGTGGCGCGGCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCTGCGCGCtt                                                                                                                          <  6:156377/135‑1 (MQ=255)
 tactaAGCGACCCGTCCGCCGCTAGAACCTTGATAGCCGGCTGCACCCTTCAGGGGTGTGGCGCGGCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCTGCGCGCtt                                                                                                                          <  7:363413/135‑1 (MQ=255)
 tactaAGCGACCCGTCCGCCGCTAGAACCTTGATAGCCGGCTGCACCCTTCAGGGGTGTGGCGCGGCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCTGCGCGCtt                                                                                                                          <  5:293368/135‑1 (MQ=255)
 tactaAGCGACCCGTCCGCCGCTAGAACCTTGATAGCCGGCTGCACCCTTCAGGGGTGTGGCGCGGCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCGgcgcgcgt                                                                                                                          >  2:17688/1‑133 (MQ=255)
 tactaAGCGACCCGTCCGCCGCTAGAACCTTGATAGCCGGCTGCACCCTTCAGGGGTGTGGCGCGGCCGTCGGCGGGCATGAATTACCTGCCAGGCCAGCCGGGGGGGGCGCGGGCGGCCCCCGCCCgggcggcg                                                                                                                          >  1:335275/1‑132 (MQ=255)
      ccctaCCCGTCCGCCGCTAGAACCTTGATAGCCTGCTGCACCCTTCAGGGGTGTGGCGCGGCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCTGCGCGCTTGTCAg                                                                                                                     <  1:159234/131‑1 (MQ=255)
          aCCCGTCCGCCGCTAGAACCTTGATAGCCGGCTGCACCCTTCAGGGGTGTGGCGCGGCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCTGCGCGCTTGTCAGGGGc                                                                                                                 <  7:248170/135‑1 (MQ=255)
               tCCGCCGCTAGAACCTTGATAGCCGGCTGCACCCTTCAGGGGTGTGGCGCGGCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCTGCGCGCTTGGCAGGGGCAAACGt                                                                                                           <  3:285355/136‑1 (MQ=255)
                                  tcgCCGGCTGCACCCTTCAGGGGTGGGGCGCGGCCTGCGGCGGCCTTGAATGCCCGGACTGCCCTGCCGGCGGGCGCGCTGGCGGCCCGCGCCTGCGCGCGTGTCCGGGGCAAACGTGCCCACGGCCAGGGGAAt                                                                                         <  1:215225/133‑1 (MQ=255)
                                   aGCCGGCTGCACCCTTCAGGGGTGTGGCGCGGCCGTCGGCTGCCTTGAATTACCTGACTGCCCTGCCGGCGTGCGCGCTGGCGGCCCGCCCCGGCGCGCGTGGCAGGGGGAAAAGGGCCCCCGGCCACGGGaaat                                                                                        >  5:135151/1‑133 (MQ=255)
                                     ccGGCTGCACCCTTCAGGGGGGTGGCGCGGCCTTCGGCTGCCTTGAATTCCCTGAGTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCTGCGCGCTTGTCAGGGGCAAACGTGCCCACGGCCAGGGGAATgg                                                                                       <  7:237855/134‑1 (MQ=255)
                                                 tgccggggggTGGCGCGGCATTCGGCTGGATTGATTTCCCTGACTGCCCTGCGGGGTTGCTGGCTGGCGGCCCGCGCCGGCGGTCTTGTCTGGGGCAAACTTGCCCACGGCCAGGGGAATGGCCGTGCGACTaga                                                                          <  7:345063/131‑1 (MQ=255)
                                                          gtgGGGGGGCCTTCGGCTGCTTTTATTTCCCTGAGTGCCTGGCGGGGGGGCGGGGGGGGGGCCGGGGCGGGGGGGTTTTTCCGGGGAATACGTGCCCACGGCCAGGGGAATGGCCGTGCGACTAGAGAAAGgaga                                                                 <  6:227284/135‑1 (MQ=255)
                                                          gtgGCGCGGCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCTGCGCGCTTGTCAGGGGCAAACGTGCCCACGGCCAGGGGAATGGCCGTGCGACTAGAGAAAGgaga                                                                 <  8:229638/135‑1 (MQ=255)
                                                              cggggCCTTCGGGTGCCTTGAATTCCCTGACGGCCGGGCCGGCGGGCGCGCGGGCGGCCCGCGCCTGCGCGCTTTCCCGGGGCAAACGTCCCCACGGCCAGGGGAATGGCCGTGCGACTAGAGAAAGGAGAAacac                                                            <  4:12194/133‑1 (MQ=255)
                                                              cgcgGCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCTGCGCGCTTGTCAGGGGCAAACGTGCCCACGGCCAGGGGAATGGCCGTGCGACTAGAGAAAGGAGAAaca                                                             >  7:208962/1‑135 (MQ=255)
                                                              cgcgGCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCGGCGCGCTTGTCAGGGGCAAACGTGCCCACGGCCAGGGGCATGGCCGTGGGGCTAGAGAAAGGAGAAacac                                                            >  3:352198/1‑136 (MQ=255)
                                                              cgcgGCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCGGCGCGCGTGTCAGGGGCAAACGTGCCCACGGCCAGGGGAATGGCCGGGGGGCGAGAGAAAGGAGAAaccc                                                            >  3:268322/1‑134 (MQ=255)
                                                              cgcgGCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCGGCGCGCGTGGGCGGGGGCAACGTGCCCACGGCCAGGGGGATGGCCGGGGGGCGAGAGAAAGGGGAAaca                                                             >  5:184673/1‑135 (MQ=255)
                                                              cgcgGCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCCCCTGCGCGCTTGTCAGGGGCAAACGTGCCCACGGCCAGGGGAATGGCCGTGCGGCTAGAGAAAGGAGAAacaa                                                            >  3:53313/1‑135 (MQ=255)
                                                               gcgGCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCTGCGCGCTTGTCAGGGGCAAACGTGCCCACGGCCAGGGGAATGGCCGTGCGACTAGAGAAAGGAGAAacac                                                            <  7:227510/135‑1 (MQ=255)
                                                               gcgGCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCGGGCGGCCCGCGCCTGCGCGCTTGTCAGGGGCAAACGTGCCCACGGCCAGGGGAATGGCCGTGCGACTAGAGAAAGGAGAAacac                                                            <  7:342398/135‑1 (MQ=255)
                                                                 ggCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCTGCGCGCTTGTCAGGGGCAAACGTGCCCACGGCCAGGGGAATGGCCGTGCGACTAGAGAAAGGAGAAACACGGt                                                         <  4:273924/136‑1 (MQ=255)
                                                                   ccTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCGGCGCGCGTGTCAGGGGCAAACGTGCCCACGGCCCGGGGTATGGCCGGGCGACTAGAGAAAGGAGCAACACgggg                                                        >  8:215139/1‑133 (MQ=255)
                                                                        ggCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCGGGCGGCCCGCGCCGGCGCGCTTGTCAGGGGCAAACGTGCCCACGGCCAGGGGAATGGCCGTGCGACTAGAGAAAGGAGAAACACGGTGAACAtt                                                  <  3:300640/136‑1 (MQ=255)
                                                                         gCTGCGTTGAATGCCCTGAGGGCCCGGCGGGGGGGCGCGCGGGCGGCCCGCGCCGGCGGGCTTGGCTGGGCCAAAGGGGCCCACGGCCAGGGGAATGGCCGTGCGACTAGAGAAAGGAGAAACACGGTGAACAtt                                                  <  1:171377/135‑1 (MQ=255)
                                                                         gCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCTGCGCGGTTGGCAGGGGCAAACGTGCCCACGGCCAGGGGAATGGCCGTGCGACTAGAGAAAGGAGAAACACGGTGAACAtt                                                  <  6:87614/135‑1 (MQ=255)
                                                                         gCTGCCTTGAATTCCCTGACCGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCGTGCGCGCTTGTCAGGGTCAAACGTGCCCACGGCCAGGGGAATGGCCGTGCGACTAGAGAAAGGAGAAACACGGTGAACAtt                                                  <  6:170139/135‑1 (MQ=255)
                                                                                    ttCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCTGCGCGCTTGTCAGGGGCAAACGTGCCCACGGCCAGGGGAATGGCCGTGCGACTAGAGAAAGGAGAAACACGGTGAACATTACACAGTTGa                                        <  8:208962/134‑1 (MQ=255)
                                                                                        cTGACGGCCCGGCCGGGGGGCGCGCGGGGGGCCCGCGCCGGCGCGCTGGTCCGGGGAAACCGGGCCCACGGCCAGGGGAATGGCCGTGCGACTAGAGAAAGGAGAAACACGGTGAACATTACACAGTTGAAAGCgg                                  <  4:268322/136‑1 (MQ=255)
                                                                                                                       ccggcgcCGGCGCGCTTGTCGGGGCCAAACGGGCCCACGGCCTGGGGAATGGCCGTGCGACTAGAGAAAGGAGAAACCCGGTGAACATTACACAGTTGAAAGCGGCCCTGCCGCGTGAGTTGCTGGCGTCggtggt   <  3:304602/133‑1 (MQ=255)
                                                                                                                         cgcgCCTGCGCGCTTGTCAGGGGCAAACGTGCCCACGGCCAGGGGAATGGCCGTGCGACTAGAGAAAGGAGAAACACGGTGAACATTACACAGTTGAAAGCGGCCCTGCCGCGTGAGTTGCTGGCGTCggtggtg  >  1:69182/1‑135 (MQ=255)
                                                                                                                               |                                                                                                                                
CTACTAAGCGACCCGTCCGCCGCTAGAACCTTGATAGCCGGCTGCACCCTTCAGGGGTGTGGCGCGGCCTTCGGCTGCCTTGAATTCCCTGACTGCCCTGCCGGCGTGCTCGCTGGCGGCCCGCGCCTGCGCGCTTGTCAGGGGCAAACGTGCCCACGGCCAGGGGAATGGCCGTGCGACTAGAGAAAGGAGAAACACGGTGAACATTACACAGTTGAAAGCGGCCCTGCCGCGTGAGTTGCTGGCGTCGGTGGTG  >  KT2440/96353‑96608

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: