Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 2,737,398 G→C 33.8% Q241H (CAG→CAC)  CJ019_02580 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,737,3980GC33.8% 26.4 / 3.5 26Q241H (CAG→CAC) CJ019_02580hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (14/3);  new base C (3/6);  total (17/9)
Fisher's exact test for biased strand distribution p-value = 2.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTACAGCCGACGGTGCTTGAACGTGATTTCAACTATCCGCCC‑TTGGCCGAACTGTTGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGT  >  ACB143/2737272‑2737519
                                                                                                                              |                                                                                                                          
ggCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTACAGCCGACGGTg                                                                                                                   >  4:132819/1‑136 (MQ=255)
                   gggCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTACAGCCGACGGTGCTTGAACGTGATTTCCAc                                                                                                 >  8:63620/1‑135 (MQ=255)
                            gATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTTTACGCCAGCCCTTGTGCGCGCCTTGGCGTACCCCCCCCCGGGCTTGAACGTGATTTGAACAAAACGCCC‑t                                                                                      >  1:311312/1‑136 (MQ=255)
                                      cGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTACAGCCGACGGTGCTTGAACGTGATTTCAACTATCCGCCC‑TTGGCCGAACt                                                                            >  2:275982/1‑136 (MQ=255)
                                        gACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTAAGCCACCGCTTGTGCGCGCTTTGGCGTTCACCCGACGGTGGTTGAACGTGATTTAAACTAACCCCCC‑TTGGCCGCACGGt                                                                          >  7:90119/1‑136 (MQ=255)
                                           cTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTACAGCCGACGGTGCTTGAACGTGATTTCAACTATCCGCCC‑TTGGCCGAACTGTTgg                                                                       >  5:372520/1‑136 (MQ=255)
                                                    taTACCCATGGGGCAGCCGTGAAGGAGGAGGGCTGGGGGCTGTACGCCAGCGCTTGTGCGCGCTTGGGCGTACAGCCGACGGTGCTTGAACGTGATTTCAACTATCCGCCC‑TTGGCCGAACTGTTGGCCGAGACg                                                               <  8:212683/134‑1 (MQ=255)
                                                        cccATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACCCCCCCCCGTGTGCCCGCCTTGGCGTACAGCCCCCGGTGCGTGAACCTGATTTCAACAATCCCCCC‑TTGGCCGCACTGTTGGCCGCGACcgcgcg                                                          >  3:46384/1‑136 (MQ=255)
                                                             ggggCAGCCTTGACGGCGGTTGTGGGGGCGCTGACCGCCAGCGCTTGGGCGCGCTTTGGCTAACACCCGACGGTGCTTGAACGTGATTTCAACTATCCGCCC‑TTGGCCGAACTGTTGGCCGAGACGGCGCGTATGc                                                     <  5:186491/136‑1 (MQ=255)
                                                              gggCAGCCGTGAAGGAGGATGTCTGGGCGCTGTTAGCCACCGCTTTTGGCCGCGTTTGCGGTCACCCGACCGGGCGTTAACGTGATTTAAACTATCCGCCCATTGGCC‑CACTGTTTGCCGACACCGCGCCTATTcg                                                    >  6:20175/1‑136 (MQ=255)
                                                               ggCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTACAGCCGACGGTGCTTGTACGTGATTTCAACTATACGCCC‑TTGGCCGAACTGTTGGCCGAGACGGCGCGTATGcgt                                                   >  3:169115/1‑136 (MQ=255)
                                                                 cAGCCGTGAAGGAGGATGTCTGGGCGCTTTACGCCCTCGCTTGTGAGCGCGTTGTCGTACAGCCGCCGGTTCATGAACATGATTTCCACTATCCACCC‑TTTTCCTAAATGGTGGCCGACACGGCGCGGATTcgtgc                                                 >  1:55029/1‑136 (MQ=255)
                                                                      gTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTACAGCCGACGGTGCTTGAACGTGATTTCAACTATCCGCCC‑TTTGCCGCACTGTTGGCCGAGACGGCGCGTATGcgtgccgtg                                             >  8:178390/1‑135 (MQ=255)
                                                                      gTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTACAGCCGACGGTGCTTGAACGTGATTTCAACTATCCGCCC‑TTGGCCGAACTGTTGGCCGAGACGGCGCGTATGcgtgccgtg                                             >  1:112412/1‑135 (MQ=255)
                                                                            gaggaTGTCTGGGCGCTGTACGCCAGCCCTTGTGCGCGCATTGGCGTACAGCCGCCGGTGCTTGAACGTGATTTAAACTATCCGCCC‑TTGGCCCAACTTTTGGCCGAGACCGCGCGTATGCGTGCCGTGCAATtcg                                      >  8:220519/1‑136 (MQ=255)
                                                                            gaggaTGTCTGGGCGCTGTAAGCCAGCGCTTGTTCGCGCTTTGGCGTACAGCCGACGGTGGGTGAACGTGATTTGAACTATCCGCCC‑GTGGCCGAACGGTTGGCCGAGACGGCGCGTATGCGTGCCGGGCAATgcg                                      >  4:255234/1‑136 (MQ=255)
                                                                                     tGGGCGCTGTACCACAGCGATTGGGCGCGTTTTGGATTACCCCCGACGGTGCTTGAACGTGATTTCAACTATCCGCCC‑TTGGCCGAACTGTTGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAg                             <  6:118690/136‑1 (MQ=255)
                                                                                          gcTGTACGCCACCCCCTGTGCGCGCGTTGGCGTTCAGCCGCCCGTGGATGCACCTGATTTGAACTATCCCCCC‑TTGGCCCCACAGTTGGCCGAGACGGCGCGGAATCGTGCCGTGCAATGCCCGGCCGGAGGGcca                        >  6:140385/1‑134 (MQ=255)
                                                                                           cTGTCCGCCTGGGCTGGTGGGCGGTTTGGGCTACAGCCGACGGTGCTTGAACGTGATTTCAACTATCCGCCC‑TTGGCCGAACTGTTGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAg                       <  1:260665/136‑1 (MQ=255)
                                                                                             gtccGCCGGCGGTGGGGCGGGCTTGGGCTAACCCCCGACGGTGCTTGAACGTGATTTCAACTATCCGCCC‑TTGGCCGAACTGTTGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCa                     <  3:404300/133‑1 (MQ=255)
                                                                                               aCGCCACCGCTTGTGCGCGCTTTGGCGTACAGCCGACGGGGCTTGAACCTGATTTCAACTATACGCCC‑TTGGCCGAACTGTTGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCAGCCGGAGGGCCAg                       >  2:148384/1‑132 (MQ=255)
                                                                                                         ttGTGCGCGCTTTGGCGTACAGCCGACGGTGCTTGAACGTGATTTCAACTATCCGCCC‑TTGGCCGAACTGTTGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGAc         >  5:352208/1‑136 (MQ=255)
                                                                                                                ggttttGGGCTACACCCGACGGTGCTTGAACGTGATTTCAACTATCCGCCC‑TTGGCCGAACTGTTGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCtgcgt  <  7:85571/133‑1 (MQ=255)
                                                                                                                ctttttGGCTTACCCCCGACGGTGCTTGAACGTGATTTCAACTATCCGCCC‑TTGGCCGAACTGTTGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCtgcgt  <  3:258571/133‑1 (MQ=255)
                                                                                                                cgtttgggCGAACAGCCGACGGTGCTTGAACGTGATTTCAACTATCCGCCC‑TTGGCCGAACTGTTGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCtgcgt  <  3:132819/130‑1 (MQ=255)
                                                                                                                cggtttGGCGTACCCCCGACGGTGCTTGAACGTGATTTCAACTATCCGCCC‑TTGGCCGAACTGTTGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCtgcgt  <  2:79537/133‑1 (MQ=255)
                                                                                                                              |                                                                                                                          
GGCCGGGATGCACGTTGCCGGGCATTACGATGAAGCGCCGGACCTCAAAGTGGATACCCATGGGGCAGCCGTGAAGGAGGATGTCTGGGCGCTGTACGCCAGCGCTTGTGCGCGCTTTGGCGTACAGCCGACGGTGCTTGAACGTGATTTCAACTATCCGCCC‑TTGGCCGAACTGTTGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGT  >  ACB143/2737272‑2737519

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: ‑