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.5% 199.9 / ‑5.5 61D233N (GAT→AAT) ydeMputative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (28/31);  total (28/33)
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

CTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTA  >  NC_000913/1579999‑1580188
                                                                                               |                                                                                              
cTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGttt                                                                                            >  2:1671947/1‑100 (MQ=255)
     aGGACTCCTGGAAAATACAACTGGTATGCCCATACCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAAc                                                                                       <  2:855367/100‑1 (MQ=255)
       gtcTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAAccc                                                                                     <  2:162655/98‑1 (MQ=255)
        aCTCCTGGAAAATACAACTGGTATAAAAAAAAAAAAAAAAACGGTTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCa                                                                                    <  2:311066/100‑1 (MQ=255)
          tCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCAtt                                                                                  <  1:174346/100‑1 (MQ=255)
           ccTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTg                                                                                 >  2:37093/1‑100 (MQ=255)
            cTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACTGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGc                                                                                >  2:1474569/1‑100 (MQ=255)
            cTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGc                                                                                <  2:1460716/100‑1 (MQ=255)
            cTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGc                                                                                >  1:1319092/1‑100 (MQ=255)
              ggAAAATACAACTGGTATGCCCATTCCACAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCAACATTGTTTTTAACACATTGc                                                                                >  2:887980/1‑98 (MQ=255)
                aaaaTACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATaa                                                                            <  2:1144776/100‑1 (MQ=255)
                   aTACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATaaaaa                                                                         >  1:1413415/1‑100 (MQ=255)
                   aTACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATaaaaa                                                                         >  1:499225/1‑100 (MQ=255)
                       aaCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATg                                                                       <  1:456382/98‑1 (MQ=255)
                        aCTGGTATGCTCATTCCACAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTAACCAACATTGTTTTTAACCCATTGCATAAAAATgg                                                                      >  2:423353/1‑98 (MQ=255)
                        aCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGtt                                                                    <  2:609951/100‑1 (MQ=255)
                        aCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGtt                                                                    >  1:220370/1‑100 (MQ=255)
                         cTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTg                                                                   >  1:1654011/1‑100 (MQ=255)
                         cTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTg                                                                   <  2:2126411/100‑1 (MQ=255)
                           ggTATGCCCATTCCCCAAAAATCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGAc                                                                 <  1:1854557/100‑1 (MQ=255)
                            gTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACa                                                                <  1:2185341/100‑1 (MQ=255)
                            gTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACa                                                                <  1:441022/100‑1 (MQ=255)
                               tGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATaa                                                             >  1:1227210/1‑100 (MQ=255)
                               tGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATaa                                                             >  2:1018087/1‑100 (MQ=255)
                                gCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATaaa                                                            >  2:1649177/1‑100 (MQ=255)
                                 cccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAAc                                                           >  1:2308714/1‑100 (MQ=255)
                                 cccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAAc                                                           >  1:2342441/1‑100 (MQ=255)
                                  ccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACt                                                          <  2:233100/100‑1 (MQ=255)
                                  ccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACt                                                          >  1:775137/1‑100 (MQ=255)
                                    aTTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTg                                                        <  1:1334701/100‑1 (MQ=255)
                                     ttCCTCAAAAAACGGTTGTCAAAGCTTTCCAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGc                                                       <  1:1191053/100‑1 (MQ=255)
                                     ttCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGc                                                       <  2:499222/100‑1 (MQ=255)
                                     taccccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGc                                                       <  1:1058985/98‑1 (MQ=255)
                                       ccccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAGTTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCa                                                     >  1:1003861/1‑100 (MQ=255)
                                       ccccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCa                                                     <  2:2391293/100‑1 (MQ=255)
                                               aaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGt                                             >  2:2366686/1‑100 (MQ=255)
                                                  ggCTGACAAAGCTTTCAAACTGACGTATGAAAATTTCACCCACATTGTTTTTAACACATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCg                                              >  2:2354105/1‑96 (MQ=255)
                                                  ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTggg                                          <  2:1604931/100‑1 (MQ=255)
                                                     tGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAgg                                       <  2:761078/100‑1 (MQ=255)
                                                     tGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAgg                                       <  2:157110/100‑1 (MQ=255)
                                                      gACAAAGCTTTCAAACTGACGGATGAAAATTTCACCTACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGc                                      <  2:1536662/100‑1 (MQ=255)
                                                      gACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGc                                      <  1:884988/100‑1 (MQ=255)
                                                      gACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGc                                      >  1:1120806/1‑100 (MQ=255)
                                                       aCAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGca                                     >  1:240725/1‑100 (MQ=255)
                                                               ttCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGaaaaa                             <  1:181535/100‑1 (MQ=255)
                                                               ttCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGaaaaa                             <  2:72810/100‑1 (MQ=255)
                                                                  aaatcaaCGGATGACAATTTCACCTACATGGTTTTTAACCCATTTCATAAAAATAGACGACATAATCTTGCCATACGCCGTGGTAGGCACAGAAAAATc                           <  2:944131/93‑1 (MQ=255)
                                                                    aCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                        >  1:1864028/1‑100 (MQ=255)
                                                                     cTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATg                       <  2:246972/100‑1 (MQ=255)
                                                                     cTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATg                       <  1:1391141/100‑1 (MQ=255)
                                                                      tGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGCGGCACAGAAAAATCAATGa                      >  1:2235550/1‑100 (MQ=255)
                                                                      tGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATc                           <  1:733947/95‑1 (MQ=255)
                                                                          ggATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAACCGTGGGGGGCACAGAAAAATCAATGATACg                  >  2:262417/1‑100 (MQ=255)
                                                                                 aaTTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAtgtg           <  1:1483734/100‑1 (MQ=255)
                                                                                 aaTTTCACCCACATCTTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGTCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAtgtg           <  2:1110698/100‑1 (MQ=255)
                                                                                       aCCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCa     >  1:477134/1‑100 (MQ=255)
                                                                                        cccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCAc    >  1:1512100/1‑100 (MQ=255)
                                                                                        cccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCAc    >  1:1919861/1‑100 (MQ=255)
                                                                                         ccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACt   <  2:2031200/100‑1 (MQ=255)
                                                                                         ccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGAGTTCGCACt   <  2:877985/100‑1 (MQ=255)
                                                                                          cacaTTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTa  >  1:2161558/1‑100 (MQ=255)
                                                                                               |                                                                                              
CTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTA  >  NC_000913/1579999‑1580188

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