Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,491,319 A→C 35.7% intergenic (‑681/+40) CJ019_03306 ← / ← CJ019_03307 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,491,3190AC35.7% 10.0 / 2.6 31intergenic (‑681/+40)CJ019_03306/CJ019_03307hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (5/13);  new base C (9/2);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 7.78e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

GCCCGCTCCTACATCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCC  >  ACB122/3491194‑3491453
                                                                                                                             |                                                                                                                                      
gCCCGCTCCTACATCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCGGCCTGGCcg                                                                                                                               >  6:174743/1‑135 (MQ=255)
       ccaaCATCGACCGGGGGTGCCGCCACGTGGTCGGAGCTGGCTAACCCGCGAATAGGCCGGACCGGCCGCCCACTACCCCTGGTCATCCCCCGACACCCCCGCCCCCCCATTCGCGGCCAGCCTGCCCGCGAacac                                                                                                                        <  6:278513/132‑1 (MQ=19)
             tCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGGCGCCGGCCCGCGCAcccccccgc                                                                                                                  >  3:211312/1‑130 (MQ=21)
              cGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCGGCCCCCGAACACCCGCGCa                                                                                                                 >  7:241269/1‑135 (MQ=255)
                   ggtggtGTCGCCACGTTGTAGGCGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCCCCGTCACAGCCGCCCCCGCATTCGCGGGCCGCCGGCCCGCGCACCCCCGGGCCGCCgg                                                                                                            >  6:29353/1‑135 (MQ=18)
                   ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCGCTATCCCTGGTCATTCACCGTCACAGCCGCCCCCGCATTAGCCGGCCGCCGGGCCGCGCCCACCCGCgccgcccg                                                                                                            >  7:339651/1‑133 (MQ=18)
                   ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGGCGCCGGCCCGCCAACACCGgcgccccccg                                                                                                            >  4:385874/1‑129 (MQ=21)
                   ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGGCGCCGGCCCCCGAACCCCCGCGCCGCCgg                                                                                                            >  2:138148/1‑135 (MQ=21)
                   ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGGAGCCAGCCCGCGCACACCGGCGCAGCCgg                                                                                                            >  1:153419/1‑135 (MQ=255)
                   ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCa                                                                                                                                        >  8:144402/1‑107 (MQ=255)
                   ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCa                                                                                                                                        <  7:144402/107‑1 (MQ=255)
                   ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCGGCCCGCGAACACCGGCGCAGCCgg                                                                                                            >  4:110454/1‑135 (MQ=255)
                   ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCCGGCAGCCTGCCCGCCAACACCGGCCCCGCCgg                                                                                                            >  4:94133/1‑135 (MQ=255)
                        tgtCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCCGCCCCCGCACACCCGCCCCGCCGGGTTCa                                                                                                       >  3:342768/1‑135 (MQ=18)
                            gCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCGCGCACACCGGCGCCGCCGGGTTCCTCCa                                                                                                   >  2:24002/1‑135 (MQ=21)
                                          gcgggTTTCCCCGCTAATCGGCCGGACCGGCCGGCAACTAGCCCTGTGCATGCCCCGTCCCCCCCGCCCCCCCATTCGGGGGCCGCCGGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATc                                                                                     <  2:45145/132‑1 (MQ=255)
                                           cTGGTTTACCCGCGAATAGGCCGTACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCg                                                                                    <  4:401878/135‑1 (MQ=255)
                                             gttttCCCCGCTATTAGGCCGGCCCTGCCGGCCACTATCCGTGGTCATTCACGGCCCCAGCCCCCACCGCATTCGGGGGCAGCGTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGAc                                                                                  <  6:457540/133‑1 (MQ=255)
                                             ggTTTACCCGCTAATTGGCCGGACCTCCCGCCAACTTTCCCTGGCAATTCACCGTCACTCCCGCCACCGCTTTGGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGAc                                                                                  <  3:313995/135‑1 (MQ=255)
                                             ggTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCTGCCCGCGCACACCGGCGCCGCCGGGTTCATGCAACTCAACGccccacgcc                                                                                  >  7:240379/1‑129 (MQ=21)
                                                   cccGCTAATCGGCCGGCCCTGCCGGCAAATAGCCTTGGTATTCCCCCGACCCCGCCGCCACCGCATGCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTAc                                                                            <  7:258106/135‑1 (MQ=255)
                                                     cgggattAGGCCGGACGTCCCGGCCAGTTGACCTGGGCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCACCGTGCCCGGGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGc                                                                          <  8:116419/129‑1 (MQ=37)
                                                         aaTAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCTGCCCGCGAACACCGGCCCCGCCGGTTTCATACAACTCAACGACCCACGACCGGTACGCCAAc                                                                      >  6:347505/1‑135 (MQ=25)
                                                                    ccTGCCGGCCAATACCCCTGGTCATACCCCGCAACCCCCCCCACCGCATGCGCGGGCAGCCTCCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCAt                                                           <  1:150551/135‑1 (MQ=255)
                                                                             cACTAACCCTGGACATTCACCGACACAGCGGCCCCCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCt                                                  <  4:69363/135‑1 (MQ=255)
                                                                              aCTACCCCTGGCAATCCCCCGCCACCCCCGCCCCCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCtt                                                 <  5:174743/135‑1 (MQ=255)
                                                                                               cccGTCACCGCCGCCCCCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCg                                 <  4:369363/133‑1 (MQ=255)
                                                                                               aCCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGGTTCATCCAACGAAACGCCCAACGCCCTGTACGCAAACGTCTCTGCCATCGCTTGCCGTGCATGGCCAAATGTCgg                                >  2:365674/1‑134 (MQ=255)
                                                                                                              aCCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCg                 <  8:241269/135‑1 (MQ=255)
                                                                                                                      cgcgGGCCCCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGaa          <  3:385874/134‑1 (MQ=255)
                                                                                                                       gcgGGCTGCCTGCCGGCCAACACCGGCGCAGCCGGTTTCATCCAACTCCACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGc        <  5:164543/135‑1 (MQ=255)
                                                                                                                        cgGGCAGCGTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGca       <  3:94133/135‑1 (MQ=255)
                                                                                                                        cgGGCAGCCTCCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGGCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGca       <  7:362725/135‑1 (MQ=255)
                                                                                                                             aGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACAccc  <  4:110403/135‑1 (MQ=255)
                                                                                                                             |                                                                                                                                      
GCCCGCTCCTACATCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCC  >  ACB122/3491194‑3491453

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: