Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 5,644,213 T→C 9.1% D188G (GAC→GGC)  CJ019_05296 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,644,2130TC 9.1% 129.8 / 5.0 44D188G (GAC→GGC) CJ019_05296hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/15);  new base C (2/2);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

TCGGCCTTCATCGCCGCTATCGCCTGCTCGAAGCGCTCGACGATCCGATCATGCTCCGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCGCAGAAGAAAAGTTGCGCACCGGCACTTTGTGCAGACGGGGATCCCCGT  >  ACB143/5644078‑5644348
                                                                                                                                       |                                                                                                                                       
tCGGCCTTCATCGCCGCTATCGCCTGCTCGAAGCGCTCGACGATCCGATCATGCTCCGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGt                                                                                                                                         <  5:97626/136‑1 (MQ=255)
tCGGCCTTCATCGCCGCTATCGCCTGCTCGAAGCGCTCGACGATCCGATCATGCTCCGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGt                                                                                                                                         <  5:97629/136‑1 (MQ=255)
     cTTCATCGCCGCTATCGCCTGCTCGAAGCGCTCGACGATCCGATCATGCTCCGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGc                                                                                                                                    <  6:227487/136‑1 (MQ=255)
        cATCGCCGCTATCGCCTGCTCGAAGCGCTCGACGATCCGATCATGCTCCGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCg                                                                                                                                   <  5:212225/134‑1 (MQ=255)
          tCGCCGCTATCGCCTGCTCGAAGCGCTCGCCGATCCGATCATGCTCCGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTc                                                                                                                               >  1:277580/1‑136 (MQ=255)
            gccgcTATCGCCTGCTCGAAGCGCTCGACGATCCGATCATGCTCCGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCCCGCACGTCCTGCGGCTCAc                                                                                                                             <  2:210418/136‑1 (MQ=255)
             ccgcTATCGCCTGCTCGAAGCGCTCGACGATCCGATCATGCTCCGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACg                                                                                                                            <  5:199163/136‑1 (MQ=255)
             ccgcTATCGCCTGCTCGAAGCGCTCGACGATCCGATCATGCTCCGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACg                                                                                                                            <  3:316784/136‑1 (MQ=255)
                   tCGCCTGCTCGAAGCGCTCGACGATCCGATCATGCTCCGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACAAGCATCACCCCATCACGCACGTCATGCGGATCACGCTGCAg                                                                                                                      >  2:232645/1‑136 (MQ=255)
                             gAAGCGCTCGACGATCCGATCATGCTCCGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGc                                                                                                              <  2:286606/134‑1 (MQ=255)
                             gAAGCGCTCGACGATCCGATCATGCTCCGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCgg                                                                                                            <  6:94573/136‑1 (MQ=255)
                                   cTCGACGATCCGATCATGCTCCGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACg                                                                                                      <  2:318034/136‑1 (MQ=255)
                                        cgatccgatcATGCTCCGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTc                                                                                                 <  6:295306/136‑1 (MQ=255)
                                           tccgatcATGCTCCGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCTGGCCACGTACTCGTc                                                                                              <  8:19200/136‑1 (MQ=255)
                                                tcATGCTCCGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCa                                                                                         <  5:119043/136‑1 (MQ=255)
                                                     cTCCGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGCCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGcc                                                                                    <  3:207548/136‑1 (MQ=255)
                                                        cGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTc                                                                                                                                        <  2:213996/81‑1 (MQ=255)
                                                        cGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTc                                                                                                                                        >  1:213996/1‑81 (MQ=255)
                                                                 cAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGc                                                                        >  3:65432/1‑136 (MQ=255)
                                                                  aGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGcc                                                                       >  4:123339/1‑136 (MQ=255)
                                                                   ggCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGccc                                                                      >  1:208481/1‑136 (MQ=255)
                                                                           aGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTTGCGGGCGACGTACTCGTCTTCCACCGCCAGGCTTACCAGCCCTgccgcc                                                               >  2:225168/1‑135 (MQ=255)
                                                                           aGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTgccgcc                                                               >  4:74283/1‑135 (MQ=255)
                                                                                atGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATg                                                         >  3:271261/1‑136 (MQ=255)
                                                                                     gggTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTAc                                                    <  5:281534/136‑1 (MQ=255)
                                                                                     gggTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGCCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTAc                                                    >  3:295992/1‑136 (MQ=255)
                                                                                         gtTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCg                                                  >  5:51613/1‑134 (MQ=255)
                                                                                          ttCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCg                                                  >  6:277031/1‑133 (MQ=255)
                                                                                          ttCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCg                                                  >  8:89782/1‑133 (MQ=255)
                                                                                          ttCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGCCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCg                                                  >  6:48510/1‑133 (MQ=255)
                                                                                              cccccAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCGCagaag                                            >  4:303703/1‑135 (MQ=255)
                                                                                                   aGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCGCAGAAGAAAAGt                                      >  2:98208/1‑136 (MQ=255)
                                                                                                     cggcggATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCGCAGAAGAAAAGTTg                                    >  4:73623/1‑136 (MQ=255)
                                                                                                     cggcggATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCGCAGAAGAAAAGTTg                                    >  6:120579/1‑136 (MQ=255)
                                                                                                       gcggATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCGCAGAAGAAAAGTTgcg                                  >  8:299586/1‑136 (MQ=255)
                                                                                                              caccaGCATCACCCCCTCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCGCAGAAGAAAAGTTGCGCACCGGc                           <  5:240436/136‑1 (MQ=255)
                                                                                                                   gCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCGCAGAAGAAAAGTTGCGCACCGGCACttt                      <  7:153050/136‑1 (MQ=255)
                                                                                                                           ccctcactcacGCCCTGCGGCTCACGCTGCAGGGTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCGCAGAAGAAAAGTTGCGCACCGGCACTTTGTGCAGAc              <  6:209749/133‑1 (MQ=255)
                                                                                                                           ccATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCGCAGAAGAAAAGTTGCGCACCGGCACTTTGTGCAGAc              >  1:84220/1‑136 (MQ=255)
                                                                                                                            cATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCGCAGAAGAAAAGTTGCGCACCGGCACTTTGTGCAGACg             >  1:163893/1‑136 (MQ=255)
                                                                                                                               cacgcacgTCCTGCGGCTCACGCTGCAGGTAGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACAAGCGTACCGCAGAAGAAAAGTTGCGCACCCGCACTTAGTGCAGACgggg          >  4:13816/1‑136 (MQ=255)
                                                                                                                                 cgcacgTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCTCCCTGCCGCCAACATGCGTACCGCAGAAGAAAAGTTGCGCACCGGCACTTTGTGCAGACGGGGAt        >  2:232019/1‑136 (MQ=255)
                                                                                                                                 cgcacgTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCGCAGAAGAAAAGTTGCGCACCGGCACTTTGTGCAGACGGGGAt        >  2:539/1‑136 (MQ=255)
                                                                                                                                 cgcacgTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACAGCAGAAGAAAAGTTGCGCACCGGCACTTTGTGCAGACGGGGAt        >  8:196908/1‑136 (MQ=255)
                                                                                                                                    acgTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCGCAGAAGAAAAGTTGCGCACCGGCACTTTGTGCAGACGGGGATccc     >  2:6273/1‑136 (MQ=255)
                                                                                                                                      gTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCGCAGAAGAAAAGTTGCGCACCGGCACTTTGTGCAGACGGGGATCCCCg   >  6:173936/1‑136 (MQ=255)
                                                                                                                                       tCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTCCCGCAGAAGAAAAGTTGCGCACCGGCACTTTGTGCAGACGGGGATCCCCGt  >  3:72509/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                                       
TCGGCCTTCATCGCCGCTATCGCCTGCTCGAAGCGCTCGACGATCCGATCATGCTCCGGGTGCTTCAGGCTGACCAGAATATGCAGGGTGTTCTCCCCCAGCGGCGGATCCACCAGCATCACCCCATCACGCACGTCCTGCGGCTCACGCTGCAGGTTGTAGCGGGCCACGTACTCGTCTTCCACCGCCAGGCTTACCCGCCCTGCCGCCAACATGCGTACCGCAGAAGAAAAGTTGCGCACCGGCACTTTGTGCAGACGGGGATCCCCGT  >  ACB143/5644078‑5644348

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