Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 740,674 T→G 50.5% L38R (CTG→CGG)  CJ019_00704 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122740,6740TG50.5% 1.5 / 10.3 21L38R (CTG→CGG) CJ019_00704hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (7/2);  new base G (8/3);  total (16/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

GGATCTCCGTATGAGAACGGTTTTGTTGTGTCTTGGGCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTGGGGGGGGCGCTGTGGGTTCTGGCTGAATGGCTGGAGCGGAGGGGGCGCAATCAATGATGCGCGCCAATCTAGACGAGCTGGTGATCGTGAGCTGTTTGTTTCCAGCTATGAGGTGGTCAGCTTCTATTGCGCGATCCGT  >  ACB122/740552‑740802
                                                                                                                          |                                                                                                                                
ggaTCTCCGTATGAGAACGGTTTTGTTGTGTCTTGGGCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGGTGTGGGTGGTGGGGGGGGGGCGGGGGGTTGTggg                                                                                                                      >  7:149481/1‑134 (MQ=255)
            gagaACGGTTTTGTTGTGTCTTGGGCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGGTGTGGGTGGTGGGGGGGGCGCGGTTGGGTGTGGGTGAAAGGGTGGa                                                                                                          >  6:294339/1‑135 (MQ=255)
                  ggTTTTGTTGTGTCTTGGGCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTTGTGGGGGGGGCGCTGTGGGGTCTGGGGGAATGGCGGGAGCGGAg                                                                                                    >  1:43962/1‑135 (MQ=255)
                       tgttgtGTCTTGGGCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCAGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTGGGGGGGGAGCTGTGGGGTGTTGATGAATGGCGGGAGCGGAGGGGgc                                                                                               >  4:84116/1‑135 (MQ=255)
                       ggttgtGGCTGGTGCTTGCACTGGCTGGGTGGTGTGTGGGCTCTCATGGGCCCACACCACTCTGGGGCTTCGGGGTGGTGTGGGTGGGGGGGGGGGCGCGGTGGGTGCTGGCTGAATGGCTGGAGCGGAggggg                                                                                                <  6:237774/133‑1 (MQ=255)
                                 tGGGCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGGTGTGGGTGGTGGGGGGGGGGGGGTGGGTGGTTGCGGATTGGGGGGTTGGGCGGGGGCGCTATCCatg                                                                                     >  1:186746/1‑135 (MQ=255)
                                 tGGGCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGGTGTGGGGGGTGGGGGGGGGGCGGGGGGTTGGGGGGGGATGGGGGGGGCGGCGGGGGGGccatcaatt                                                                                     >  5:45815/1‑134 (MQ=255)
                                 tGGCCTTGCAAGTGCACGCTTCTGTCTCGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGAGGGAGGTGGGGGGGGCGCGGTGGGTTGATGCTGAATGAGTGGACCGGATGGGGCGGAATCAatg                                                                                     >  4:11556/1‑135 (MQ=255)
                                                     tCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTGGGGGGGGCGCAGTGGGTTCGGGCTGAATGGCGGGGGCGGAGGGGGCGCAATCAATGATGCGCGCCCATCTAGACGa                                                                 >  3:129242/1‑135 (MQ=255)
                                                                tctcACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTGGGGGGGGCGCTGTGGGt                                                                                                                            <  6:150875/65‑1 (MQ=255)
                                                                tctcACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTGGGGGGGGCGCGGTGGGt                                                                                                                            >  5:150875/1‑65 (MQ=255)
                                                                      tGGCCGACACCACTCTGGGTCATCGTCGTGGTGTGGGTGGTGGGGGGGGCGCGGTGGGTTGTGGGGGATTGGGGGGCGCCGGGGGGGCGGCATCCATGATGGGCGCCCATAGAGACGAGCTGgtggtggtggggt                                                >  1:230372/1‑131 (MQ=255)
                                                                           gACACCACTCTGGGTCATCGTCGTGGTGGGGGTGGTGGGGGGGGGGCGGGGGGGTGTGGCTGGATTGGTGGTGGGGGGGGGGGGGCATCACTGATGCGCGCCCATCTTGACGCGCTGGTGAGCGTGAGCggtgtg                                           >  6:5842/1‑135 (MQ=255)
                                                                                     tGGGTCATCGTCGTGTTGTGGGTGGTGGGGGGGGCGCTGTTGGTTATGGCTGCATGGCTGGAGCCGCCGGGGCCCAATCAATTATGCGCGCCCAAAAATAAGAGCTGGTGTTATTTAGCTGTTTTTTTCCAGCTa                                 >  2:58434/1‑135 (MQ=255)
                                                                                               tcgtGTTGTGGGTGGTGGGGGGGGCGCTGTGGGGTCTGGCTGAATGGCTGGCGCGGAGGGGGCGCAATCAATGATTCGCGCCAATCTAGACGAGCGGGTGATCGTGAGCTGTTTGTTTACAGCTATGAggtgg                         >  6:59086/1‑133 (MQ=255)
                                                                                                cgtGTTGTGGGTGGTGGGGGGGGCGCTGTGGGTTCTGGCTGAATGGCTGGAGCGGCGGGGGCGCAAACAATTATGCGCGCCAATATAGAAGAGCTGGTGTTCGTGAGCTGTTTTTTTCCAGCTATGAGGTGGTCa                      >  8:11738/1‑135 (MQ=255)
                                                                                                            ggtggGGGGGGCGCGGGGGGGGCTGGCTGAAGGGCTGGAGCGGCGGGGGCGCAATCAATGATGCGCGCCAATCTAGACGAGCTGGTGATCGTGAGCTGTTTGTTTCCAGCTATGAGGTGGTCAGCTTCTATTgcg          <  8:149481/135‑1 (MQ=255)
                                                                                                               ggggggggCGCTGTGGGTTCTGGCTGAATGGCTTGAGAGGAGGGGGCGCAATCAATGATGCGCGCCAATCTAGACGAGCTGGTGATCGTGAGCTGTTTGTTTCCAGCTATGAGGTGGTCAGCTTCTATTg            >  8:86704/1‑130 (MQ=255)
                                                                                                               ggggggggCGCTGTGGGTTCTGGCTGAATGGCTGGAGCGGAGGGGGCGCAATCAATGATGCGCGCCAATCTAGACGAGCTGGTGATCGTGAGCTGTTTGTTTCCAGCTATGAGGTGGTCAGCTTCTATTGCGCGa       >  2:288649/1‑135 (MQ=255)
                                                                                                                    gggCGGTGTGGGTGCGGGCTCAACGGCTGGAGCGCAGGGGGCGCAATCAATGATGCGCGCCAATCTAGACGAGCTGGTGATCGTGAGCTGTTTGTTTCCAGCTATGAGGTGGTCAGCTTCTATTGCGCGATCCGt  <  3:11556/135‑1 (MQ=255)
                                                                                                                    cggCGCGGGGGGTGCTGGCTCACTGGCGGGAGCGGCGGGGGCGCAATCAATGATGCGCGCCAATCTAGACGAGCTGGTGATCGTGAGCTGTTTGTTTCCAGCTATGAGGTGGTCAGCTTCTATTGCGCGATCCGt  <  7:11738/134‑1 (MQ=255)
                                                                                                                          |                                                                                                                                
GGATCTCCGTATGAGAACGGTTTTGTTGTGTCTTGGGCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTGGGGGGGGCGCTGTGGGTTCTGGCTGAATGGCTGGAGCGGAGGGGGCGCAATCAATGATGCGCGCCAATCTAGACGAGCTGGTGATCGTGAGCTGTTTGTTTCCAGCTATGAGGTGGTCAGCTTCTATTGCGCGATCCGT  >  ACB122/740552‑740802

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