Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,490,632 G→A 8.3% Q3* (CAG→TAG)  CJ019_03306 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,490,6320GA 8.3% 149.5 / 2.7 49Q3* (CAG→TAG) CJ019_03306hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/19);  new base A (2/2);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGCCAAACAGGGTTTGCGGACGGATCAGCGGTGCCACGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCGGCTGCAGCGAATTTGAGTGCTGGCATATGAATATTGTGTTGCCAGCAATGCCGA  >  ACB122/3490498‑3490764
                                                                                                                                      |                                                                                                                                    
ttCGCCAAACAGGGTTTGCGGACGGATCAGCGGTGCCACGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCtg                                                                                                                                      <  4:367272/135‑1 (MQ=255)
   gCCAAACAGGGTTTGCGGACGGATCAGCGGTGCCACGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTgtgt                                                                                                                                   >  4:108527/1‑135 (MQ=255)
   gCCAAACAGGGTTTGCGGACGGATCAGCGGTGCCACGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTgtgt                                                                                                                                   >  3:254447/1‑135 (MQ=255)
     cAAACAGGGTTTGCGGACGGATCAGCGGTGCCACGGTGGCAAGGTAATCCGCCGCCTGCTCCGTCAGCGCCTAGCACAACGCCCCCGATTGGCTACGCAGCTTGAAGAGTGCTTCGACCTCGAATCACTGTGTGAt                                                                                                                                >  8:141552/1‑136 (MQ=255)
          aGGGTTTGCGGAGGGATCAGCGGTGCCACGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACTTCGAATGTCTGTGTCATTGGt                                                                                                                            <  2:144080/135‑1 (MQ=255)
          aGGGTTTGCGGACGGATCAGCGGTGCCACGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGt                                                                                                                            <  2:259553/135‑1 (MQ=255)
           gggTTTGCGGACGGATCAGCGGTGCCACGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGCCATCGAATGTCTGTGTCATTGGt                                                                                                                            >  2:51706/1‑134 (MQ=255)
           gggTTTGCGGACGGATCAGCGGTGCCACGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGc                                                                                                                          <  7:141552/136‑1 (MQ=255)
                     acggaTCAGCGGTGCCACGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCAt                                                                                                                 <  2:92875/135‑1 (MQ=255)
                     acggaTCAGCGGTGCCACGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCAt                                                                                                                 <  2:176786/135‑1 (MQ=255)
                         aTCAGCGGTGCCACGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGcc                                                                                                            >  4:172866/1‑136 (MQ=255)
                                gTGCCACGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCGCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGa                                                                                                     <  8:137/136‑1 (MQ=255)
                                     aCGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTc                                                                                                >  7:254607/1‑136 (MQ=255)
                                      cGGTGGCAAGGTAATCCGCCGCCTGCGGCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCAGTGTCATTGGTGCATCGTCCATGTGGCTGCTTGATCGTcc                                                                                               >  7:235312/1‑136 (MQ=255)
                                      cGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTTGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCCTGGTCCATGTGCCTGCTTGATCGTcc                                                                                               >  3:69620/1‑136 (MQ=255)
                                      cGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTc                                                                                                                                  <  3:22200/101‑1 (MQ=255)
                                      cGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTc                                                                                                                                  >  4:22200/1‑101 (MQ=255)
                                      cGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTcc                                                                                               >  7:241391/1‑136 (MQ=255)
                                      cGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGGGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTcc                                                                                               >  7:23086/1‑136 (MQ=255)
                                      cGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATAGGAGATAGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGAACGTcc                                                                                               >  7:47632/1‑136 (MQ=255)
                                       ggTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTTTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTcc                                                                                               <  2:243113/135‑1 (MQ=255)
                                              aGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACGc                                                                                       >  3:301512/1‑136 (MQ=255)
                                                      cgccgcCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACGCCGGTGGCg                                                                               >  3:189089/1‑136 (MQ=255)
                                                      cgccgcCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCGGCTGGATCGTCCTGAAACGCCGGTGGCg                                                                               >  7:372995/1‑136 (MQ=255)
                                                      cgccgcCTGCGCCTTAAGCGCAAAGGAGATAGCCCGCGTTTTGATACGCAGCGTTAACAGTTCTTAGTAAGAGAATGTGTGTGTGCTTTGTTCTTAGTACATGTTGCTGCTTGATCGTCCTGAAACGCCGGTgggg                                                                               >  8:163529/1‑134 (MQ=255)
                                                        ccgccTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGGATCGTCCATTTGCCTGCTTGATAGTCCTGAAACGCCGGTGGCGt                                                                              >  6:325130/1‑135 (MQ=255)
                                                         cgccTGCGCCTTCAGCGCATCGGAGAGCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGGCATTGGGGCAGCGTCCATGTGCCGGCTTGATCGGCCTGAAACGCCGGTGGCGTgg                                                                            <  4:269436/136‑1 (MQ=255)
                                                            cTGCGCCTTCAGCGCATAGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATAGTCCATGTGCCTGCATGATCGTCCTGAAACGCCGGTGGCGTGGAc                                                                          >  2:299830/1‑135 (MQ=255)
                                                                 ccTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGAGCGGCCTGAAACGCCGGTGGCGTGGACTTGaaa                                                                    <  7:33610/136‑1 (MQ=255)
                                                                  cTTCCGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACGCCGGTGGCGTGGACTTGaa                                                                     >  5:196006/1‑134 (MQ=255)
                                                                              ggAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTATGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCgg                                                       >  4:236761/1‑136 (MQ=255)
                                                                              ggAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTATGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCgg                                                       >  4:241510/1‑136 (MQ=255)
                                                                               gagaTCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCgg                                                       >  2:38073/1‑135 (MQ=255)
                                                                               gagaTCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCGGc                                                      >  3:262881/1‑136 (MQ=255)
                                                                               gagaTCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCGGc                                                      >  4:311924/1‑136 (MQ=255)
                                                                                   tCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTCGATCGTCCTGAAACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCGGCTg                                                    >  7:93179/1‑134 (MQ=255)
                                                                                             ttGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGt                                                                                                                            >  5:149519/1‑52 (MQ=255)
                                                                                             ttGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGt                                                                                                                            <  6:149519/52‑1 (MQ=255)
                                                                                                gcgaCGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATGGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGACACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCGGCTGCAGCGAATTTGAGt                                      <  1:142654/132‑1 (MQ=255)
                                                                                                 cTACGCTGCTTGATCCGTGCTTCGACATCTAATGTCTGTGTAATTGGTGCATCGTCCATGTGCCTGCTTGATCGGCCTTAAACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCGGCTGCAGCGAATTTGAGTGc                                    <  4:262881/136‑1 (MQ=255)
                                                                                                  tACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCGGCTGCAGCGAATTTGAGTGc                                    <  5:288384/135‑1 (MQ=255)
                                                                                                    cGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCGGCTGCAGCGAATTTGAGTGCTg                                  <  7:301650/135‑1 (MQ=255)
                                                                                                      ctgCTTGATCAGGGCTTCTACATCGAATGTGTTTGTCATTGGGGCATCTGCCATGTGCCTGCTTGAGCGCCCTGAAACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCGGCTGCAGCGAATTTGAGTGCTGGCa                               <  3:8790/134‑1 (MQ=255)
                                                                                                               cAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCGGCTGCAGCGAATTTGAGTGCTGGCATATGAatat                      >  8:204134/1‑136 (MQ=255)
                                                                                                                  tGCTTCGACATCGAATGTCTATGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCGGCTGCAGCGAATTTGAGTGCTGGCATATGAATATtgt                   <  3:236761/136‑1 (MQ=255)
                                                                                                                  tGCTTCGACATCGAATGTCTATGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCGGCTGCAGCGAATTTGAGTGCTGGCATATGAATATtgt                   <  3:241510/136‑1 (MQ=255)
                                                                                                                   ccTTCGACAGCGAATGTCTGGGTCATTGGTGCATCGCCCATGGCCCTGCTTGCTCGTCCTCACACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCGGCTGCAGCGAATTTGAGTGCTGGCATATGAATATtgt                   <  5:100372/134‑1 (MQ=255)
                                                                                                                    cTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACg                                                                                        >  7:386682/1‑65 (MQ=255)
                                                                                                                    cTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACg                                                                                        <  8:386682/65‑1 (MQ=255)
                                                                                                                    cTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCGGCTGCAGCGAATTTGAGTGCTGGCATATGAATATTgtgt                 <  8:254607/136‑1 (MQ=255)
                                                                                                                      tCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACGCCGGTGGCGTGGACTTGAAAAATCGCATTTCGGCTGCAGCGAATTTGAGTGCTGGCATATGAATATTGTGtt                >  5:158023/1‑135 (MQ=255)
                                                                                                                                   tCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCGGCTGCAGCGAATTTGAGTGCTGGCATATGAATATTGTGTTGCCAGCAATGCCGa  <  4:41591/136‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                    
TTCGCCAAACAGGGTTTGCGGACGGATCAGCGGTGCCACGGTGGCAAGGTAATCCGCCGCCTGCGCCTTCAGCGCATCGGAGATCGCCCGCGTTTGGCTACGCAGCTTGATCAGTGCTTCGACATCGAATGTCTGTGTCATTGGTGCATCGTCCATGTGCCTGCTTGATCGTCCTGAAACGCCGGTGGCGTGGACTTGAAAAATCGCCTTTCGGCTGCAGCGAATTTGAGTGCTGGCATATGAATATTGTGTTGCCAGCAATGCCGA  >  ACB122/3490498‑3490764

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: