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,9390TC97.8% 152.7 / ‑4.0 46D238G (GAC→GGC) yegXputative family 25 glycosyl hydrolase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (21/24);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

TGCCGGGGAGTCTGGTGGGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTT  >  NC_000913/2182794‑2183082
                                                                                                                                                 |                                                                                                                                               
tGCCGGGGAGTCTGGTGGGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGa                                                                                                                                             >  1:182071/1‑150 (MQ=255)
  ccGGGGAGTCTGGTGGGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAAtg                                                                                                                                          <  2:54534/151‑1 (MQ=255)
   cGGGGAGTCTGGTGGGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAAtgc                                                                                                                                         <  2:652098/151‑1 (MQ=255)
             ggtggGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAg                                                                                                                                <  2:633047/150‑1 (MQ=255)
             agagggAGTATATCACTGCTTCTTATTTAAGGCGTTTCCTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAgc                                                                                                                               >  1:502189/4‑151 (MQ=255)
              gtggGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAgc                                                                                                                               >  1:134691/1‑150 (MQ=255)
                 ggAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGTCCGGAATGCTGCCAGAAGCGcc                                                                                                                            >  1:538931/1‑150 (MQ=255)
                  gAGTATATAACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAg                                                                                                                          <  1:255754/151‑1 (MQ=255)
                           aCTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCg                                                                                                                 >  1:183119/1‑151 (MQ=255)
                                    tATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGt                                                                                                          <  1:647071/149‑1 (MQ=255)
                                         aaGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGgacg                                                                                                   <  1:76053/151‑1 (MQ=255)
                                         aaGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGgacg                                                                                                   <  2:435219/151‑1 (MQ=255)
                                            gCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGacgacg                                                                                                >  2:486382/1‑151 (MQ=255)
                                                  cTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGAt                                                                                           >  2:496116/1‑150 (MQ=255)
                                                     tAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGAtagt                                                                                        >  1:183348/1‑150 (MQ=255)
                                                       aTCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATagtagt                                                                                     >  2:594496/1‑151 (MQ=255)
                                                       aTCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATagtagt                                                                                     <  1:158192/151‑1 (MQ=255)
                                                               aaCGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGccacc                                                                             >  2:658865/1‑151 (MQ=255)
                                                               aaCGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGccacc                                                                             <  2:302839/151‑1 (MQ=255)
                                                                 cGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGccaccc                                                                            >  2:2214/1‑150 (MQ=255)
                                                                  gAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGccacccac                                                                          <  1:479585/151‑1 (MQ=255)
                                                                       ccTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACgg                                                                     >  1:501271/1‑151 (MQ=255)
                                                                               tcctccACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCAt                                                             >  2:504003/1‑151 (MQ=255)
                                                                                cctccACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCAt                                                             <  1:95016/150‑1 (MQ=255)
                                                                                 ctccACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTg                                                           >  2:66303/1‑151 (MQ=255)
                                                                                    cACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGaa                                                         >  1:254571/1‑150 (MQ=255)
                                                                                      ccGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAAt                                                       <  2:709096/150‑1 (MQ=255)
                                                                                           ccATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCg                                                 <  1:324903/151‑1 (MQ=255)
                                                                                                  aTACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGAtt                                          <  1:286899/151‑1 (MQ=255)
                                                                                                          aaaTCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGAt                                  >  2:576384/1‑151 (MQ=255)
                                                                                                           aaTCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGAt                                  >  1:444481/1‑150 (MQ=255)
                                                                                                           aaTCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATa                                 <  2:148527/151‑1 (MQ=255)
                                                                                                           aaTCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATa                                 <  2:537729/151‑1 (MQ=255)
                                                                                                              ccACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATaaaa                              <  2:66873/151‑1 (MQ=255)
                                                                                                              ccACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATaaaa                              <  1:349547/151‑1 (MQ=255)
                                                                                                                 ccGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAAc                           <  2:532198/151‑1 (MQ=255)
                                                                                                                 ccGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAAc                           >  2:477340/1‑151 (MQ=255)
                                                                                                                 ccGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAAc                           <  2:137033/151‑1 (MQ=255)
                                                                                                                    gACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGc                        <  1:642350/151‑1 (MQ=255)
                                                                                                                     aCCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGc                        >  2:83017/1‑150 (MQ=255)
                                                                                                                      ccATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTc                      >  1:742920/1‑151 (MQ=255)
                                                                                                                       cATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTcc                     >  1:437910/1‑151 (MQ=255)
                                                                                                                              gCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTaa              <  2:76475/151‑1 (MQ=255)
                                                                                                                                cATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAaat            >  2:576596/1‑151 (MQ=255)
                                                                                                                                     aCCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAaata           <  1:753572/147‑1 (MQ=255)
                                                                                                                                           ccACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCttt  <  2:584655/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TGCCGGGGAGTCTGGTGGGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTT  >  NC_000913/2182794‑2183082

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