Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 3,004,952 T→G 27.1% M319R (ATG→AGG)  bcsG → Cellulose biosynthesis protein BcsG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24403,004,9520TG27.1% 57.3 / 2.3 49M319R (ATG→AGG) bcsGCellulose biosynthesis protein BcsG
Reads supporting (aligned to +/- strand):  ref base T (25/8);  new base G (2/13);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 8.42e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCCCGGAGCAGTGCCTGTTGTTCCAGAACCTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATGCCGCAACCGGG‑CATGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAAC  >  KT2440/3004818‑3005085
                                                                                                                                        |                                                                                                                                     
ggCCCCGGAGCAGTGCCTGTTGTTCCAGAACCTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCCTGCCGCAACCGGG‑CCTg                                                                                                                                      >  7:193131/1‑136 (MQ=255)
ggCCCCGGAGCAGTGCCTGTTGTTCCAGAACCTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATGGCGCCCCCGGG‑CAt                                                                                                                                       >  5:381193/1‑135 (MQ=255)
ggCCCCGGAGCAGTGCCTGTTGTTCCAGAACCTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATGCCGCCACCGGG‑CATg                                                                                                                                      >  7:43119/1‑136 (MQ=255)
ggCCCCGGAGCAGTGCCTGTTGTTCCAGAACCTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATGCCGCAACCGGG‑CAt                                                                                                                                       >  4:362458/1‑135 (MQ=255)
ggCCCCGGAGCAGTGCCTGTTGTTCCAGAACCTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATGCCGCAACCGGG‑CAt                                                                                                                                       >  4:308112/1‑135 (MQ=255)
ggCCCCGGAGCAGTGCCTGTTGTTCCAGAACCTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATGCCGCAACCGGG‑CATg                                                                                                                                      >  8:319063/1‑136 (MQ=255)
ggCCCCGGAGCAGTGCCTGTTGTTCCAGAACCTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACACAGCAGCGCATGCCGCAACCGGG‑CAt                                                                                                                                       >  3:189612/1‑135 (MQ=255)
     cGGAGCAGTGCCTGTTGTTCCAGAACCTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCCCAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATGCCGCACCCCGG‑CATGCGCa                                                                                                                                  >  8:326619/1‑135 (MQ=255)
     cGGAGCAGTGCCTGTTGTTCCAGAACCTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCCTGCCGCAACCGGG‑CATGCGCaa                                                                                                                                 >  7:64669/1‑136 (MQ=255)
     cGGAGCAGTGCCTGTTGTTCCAGAACCTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCCTGCCGCAACCCGG‑CATGCGCa                                                                                                                                  >  7:372676/1‑135 (MQ=255)
     cGGAGCAGTGCCTGTTGTTCCAGAACCTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATGCCGCACCCGGG‑CAAgcgc                                                                                                                                   >  2:93391/1‑134 (MQ=255)
                cTGTTGTTCCAGAACCTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGCCAACTTCATCGGCGACATTCCCCAGCAGCGCATGCCGCACCCGGG‑CATGCGCAATACCGACttg                                                                                                                       >  5:52531/1‑134 (MQ=255)
                cTGTTGTTCCAGAACCTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCCTGCCGCCCCCCGG‑CCTGCGCCAAACCgccgtc                                                                                                                       >  6:381982/1‑130 (MQ=255)
                cTGTTGTTCCAGAACCTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCAAGCCGCCACCCGG‑CAGGCGCAATACCGCCGTcc                                                                                                                      >  7:184756/1‑136 (MQ=255)
                     gttCCAGAACCTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATTCCGCGACCCGG‑CCTCGGCGATAACGCCTTcccccc                                                                                                                  >  2:207916/1‑132 (MQ=255)
                         cAGAACAGGGCGCGGGGGGGTTT‑CCATACCAACACGTTGCTAAACCACAGCGGCGGGTTCTACAACTTAATGGCCTACATAACCCCCCCGCGCATGCCGCAACCGGG‑AAGGCGCAATACCGACTTCCCACGGGCAc                                                                                                             <  8:24281/136‑1 (MQ=255)
                          aGAACCTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATGCCGCACCCGGG‑CATGCGCAATACCGCCTTCCCACGGGCAc                                                                                                             >  6:158498/1‑135 (MQ=255)
                           agaCCCGGCCCGGGGGGGTTTGCCAGCCCACCACCTTGCCCAACCACCGCGGCCGGTTCGACAACTTCATCGGCGACATAACCCCGCCGCGCAGGCCGCACCCGGG‑CAGGCGCAATACCGACTTCCCACGGGCAc                                                                                                             <  6:386743/133‑1 (MQ=255)
                              ccTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATGCCGCCACCGGG‑CAGGCGCAATACCGACTTCCCCCGGGCACtggtg                                                                                                        >  7:274477/1‑136 (MQ=255)
                                   cccGGCTGGGTTT‑CCAGAGCAACACCTTGCTCACCCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATGCCGCAACCGGG‑CCTGCGCAATACCGACTTCCCACGGGCACTGGTGGGCt                                                                                                    >  3:405717/1‑135 (MQ=255)
                                              tt‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACACAGCAGCGCATGCCGCAACCGGG‑CATGCGCAATACCGACTTCCCCCGGGCCCTGGTGGGCTTCGACGGCTCGc                                                                                        >  7:252812/1‑136 (MQ=255)
                                               t‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATGCCGCAACCGGG‑CCTGCGCAATACCGACTTCCCCCGGGCACTGGTGGGCTTCGACGGCTCGc                                                                                        >  1:73221/1‑135 (MQ=255)
                                                 ccaccccAACACCTGGCGCAACCCCCCCGGCCGGTTCCACAACTGCACCGGCCACATACCCCAGCAGCGCAGCCCCCACCCGGG‑CAGGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGcc                                                                                       <  4:91127/129‑1 (MQ=255)
                                                        aaCAGCTTGCACACCCACCGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCCGCAGCGCATGCCGCAACCGGG‑CATGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGcc                                                                                <  1:266931/135‑1 (MQ=255)
                                                                 caaaaCCCCCGCGGCCGGTTCTACAACTACATGGGCCACATACCCCAGCCGCGCACCCCCCACCCGGG‑CTGGCCCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGAcc                                                                         <  1:58186/130‑1 (MQ=255)
                                                                  tCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATTCCGCCACCGGG‑CATGCGCAATACCGACTTCCCACGGGCCCCGGTGGGCTTCGACGGCTCGCCGCTTGCCAGTGCCCTgg                                                                      >  2:416463/1‑135 (MQ=255)
                                                                  tCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATGCCGCAACCGGG‑CATGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTgg                                                                      >  4:53063/1‑135 (MQ=255)
                                                                  tCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATGCCGCAACCGGG‑CATGCGCAATAACGACTTCCCACGGGCACTGGGGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGa                                                                     >  8:283992/1‑136 (MQ=255)
                                                                    aaCCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATGCCGCAACCGGG‑CCTTCGCAATACCGACTTCCCACGGGCACTGGTGGGCGTCGACGGCTCGCCGATTGCCAGTGCCCTGGTTg                                                                   >  8:191033/1‑136 (MQ=255)
                                                                       caccGCGGCGGGTTCGACAACTTAATCGGCGACTTTCCCCACCAGCGCATGCCCCAACCGGG‑CAGGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTc                                                                 <  6:188734/135‑1 (MQ=255)
                                                                           ggggCGGGTTCTACAACTTCATCGGCGACATAACCCCGCCGCGCAGGCCCCAACCGGC‑CAGGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTgcg                                                            <  7:144877/134‑1 (MQ=255)
                                                                           gCGGCCGGTTGGAAAAATTCATCGGCTAAATAACCCAGCAGCGCATGCCGCAACCGGG‑CATGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTgcg                                                            <  8:331431/136‑1 (MQ=255)
                                                                            cGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATGCCGCAACCGGG‑CATGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTgcgc                                                           >  8:189855/1‑136 (MQ=255)
                                                                            cGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATGCCGCAACCGGG‑CATGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTgcg                                                            >  4:338332/1‑135 (MQ=255)
                                                                                 ggtgccaCAACTCCAGGGGCCACATACCCCCCCCGCGCAGCCCGCACCCGGG‑CAGGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGtt                                                       <  2:311140/129‑1 (MQ=255)
                                                                                 ggTTCGACAACTTCATCGGCGACATTACACAGCAGCGCATGCCGCAACCGGG‑CATGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGCCCTGGATGTCTTGCGCCGTtg                                                      >  7:280695/1‑136 (MQ=255)
                                                                                        cAACTTCATCGCCGACATTACCAAGCAGCGCATGCCGCAACGGTG‑CATGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTggggg                                               <  8:280695/136‑1 (MQ=255)
                                                                                                tCGGCGACATTACCCAGAAGCGCATGCCGCAACCGGG‑CATGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGGCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCACCGc                                         >  4:265486/1‑134 (MQ=255)
                                                                                                             ccAGCCGCGCAGCCCGCAACGGGG‑CATGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGa                           <  1:50185/135‑1 (MQ=255)
                                                                                                             ccAGCAGCGCATGCCGCAACCGGG‑CATGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGa                           <  2:30591/135‑1 (MQ=255)
                                                                                                               agcagcGCATGCGCCACCGGGG‑CAGGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGa                        <  8:43119/136‑1 (MQ=255)
                                                                                                                    cgcATGCCGCAACCGGG‑CATGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGCCGAACAt                    >  4:393893/1‑135 (MQ=255)
                                                                                                                    cgcATGCCGCAACCGCG‑CATGCGCAATACCGACTTCCCACGGGCACAGGTGGGCTTCGACGGCTAGCCGAATGCCAGTGACCTGGATGTCTAGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGCCCCAACAtg                   >  7:124750/1‑136 (MQ=255)
                                                                                                                    cgcATCCCGCAACCGGC‑CATGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACAtg                   <  8:274477/136‑1 (MQ=255)
                                                                                                                     gcATGCGGCAACCGGG‑CAGGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCAGGACGAACAtg                   <  4:108767/135‑1 (MQ=255)
                                                                                                                       atcccccAACGGGC‑CAGGCCCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACAtgtg                 <  2:173645/129‑1 (MQ=255)
                                                                                                                       aggCCGCAACCGGG‑CATGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATgtgt                <  7:189855/134‑1 (MQ=255)
                                                                                                                       agcccGCAACGGGG‑CAGGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATgtgt                <  8:287240/133‑1 (MQ=255)
                                                                                                                        ggCCGCAACCGGG‑CAGGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATgtgt                <  5:106659/134‑1 (MQ=255)
                                                                                                                         ccccccaCCGGGCCAGGCCCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATgtgt                <  1:416457/129‑1 (MQ=255)
                                                                                                                                  ggg‑CTTGCGCAATACCGACTTCCCACGGGCTCTGGTGGGCTTCGACGGCTCGCCGATTTCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAATGCACCCCGGACGAACATGTGTCACTGTTCtt      <  2:335202/135‑2 (MQ=255)
                                                                                                                                       ttGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTacaac  <  6:52531/134‑1 (MQ=255)
                                                                                                                                       aTGCGCAATACCGACTTCCCACGGGCACTGGTGTGCTTCGCCGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGAACCCCGGACGAACATGTGTCACTGTTCTaca    >  2:389638/1‑133 (MQ=255)
                                                                                                                                       aTGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTaca    >  4:1263/1‑133 (MQ=255)
                                                                                                                                        |                                                                                                                                     
GGCCCCGGAGCAGTGCCTGTTGTTCCAGAACCTGGCCCGGCTGGGTTT‑CCAGAGCAACACCTTGCTCAACCACAGCGGCCGGTTCGACAACTTCATCGGCGACATTACCCAGCAGCGCATGCCGCAACCGGG‑CATGCGCAATACCGACTTCCCACGGGCACTGGTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAAC  >  KT2440/3004818‑3005085

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