Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,182,939 T→C D238G (GAC→GGC)  yegX ← putative family 25 glycosyl hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,182,9390TC98.2% 178.9 / NA 57D238G (GAC→GGC) yegXputative family 25 glycosyl hydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (23/33);  minor base A (1/0);  total (24/33)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGCCGGGGAGTCTGGTGGGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTT  >  NC_000913/2182792‑2183084
                                                                                                                                                   |                                                                                                                                                 
ccTGCCGGGGAGTCTGGTGGGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCg                                                                                                                                                 <  1:798269/150‑1 (MQ=255)
   gCCGGGGAGTCTGGTGGGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGa                                                                                                                                               <  2:151461/149‑1 (MQ=255)
    ccGGGGAGTCTGGTGGGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAAtg                                                                                                                                            <  1:530177/151‑1 (MQ=255)
              tggttgaagaagatCACTGCTTCTTAGTAAAGGCGTTTCATTTATCCCATCACCGAACGCCCGCAGCTCCTCCACCGTGCCATTAAATTCATTAAAATCCTCCCGACCATTAATCCCATCTGCCTCTCCACTGCCGGAATGTTGTCat                                                                                                                                     <  1:462261/136‑2 (MQ=255)
              tggtggGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGACGGAATGCTGCCAGa                                                                                                                                    >  1:89847/1‑149 (MQ=255)
                       tataTCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGcc                                                                                                                         >  2:584443/1‑151 (MQ=255)
                          atCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCAt                                                                                                                       <  1:323275/150‑1 (MQ=255)
                                gCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGt                                                                                                                  <  2:414349/149‑1 (MQ=255)
                                gCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCa                                                                                                                <  1:365298/151‑1 (MQ=255)
                                        tttAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCg                                                                                                         >  2:267260/1‑150 (MQ=255)
                                         ttAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCgg                                                                                                        >  1:187541/1‑150 (MQ=255)
                                          tAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGcc                                                                                                                         >  2:610514/1‑132 (MQ=255)
                                          tAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGcc                                                                                                                         <  1:610474/132‑1 (MQ=255)
                                              gCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGacgacg                                                                                                  >  1:180615/1‑151 (MQ=255)
                                                    cTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGAt                                                                                             >  2:628263/1‑150 (MQ=255)
                                                      ttAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGAtag                                                                                           <  1:48004/150‑1 (MQ=255)
                                                      ttAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGAt                                                                                             <  2:542976/148‑1 (MQ=255)
                                                      ttAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGAt                                                                                             <  1:826365/148‑1 (MQ=255)
                                                        aaTCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGAtagtag                                                                                        <  1:781913/151‑1 (MQ=255)
                                                          tCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTg                                                                                      >  1:126644/1‑151 (MQ=255)
                                                           cccATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTg                                                                                      <  2:849580/150‑1 (MQ=255)
                                                            ccATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTg                                                                                      >  2:532318/1‑149 (MQ=255)
                                                                cAACGAATGCCTGCAACTCCTCCCCCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGcca                                                                                 <  1:640879/150‑1 (MQ=255)
                                                                 aaCGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGccacc                                                                               >  1:508005/1‑151 (MQ=255)
                                                                    gAATGCCTGCAACTCCTCCACCGTGCCATTAAATTCATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGccacccac                                                                            >  2:323845/1‑151 (MQ=255)
                                                                    gAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGccacccac                                                                            >  2:359479/1‑151 (MQ=255)
                                                                              aaCTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATAtt                                                                  <  2:88809/151‑1 (MQ=255)
                                                                                 tcctccACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCAt                                                               >  2:114021/1‑151 (MQ=255)
                                                                                  cctccACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCa                                                                <  2:638657/149‑1 (MQ=255)
                                                                                      cACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTg                                                             >  1:3792/1‑148 (MQ=255)
                                                                                       aCCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGAta                                                                                            >  2:540230/1‑116 (MQ=255)
                                                                                       aCCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGAta                                                                                            <  1:540193/116‑1 (MQ=255)
                                                                                       aCCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAAt                                                         >  1:172160/1‑151 (MQ=255)
                                                                                          gTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGc                                                      >  2:345404/1‑151 (MQ=255)
                                                                                                 tAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCa                                                <  2:688151/150‑1 (MQ=255)
                                                                                                  aaaTACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAg                                               <  2:726242/150‑1 (MQ=255)
                                                                                                   aaTACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATATCCCGcc                                                 >  2:277315/1‑147 (MQ=255)
                                                                                                    aTACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGAtt                                            <  1:779638/151‑1 (MQ=255)
                                                                                                          aaaaaTCCACCGGACCATTGATGACATCTACCTGTCCACGGCTGACATCCTGCCAGAATCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGAGAGCCACCCCCCACGGATATTCATTGAAATAGCCCGCCAGCTTGGtgtg                                      >  1:265695/2‑151 (MQ=255)
                                                                                                             aaTCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATa                                   <  2:805506/151‑1 (MQ=255)
                                                                                                                ccACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATaaaa                                <  1:616834/151‑1 (MQ=255)
                                                                                                                 cACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACTTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATaaaaa                               <  1:687727/151‑1 (MQ=255)
                                                                                                                 cACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATaaaaa                               <  2:136863/151‑1 (MQ=255)
                                                                                                                   ccGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCTCCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAAc                             <  2:615737/151‑1 (MQ=255)
                                                                                                                   ccGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAAc                             <  1:767432/151‑1 (MQ=255)
                                                                                                                     ggACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACgg                           <  2:665586/151‑1 (MQ=255)
                                                                                                                            tGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTg                     >  1:381867/1‑150 (MQ=255)
                                                                                                                             gATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGa                    <  1:525194/150‑1 (MQ=255)
                                                                                                                              aTGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGt                  >  1:259619/1‑151 (MQ=255)
                                                                                                                                 ccATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTaa                >  2:779773/1‑150 (MQ=255)
                                                                                                                                   aTCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAaat              >  2:475196/1‑150 (MQ=255)
                                                                                                                                    tCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAaata             >  2:492975/1‑150 (MQ=255)
                                                                                                                                       aCCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAaataat           <  2:244702/149‑1 (MQ=255)
                                                                                                                                        ccTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCg         <  1:45548/150‑1 (MQ=255)
                                                                                                                                             ccACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCtttt   <  1:703264/151‑1 (MQ=255)
                                                                                                                                             ccACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCttt    <  2:73780/150‑1 (MQ=255)
                                                                                                                                              cACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCttttt  <  1:226698/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCTGCCGGGGAGTCTGGTGGGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTTT  >  NC_000913/2182792‑2183084

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