Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,670,696 A→C 22.0% intergenic (+33/‑98) thrC_2 → / → CJ019_01569 Threonine synthase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,670,6960AC22.0% 79.2 / 2.2 54intergenic (+33/‑98)thrC_2/CJ019_01569Threonine synthase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/16);  new base C (12/2);  total (36/18)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCACCAGGCTCGC  >  ACB122/1670564‑1670829
                                                                                                                                    |                                                                                                                                     
ccTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCACCCCCAACg                                                                                                                                    <  7:328611/136‑1 (MQ=255)
            tGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCtttt                                                                                                                        >  8:23472/1‑136 (MQ=255)
               agagCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCCGCCCCCCCGTCGCCGCCTTTTcc                                                                                                                      >  4:393831/1‑135 (MQ=255)
               agagCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCAACGTCTCAGCCTTTTcc                                                                                                                      >  2:104228/1‑135 (MQ=255)
                 agCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTTGGGCTGCAGCGCAGCCCCCACGTATAAGCCTTTTCcctc                                                                                                                   >  7:35510/1‑136 (MQ=255)
                      tGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGt                                                                                                               >  5:160185/1‑135 (MQ=255)
                                   ccAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCAACCTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATc                                                                                                 >  7:256689/1‑136 (MQ=255)
                                     aaCGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCGTTTCCCTCCTGTGTATCAGAAACATc                                                                                                 >  6:70899/1‑134 (MQ=255)
                                          ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCCGCCCCCACGGCTCAGCCTTTTCCCTCCGGTCTATCAGAAACATCCGATTc                                                                                           >  4:379249/1‑135 (MQ=255)
                                          ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCCACGGCGACGCCGTTTACCTCCGGTGTTTAAGGAACATACTAttg                                                                                           >  7:161315/1‑134 (MQ=255)
                                          ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCCACGTCTCCGCCGTTTCCCTCCTGTCAATCATAAACATCCTAtta                                                                                           >  6:12800/1‑134 (MQ=255)
                                          ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAt                                                                                             >  6:85463/1‑133 (MQ=255)
                                          ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTc                                                                                           >  3:295783/1‑135 (MQ=255)
                                          ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCa                                                                                          >  7:111667/1‑136 (MQ=255)
                                          ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCGTTTCCCTCCTGTCTATCAGAAACATCCTATTc                                                                                           >  1:118200/1‑135 (MQ=255)
                                          ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCGGGCGCGCCGCCCCCACGGCGCCGCCTTTTACCACCTGTCGATAAGAAACAAACTAtt                                                                                            >  2:372660/1‑134 (MQ=255)
                                          ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGAGTTTGGGGCAGCAGCGCAGCCCCAACCTCTCAACCTTTTCCCTCCTGTCTATAAGAAACAAACTATTCa                                                                                          >  7:6346/1‑136 (MQ=255)
                                          ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGAAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCGTTTACCTCCTGTCTATCATAAACATCCTATTc                                                                                           >  6:257645/1‑135 (MQ=255)
                                          catGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCTCGGCAATCCGCTGTAATTGCATGGCGTTTGTGGCTGGCGCGGAGCCCCAAACTCTCCGGCATATCCCTCCGTTATAACAGAAACATACTAGTc                                                                                           >  6:112097/3‑135 (MQ=255)
                                                      accgggCTTTGCCCGCCAGCACGGCAACCGCGGAACGCCGCGGTATTTGCATGGCGTTGGGGGCTGCAGCGCAGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGgc                                                                               <  2:118200/132‑1 (MQ=255)
                                                      aCAGGGCTTTGCCCCCCCCCCCGGAAACCGCGGCACGCCGCTGAAATTGCAGGGCGTTGGGGGCGGCCGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGgc                                                                               <  7:261981/135‑1 (MQ=255)
                                                        aGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGGTGGCGCGCCGCCCCCACGGCGGAGCCGGTTCCCGCCCGTCGTTCAGCAACATCCGATTCAGGCTGCCGGgccc                                                                             >  6:42899/1‑133 (MQ=255)
                                                        aGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCCGCCCCAACGTCTCAGCCGTTTCCCTCCTGTCTATCAGCAACATCCTATTCAGCCTGCCGGgcgc                                                                             >  1:263168/1‑135 (MQ=255)
                                                        aGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGGCTCAGCCTTTTCCCTCCTGTCTATCAGAAAAATACTATTCAGCCTGGCGGgcgc                                                                             >  3:374090/1‑135 (MQ=255)
                                                        aGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCCCAGCCCCCACCGCTACCCATTTTCCCTACTTTCTATCAGAAACATCCTAATCAGCCTGCCGGgccc                                                                             >  4:217994/1‑133 (MQ=255)
                                                            cTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCCGCCCCCACGTCGCAGCCCTTTCCCGCCGGTGGATAAGAAACATCAGATTCCGCCGGCCGGGCGCCAAg                                                                         >  5:71671/1‑135 (MQ=255)
                                                             tttGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTGTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGc                                                                        >  6:291883/1‑135 (MQ=255)
                                                                gTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATa                                                                     >  8:387990/1‑135 (MQ=255)
                                                                gTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGGCGCCGCCTTTTCCCTCCTGGCTATCAGAAACATACTATTCAGCCTGCCGGGCGCCCAGCCTa                                                                     >  8:43574/1‑135 (MQ=255)
                                                                  cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGGGGCCGCGCCCCCCCCACCGCGCAGCCGTTTCCCGGCAGGCGATGAGTAACATCCGATTCAGCCAGCCCGGGGCCCAGCagta                                                                    >  5:106214/1‑131 (MQ=255)
                                                                  cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCAACGTCTCAGCCTTTTCCCTCCTGTATATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATaaa                                                                   >  3:119566/1‑135 (MQ=255)
                                                                  cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATaaa                                                                   >  2:101036/1‑135 (MQ=255)
                                                                      aagcACGGCAACCGCGGCAAGCCGGTTTATTTGCATGGCGTTTGGGGCTGCAGCGCAGCCCAAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAg                                                               <  4:119566/134‑1 (MQ=255)
                                                                       agcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGc                                                              >  4:123898/1‑135 (MQ=255)
                                                                         cACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGAGGCGGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCg                                                            <  3:183371/135‑1 (MQ=255)
                                                                           cggaacccGCGCCAAGGCGCTGAAATTGACTGGGGGTGGGGGCGGCAGCGCTGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGt                                                          <  6:81746/129‑1 (MQ=255)
                                                                           cGGCAACCGCGGCAAGCCGGTTTAATTGCTTGGCGTTTGGGGCTGCAGCGAAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGt                                                          <  4:295783/135‑1 (MQ=255)
                                                                            ggCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCCACGGCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGtcc                                                        >  7:218438/1‑136 (MQ=255)
                                                                                aCCGGGGCACGCCGCTTTAATGGATGGGCGTTGGGGGGTGCACCGCCGCCCAAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGtcctcc                                                     <  4:90665/135‑1 (MQ=255)
                                                                                     gccacgCCGCTGAAATTGCATGGCGTTTGGGGCGGCCGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCa                                                <  2:263168/130‑1 (MQ=255)
                                                                                            cGCTGTATTTGCTGGGGGGTGGGGCCGGCAGCGCACCCCAAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGc                                        <  8:101843/136‑1 (MQ=255)
                                                                                                  aaTTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATg                                   >  6:80922/1‑135 (MQ=255)
                                                                                                   atggCATGGCGTTGGGGGCGGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGc                                  <  8:111667/132‑1 (MQ=255)
                                                                                                        agggggTTGGGGGCGGACGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCa                              <  5:70899/129‑1 (MQ=255)
                                                                                                            ctttgggggCTGCAGCGCCGCCCAAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTa                        <  8:6346/131‑1 (MQ=255)
                                                                                                            cGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCg                                                            <  5:63559/100‑1 (MQ=255)
                                                                                                            cGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCg                                                            >  6:63559/1‑100 (MQ=255)
                                                                                                              tttGGGGCTGCAGCGCAGCCCCCACGTCTCAGCCGTTTGCCGCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTAt                       >  4:252188/1‑135 (MQ=255)
                                                                                                              tttGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCGATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCCCCCGCCCTAt                       >  1:223223/1‑135 (MQ=255)
                                                                                                                         agcgcagcCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAg                                               <  4:241684/100‑1 (MQ=255)
                                                                                                                         agcgcagcCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAg                                               >  3:241684/1‑100 (MQ=255)
                                                                                                                           cgcccccccCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCAcc          <  1:209128/130‑1 (MQ=255)
                                                                                                                               gcCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCACCAGGc      <  5:80922/135‑1 (MQ=255)
                                                                                                                               cccccAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCACCAGGc      <  1:44148/134‑1 (MQ=255)
                                                                                                                                   aaaCGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCACCAGGCTCGc  <  2:345883/134‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                     
CCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCACCAGGCTCGC  >  ACB122/1670564‑1670829

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