Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,580,094 C→T 100% D233N (GAT→AAT)  ydeM ← putative YdeN‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,580,0940CT98.9% 243.5 / ‑5.8 74D233N (GAT→AAT) ydeMputative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (36/37);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTAT  >  NC_000913/1579997‑1580189
                                                                                                 |                                                                                               
tCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGt                                                                                               >  2:1443379/1‑100 (MQ=255)
 ccTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGtt                                                                                              >  1:2974661/1‑100 (MQ=255)
  cTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGttt                                                                                             <  1:203451/100‑1 (MQ=255)
  cTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGttt                                                                                             <  2:542363/100‑1 (MQ=255)
   ttGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGtttt                                                                                            <  1:2105197/100‑1 (MQ=255)
   ttGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGtttt                                                                                            <  2:2313673/100‑1 (MQ=255)
   ttGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGtttt                                                                                            <  1:733417/100‑1 (MQ=255)
    tggatgtctgcTGGAAAATACAACTGGTATGCCGATTCTCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGCATGAATATTTCACCCACATTGttttt                                                                                           <  2:2337976/90‑1 (MQ=255)
       aGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAAc                                                                                        <  1:2089806/100‑1 (MQ=255)
           cTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCAt                                                                                    <  2:2005749/100‑1 (MQ=255)
            tCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCAtt                                                                                   <  1:1392811/100‑1 (MQ=255)
            tCCTGGAAAATACAACTAGTATGCCCATTTACCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCATCCACATTGTTTTTAACCCa                                                                                     <  2:563574/98‑1 (MQ=255)
             ccTGGAAAATACAACTGGTATGTCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTg                                                                                  >  1:912877/1‑100 (MQ=255)
              cTGTAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTg                                                                                  <  2:2669234/99‑1 (MQ=255)
              cTGGGAAATACAACAGGTATGCCCATGCCCCAAAAAACGGCTGACAACGCTGTCAAACTGACTGATAAAAAAATCACCCACATCGTTTTTAACCCttt                                                                                   >  2:241380/1‑98 (MQ=255)
              cTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGc                                                                                 >  1:1221639/1‑100 (MQ=255)
               tGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCa                                                                                >  1:2607452/1‑100 (MQ=255)
                ggAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCAt                                                                               >  1:197111/1‑100 (MQ=255)
                    aaTACAACTTGTATGCCTATTCCCCAAAAAACGGCTGACAATACTTTCAAACTGACTGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATaaaa                                                                           <  2:1298950/100‑1 (MQ=255)
                       aCAACTGGTATGCCCATTCCCCAAAAAACGTCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATAGCATAAAAATg                                                                        <  2:2839135/100‑1 (MQ=255)
                        cAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACACACATTGTTTTTAACCAATTGAATAAAAATg                                                                        >  1:608831/1‑99 (MQ=255)
                          aCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGtt                                                                     >  1:1905883/1‑100 (MQ=255)
                          aCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGtt                                                                     >  2:1192491/1‑100 (MQ=255)
                          aCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGtt                                                                     >  2:2750718/1‑100 (MQ=255)
                           cTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTg                                                                    <  2:664823/100‑1 (MQ=255)
                           cTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTg                                                                    <  1:1761037/100‑1 (MQ=255)
                           cTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTg                                                                    <  1:1761221/100‑1 (MQ=255)
                             ggTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGAc                                                                  >  1:834138/1‑100 (MQ=255)
                              gTATGCCTATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGAc                                                                  <  1:1347929/99‑1 (MQ=255)
                              gTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACa                                                                 <  1:2192854/100‑1 (MQ=255)
                              gTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACa                                                                 >  1:568467/1‑100 (MQ=255)
                                ttccgatctCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATa                                                               <  1:2471026/92‑1 (MQ=255)
                                  gCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATaaa                                                             >  1:1937080/1‑100 (MQ=255)
                                      aTTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTg                                                         <  1:176022/100‑1 (MQ=255)
                                        tCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGcc                                                       <  1:2795935/100‑1 (MQ=255)
                                         ccccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATaagatcggaa                                                      >  2:2470997/1‑92 (MQ=255)
                                         ccccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCa                                                      <  2:280279/100‑1 (MQ=255)
                                         ccccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCa                                                      >  1:1689074/1‑100 (MQ=255)
                                          cccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCAt                                                     >  1:2844603/1‑100 (MQ=255)
                                          cccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCAt                                                     >  1:1636217/1‑100 (MQ=255)
                                            cAAAAAACTGCTGACAAAGCTTTCAAGCTGACGTATGAAAATTTCACCCACATTGTTTTTAACCCCTTGCATAAACATGGTTGACATAAACGTGCCAt                                                     <  2:2643472/98‑1 (MQ=255)
                                              aaaaaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGc                                                 <  2:1122046/100‑1 (MQ=255)
                                                   cGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTgg                                            >  2:2433711/1‑100 (MQ=255)
                                                    ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAATCTTGCCATAAGCCGTggg                                           <  2:1483362/100‑1 (MQ=255)
                                                    ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTggg                                           <  1:2685639/100‑1 (MQ=255)
                                                    ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTgg                                            >  2:2698660/1‑99 (MQ=255)
                                                    ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTgg                                            <  2:1247445/99‑1 (MQ=255)
                                                    ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACACAAACTTGCCATAAGCCGTgg                                            >  2:2698842/1‑99 (MQ=255)
                                                      cTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAg                                         >  2:2312305/1‑100 (MQ=255)
                                                       tGACAAAGCTTTCAAACTGACGGATGAAAATTTCATCCATATTGTTTTTAATCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAgg                                        <  1:776782/100‑1 (MQ=255)
                                                        gACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGc                                       <  2:453334/100‑1 (MQ=255)
                                                        gACAAAGCTTTCAAAATGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGc                                       >  2:595660/1‑100 (MQ=255)
                                                          cAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGcac                                     >  1:1487174/1‑100 (MQ=255)
                                                           aaaGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGcaca                                    >  2:1050639/1‑100 (MQ=255)
                                                            aaGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGcaaa                                    >  2:1783004/1‑97 (MQ=255)
                                                            aaGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAg                                   <  1:1932492/100‑1 (MQ=255)
                                                            aaGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAg                                   >  2:2176097/1‑100 (MQ=255)
                                                                   cAAACTGACGGATGAAATTTTTATCCACTTTTTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATc                            <  1:465307/100‑1 (MQ=255)
                                                                    aaaCTGACGGATGAAAATTTCACCCACATTGATTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAAt                             >  2:1310279/1‑98 (MQ=255)
                                                                      aCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                         <  1:1172386/100‑1 (MQ=255)
                                                                      aCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                         <  1:2750750/100‑1 (MQ=255)
                                                                      aCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                         <  1:1287981/100‑1 (MQ=255)
                                                                        tGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGa                       >  1:2089613/1‑100 (MQ=255)
                                                                         gACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGAt                      >  1:1974744/1‑100 (MQ=255)
                                                                         gACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGAt                      >  1:1711215/1‑100 (MQ=255)
                                                                         gACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGAt                      >  1:1709166/1‑100 (MQ=255)
                                                                         gACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTACCATAAGCCATAGGAGGCACAGAAAAATCAATGAt                      >  1:1041246/1‑100 (MQ=255)
                                                                             gATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACgg                  >  1:384121/1‑100 (MQ=255)
                                                                             gATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACg                   <  1:2659599/99‑1 (MQ=255)
                                                                                gAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTTCCATAAGCCGTAGGAGGCACAGAAAAATCAATGATACGGAAt               >  1:684433/1‑100 (MQ=255)
                                                                                        cACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtctc       <  2:2197822/100‑1 (MQ=255)
                                                                                           ccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACt    >  1:2934017/1‑100 (MQ=255)
                                                                                             acaTTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTAt  <  1:906949/100‑1 (MQ=255)
                                                                                              caTTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTAt  >  1:904855/1‑99 (MQ=255)
                                                                                                 |                                                                                               
TCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTAT  >  NC_000913/1579997‑1580189

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