Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,624,111 C→A 16.4% intergenic (+117/+181) CJ019_04367 → / ← pgsA hypothetical protein/CDP‑diacylglycerol‑‑glycerol‑3‑phosphate 3‑phosphatidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,624,1110CA16.4% 70.3 / 2.7 34intergenic (+117/+181)CJ019_04367/pgsAhypothetical protein/CDP‑diacylglycerol‑‑glycerol‑3‑phosphate 3‑phosphatidyltransferase
Reads supporting (aligned to +/- strand):  ref base C (18/9);  new base A (2/4);  total (21/13)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTG  >  ACB143/4623977‑4624241
                                                                                                                                      |                                                                                                                                  
aGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCTTTGTGGTGTTTGTACTTGAATTTCGTAGATTGCGGTGTATTTAGGCATAATGCATTCTAACGAAAAAAGGCCCCGAAGGGCCTTTTTTc                                                                                                                                    <  3:82587/135‑1 (MQ=255)
aGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTc                                                                                                                                    >  1:122847/1‑135 (MQ=255)
 gCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCaa                                                                                                                                  <  3:65132/136‑1 (MQ=255)
        ttCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTTTGGTGTTTGTAGTTTAATTTCGTAGAGGGCGGGGTTTTGAGGCAGAATGCATTCGAACGAAAACAGGCCCCGAAGGGCCTTTTTTCAACGACtt                                                                                                                            <  2:269890/134‑1 (MQ=255)
          ccAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTa                                                                                                                           >  6:280416/1‑134 (MQ=255)
            aGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAACGCATTCGAACGAAAAAAGGCCCCGACGGGCCTTTTTTCAACGACTTACAGa                                                                                                                       <  1:5356/136‑1 (MQ=255)
               tAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACtt                                                                                                                    >  3:260239/1‑136 (MQ=255)
                    cACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACt                                                                                                               >  1:17213/1‑136 (MQ=255)
                         cGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAAc                                                                                                          >  2:89814/1‑136 (MQ=255)
                           tCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCt                                                                                                        <  1:82475/136‑1 (MQ=255)
                               ggCCAGGGCGCCGCGTTGTGGTGGTTGTACTTTAATTTCGTAGAGGGGGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGTAGGGCCTTTTTTGAACGACATACACACTTCGACTGACACCTg                                                                                                       >  3:32036/1‑133 (MQ=255)
                                            cgTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTAGACTGAAACCTGTAGATCTTTATTTg                                                                                         >  5:175362/1‑134 (MQ=255)
                                                         gTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACgaga                                                                          >  5:133350/1‑136 (MQ=255)
                                                         gTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCAGAAGGGCCTTTTTTCAACGACTTACAGACGTCGCCTGAAACATGTAGATCTTTATTTGGCGCTGGACACtagc                                                                          >  5:188259/1‑135 (MQ=255)
                                                                 aTTTCGTAGAGGGCGGGGTATTCAGGCAGAATGCATTCGAACGAAAAAAGGCACCGAAGGGCCTTTTTTAAACGACTTACAGACTTCGACTGAACCCTGAAGATCTTTATTTGGAGCGGGAAACGCGACTCGAACt                                                                  >  6:171564/1‑136 (MQ=255)
                                                                               ggggTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGATACGAGACTCGAACTCGCGACCCCGACCt                                                    >  5:234019/1‑136 (MQ=255)
                                                                               ggggTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCt                                                    >  4:247038/1‑136 (MQ=255)
                                                                                    aTTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGc                                                >  6:8392/1‑135 (MQ=255)
                                                                                        aGGCAGAATGAATGCGAACTAAAAACGCCCCCGAAGGGCTTTTTTTAAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGt                                           <  8:213489/136‑1 (MQ=25)
                                                                                             gattGCATTCTAACTAAAAATGGCCCCTATGGGCTTTTTTTAAACGACTTACAGACTTCGACTGAAACCTGTATATATTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGc                                      <  5:104014/133‑1 (MQ=21)
                                                                                             gAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGc                                      <  3:247038/136‑1 (MQ=255)
                                                                                                tGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAg                                                                                      <  7:276621/85‑1 (MQ=255)
                                                                                                tGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAg                                                                                      >  8:276621/1‑85 (MQ=255)
                                                                                                tGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGTAAACGAGACTCGAACTCGCGACCCCGACATTGGCAAGGTCGTGctc                                    >  6:112212/1‑135 (MQ=255)
                                                                                                tGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGCCCTTGGCAAGGTCGTGctc                                    >  1:58265/1‑135 (MQ=255)
                                                                                                tGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGctc                                    >  5:267043/1‑135 (MQ=255)
                                                                                                      cgaacgaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGCCTCGAACTCGCGACCCCGCCCTTGGCAAGGACGTGCACTACCAAc                             >  4:218190/1‑136 (MQ=255)
                                                                                                       caactaaAAAAGCCCCCTAGGGGCTTTTTTTAAACGACTTACCGACTTCGACTTAAACCTGTATATCTTTATTTGGCGCGGGAACCGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACt                            <  5:1152/135‑1 (MQ=18)
                                                                                                         actaaaaaaGGCCCCTAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTTAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTg                           <  1:208957/132‑1 (MQ=255)
                                                                                                           gaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAACCCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCACCTGAGc                        >  2:101286/1‑136 (MQ=255)
                                                                                                           gaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGc                        >  5:121431/1‑136 (MQ=255)
                                                                                                           gaaAAAAGGCCCCGAAGGGCCTATTTTAATCGACTTACAGACTTCGACTGAAACATGAAGATCTTTATATGGACCGGGAAACGAGACTCGAACTAGCGACCCCGACGTTGGCAAGGTCGTGCTCTACCAACTGAGc                        >  1:19822/1‑136 (MQ=21)
                                                                                                                          aGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGcgtcgt          <  7:94254/135‑1 (MQ=255)
                                                                                                                               cTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATggg    <  7:103222/136‑1 (MQ=255)
                                                                                                                                 tttttCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTg  >  4:89463/1‑136 (MQ=255)
                                                                                                                                 tttttAAACGACTTAAATACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATggg    <  7:214189/134‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                  
AGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTG  >  ACB143/4623977‑4624241

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