Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,816,318 T→C 8.3% intergenic (+190/+68) katE → / ← chbG catalase II/chitin disaccharide deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,816,3180TC 8.3% 114.3 / 2.4 51intergenic (+190/+68)katE/chbGcatalase II/chitin disaccharide deacetylase
Reads supporting (aligned to +/- strand):  ref base T (24/22);  new base C (2/2);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

TTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAAC  >  NC_000913/1816176‑1816458
                                                                                                                                              |                                                                                                                                            
tttGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGt                                                                                                                                              >  2:209079/1‑143 (MQ=255)
tttGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGt                                                                                                                                              <  1:554776/143‑1 (MQ=255)
tttGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGt                                                                                                                                              <  1:515937/143‑1 (MQ=255)
    gttATGCCGGATGCGGCGTGAACGCCTTATCAGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTATGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTAcg                                                                                                                                          <  1:160350/142‑1 (MQ=255)
         gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTGAATATATTGCAGGAACCACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTTCGTTCGTCAGCCGTCTCAAGCGGCGCAGTTGGgccgcc                                                                                                                                     >  2:101382/1‑143 (MQ=18)
          ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTATACAAATTCAATATATTGCAGGAAACACGTAGACCTGATAAGCGAAGCGCATCAGGCAGATTTGCGTTTGTCAGCAGTCTCAAGCAGCACAATTACGCCCCCt                                                                                                                                    >  2:336333/1‑143 (MQ=18)
            ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCttt                                                                                                                                  <  1:485879/143‑1 (MQ=255)
             gATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTACAAATTCAATATATTGCAGGAAACGCGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTg                                                                                                                                 >  2:530879/1‑143 (MQ=255)
              aTGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGt                                                                                                                                <  2:428333/143‑1 (MQ=255)
               tGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTa                                                                                                                               <  1:36177/143‑1 (MQ=255)
                gcggcgTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTCTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTGGTAg                                                                                                                              >  2:622023/1‑143 (MQ=255)
                 cggcgTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCACGCGGCGCAGTGACGCCGCCTTTGTAgg                                                                                                                             >  1:522826/1‑143 (MQ=255)
                      tGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCCGTTACGCCGCCTTTGTAGGAATTa                                                                                                                        >  1:128377/1‑143 (MQ=255)
                           gCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAGACGAATTAGGCAGTTTTGCGCTAGCGAGCTGTCTCAAGCGGCGCAGTGACGGCGGCGTTGTCGGCGTTGATCGc                                                                                                                   >  2:142692/1‑143 (MQ=18)
                                aTCCGGCCTACAAAACTGTGGAAATTCAATAGATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGAt                                                                                                              <  2:493710/143‑1 (MQ=255)
                                 tCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATg                                                                                                             <  1:207003/143‑1 (MQ=255)
                                     tcttACAAAACTGTGTAAATTTTATATATTGCAGGAAACATGTAGGCCTGATAAGCGAAGCGCATCATGTATTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAg                                                                                                         <  2:280847/140‑1 (MQ=37)
                                       cTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGt                                                                                                       >  1:125464/1‑143 (MQ=255)
                                           aaaaCTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCTGCGCAGCTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCAc                                                                                                   >  1:218107/1‑143 (MQ=255)
                                           aaaaCTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGCTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCAc                                                                                                   >  1:210827/1‑143 (MQ=255)
                                           aaaaCTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCACATCAGGCAGTGTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTGGTAGGAATTAATCGCCGGAGGCAAGGTTCAc                                                                                                   >  2:311997/1‑143 (MQ=255)
                                             tacggtgCAAATTCAATATCTTGCAGGAAACACGGTGGCCGGATATGTGTAGCGCGTCAGGCAGTTTCGCGTGTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTTAcgc                                                                                                 <  2:358097/139‑1 (MQ=255)
                                             aaCTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCAcgc                                                                                                 <  2:383264/143‑1 (MQ=255)
                                             aaCTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCAcgc                                                                                                 <  1:634376/143‑1 (MQ=255)
                                                  tgCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGATACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCAt                                                                                            >  1:20034/1‑143 (MQ=255)
                                                     aaaTTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTg                                                                                         <  1:432885/143‑1 (MQ=255)
                                                       aTTCAATATATTGCAGGAAACACGTGGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGc                                                                                       >  1:475768/1‑143 (MQ=255)
                                                        ttCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCa                                                                                      <  1:440561/143‑1 (MQ=255)
                                                               tatTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTAACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATcc                                                                               >  1:338916/1‑143 (MQ=255)
                                                                  tGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCa                                                                            >  1:366246/1‑143 (MQ=255)
                                                                  tGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGGTATGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGGAAGGTGCACGCCGGATCTGGCCAACATCCTcc                                                                            >  2:190986/1‑142 (MQ=38)
                                                                    cAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTACTCGCCGGATGCAAGGTTCACGCCGCATGTGGCAAACACCCGCACt                                                                          >  1:348066/1‑143 (MQ=255)
                                                                      ggAATCACGTGGGCGGGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTa                                                                        <  2:126055/143‑1 (MQ=255)
                                                                        aaaCACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGCTCACGCCGCATCTGGCAAACATCCTCACTTaca                                                                      >  2:252508/1‑143 (MQ=255)
                                                                          acccGTATGGCTGATAAGCGGAGCGCATCAGGCAGGTGTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTAcaca                                                                    <  1:43559/140‑1 (MQ=255)
                                                                             cGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATcc                                                                 <  1:441129/143‑1 (MQ=255)
                                                                              gTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAGTCGCCGGATGCGAGGTTCACGCCGCATCTGGCGACCACCCCCACTTACCCATccc                                                                >  1:59705/1‑143 (MQ=37)
                                                                                       aTAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGCCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTcc                                                       <  1:41023/143‑1 (MQ=255)
                                                                                           gcgaagcgCATCCGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTACTCGCCGGCTGCCAGGTTCACGCCGCATCTGGCAAACATACTCACTCGCCCATCCCGCTAACGccccca                                                   >  2:88056/1‑141 (MQ=255)
                                                                                           gcgaagcgCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCaa                                                   >  1:73509/1‑143 (MQ=255)
                                                                                              aagcgCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCGCGACAACACCCCAACCg                                                >  2:237949/1‑143 (MQ=255)
                                                                                              aagcgCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGAGGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCg                                                >  1:622015/1‑143 (MQ=255)
                                                                                                gctgaTCAGGCAGTTTCGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGAt                                              <  2:307835/139‑1 (MQ=255)
                                                                                                         gCAGTTTTGCGTTTGTCAGCAGTCGCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCt                                     <  1:343589/143‑1 (MQ=255)
                                                                                                                        tCAGCAGTCTCAAGCGGCGCAGCTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCtt                      <  2:218107/143‑1 (MQ=255)
                                                                                                                        tCAGCAGTCTCAAGCGGCGCAGCTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCtt                      <  2:210827/143‑1 (MQ=255)
                                                                                                                               tctcAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAacga               >  1:359427/1‑143 (MQ=255)
                                                                                                                               tctcAAGCGGCGCAGTTACGCCGCCTTTGTAGGAACTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCCGAGCCATAGCACCTTCCCacga               >  2:208326/1‑143 (MQ=255)
                                                                                                                                 tcAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGGCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAacgacg             <  2:56739/143‑1 (MQ=255)
                                                                                                                                   aaGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCt           >  1:346815/1‑143 (MQ=255)
                                                                                                                                    aGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTg          <  1:381615/143‑1 (MQ=255)
                                                                                                                                        gcgcAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGt      <  2:596845/143‑1 (MQ=255)
                                                                                                                                          gcAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCa    <  1:456088/143‑1 (MQ=255)
                                                                                                                                            aGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAac  >  1:34383/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                            
TTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAAC  >  NC_000913/1816176‑1816458

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