Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,348,419 T→G 22.1% D4A (GAC→GCC)  CJ019_01261 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,348,4190TG22.1% 65.8 / 3.2 41D4A (GAC→GCC) CJ019_01261hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (2/29);  new base G (7/3);  total (9/32)
Fisher's exact test for biased strand distribution p-value = 1.63e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTCGTCGGTGTAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATGCGCAATGCATGCGCATCCGCCC  >  ACB122/1348287‑1348540
                                                                                                                                    |                                                                                                                         
ttGTCGTCGGTGTACGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCTGGTGCGGTGGCAGCAAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCGGGGCCTGCGTGTTCTCGCGAGTCGTGTc                                                                                                                          <  5:199343/134‑1 (MQ=255)
 tGTCGTCGGTGTAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCgg                                                                                                                        <  5:145918/135‑1 (MQ=255)
   tcgtcgGTGTAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTc                                                                                                                          <  1:49697/131‑1 (MQ=255)
   tcgtcgGTGTAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCgg                                                                                                                        <  1:166586/133‑1 (MQ=255)
    cgtcgGTGTAAGGGACGACAGTGCCGGTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGtt                                                                                                                      <  1:164317/134‑1 (MQ=255)
          tgtAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGGGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCAtgg                                                                                                               <  5:100330/135‑1 (MQ=255)
            tAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCAtgggt                                                                                                             <  5:345470/135‑1 (MQ=255)
                    gacAGTGCCGTTTCGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGtct                                                                                                     <  1:259055/135‑1 (MQ=255)
                     acAGTGCCGTTACGCTTGATCACGCGGAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGtctc                                                                                                    <  6:194097/135‑1 (MQ=255)
                     acAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGtctc                                                                                                    <  5:69017/135‑1 (MQ=255)
                     acAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGtctc                                                                                                    <  1:109601/135‑1 (MQ=255)
                       aGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTGGGGTTGGAAAGGTGGGTGgctcgcc                                                                                                 >  8:345535/1‑133 (MQ=255)
                       aGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGGTTTGGTGGGTGGgggtggggc                                                                                                  >  1:273704/1‑130 (MQ=255)
                       aGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGGTTTGATGGGGGGgtgtgtgca                                                                                                  >  8:253136/1‑130 (MQ=255)
                       aGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGCGGGTTGGGTTGGTGGgtgtggaccc                                                                                                 >  7:194560/1‑130 (MQ=255)
                        gTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCcac                                                                                                 <  5:356856/135‑1 (MQ=255)
                              ttACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGCGGGTTGGATGGGTGGGTGTGTGCACAGtttg                                                                                           >  2:324599/1‑134 (MQ=255)
                                    ttGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATAtt                                                                                     <  6:197821/135‑1 (MQ=255)
                                       aTCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACGTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTc                                                                                 <  8:130210/136‑1 (MQ=255)
                                       aTCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTc                                                                                 <  8:91901/136‑1 (MQ=255)
                                        tCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGtt                                                                                  <  6:215621/134‑1 (MQ=255)
                                        tCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTc                                                                                 <  6:230565/135‑1 (MQ=255)
                                        tCACGCGCAGTTGGCCGGGGGCGGGGGCAGCCAGATCCTGGTTGGAAGCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTc                                                                                 <  6:239998/135‑1 (MQ=255)
                                         cACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGGTTGGATGGGTGGGTGTGTCCACAGTTTCTTTATttttact                                                                               >  4:152854/1‑133 (MQ=255)
                                          aCGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTGGGTTTGCCTTGGGGGGTGTGTCCCCAGTTTCTATATTTTTAAgg                                                                              >  3:162945/1‑136 (MQ=255)
                                                    ggCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGTTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATg                                                                     <  4:379487/135‑1 (MQ=255)
                                                    ggCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATg                                                                     <  4:60866/135‑1 (MQ=255)
                                                    ggCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATg                                                                     <  7:116146/135‑1 (MQ=255)
                                                    ggCCGGGGGCGGTGGCAGCCAGATCCTGGGTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATg                                                                     <  7:253136/135‑1 (MQ=255)
                                                         ggggCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGCGGGTTGGCTCGGTGGGGGGGTCCCCAGGTTCTATATTTGTGAGGCGGCTTATGGgcgc                                                                >  2:357758/1‑135 (MQ=255)
                                                                  gCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGGTTGGATGGGTGGGTGTGTGCCCAGGTTCTATATTTTTAAGGCGCCGTATGGGGGcccccccgtcc                                                      >  8:207143/1‑131 (MQ=255)
                                                                      ccAGATCCGGGTTGGAATCGGCGGCCGGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTGGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTGCAGGCGCCTTATGGGCGCCCACCGTTCCGTcc                                                  <  7:364358/136‑1 (MQ=255)
                                                                                     aaTCGGCGGCCTGCGGCACCTGGGCCGGCGGGTTCTCGCGAGTCGTGGCGGTTTGCATGGGTGGGTGTCGCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAAcc                                   <  4:301953/136‑1 (MQ=255)
                                                                                          gcggcCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGGGTCTCCACAGTTTCTATATTTTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTcgcg                              <  7:345535/136‑1 (MQ=255)
                                                                                            agcgaccGGCCCCTGGGCCGGCGGGTTCTCGCGCGTGGTGGCGGTTTGCATGGGGGGGTGTCTCCACAGTTTCTATATTGTTCCGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCgg                             <  1:344513/129‑1 (MQ=255)
                                                                                               agtcGGCACCGGGGCCTGCGGGTTATCGCGAGTCGTGTCGGTTTGCATGGGTGGGTTGCACAACAGTTTCTATATTGTTCCGGCGACTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGAt                         <  7:211323/133‑1 (MQ=255)
                                                                                                ggCGGCACGTGGGCCGGCGGGTTCTCGCGAGGCGGGTCGGTTTGCAGGGGTGGGGGGCGCCACAGTTTCTATATTGTTCCGGCGCCTTAGGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGAt                         <  5:29274/134‑1 (MQ=255)
                                                                                                 gCGGCACCTGGGCCGGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATgc                       <  4:13191/136‑1 (MQ=255)
                                                                                                       ccggggCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATGCGCAatg                 <  8:194560/133‑1 (MQ=255)
                                                                                                              cggCGGGTTCTCGCGAGTCGGGGCGGTTGGCCGGGGGGGGGGTCCCCACAGTTTCTATATTGTGCCGGCGCCTTAGGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATGCGCAATGCATgcg           <  5:377136/133‑1 (MQ=255)
                                                                                                                       ctcGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATGCGCAATGCATGCGCATccgccc  <  4:315826/135‑1 (MQ=255)
                                                                                                                                    |                                                                                                                         
TTGTCGTCGGTGTAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATGCGCAATGCATGCGCATCCGCCC  >  ACB122/1348287‑1348540

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