Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,580,094 C→T 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,0940CT94.1% 158.9 / 0.9 51D233N (GAT→AAT) ydeMputative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (22/26);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

TGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATG  >  NC_000913/1580001‑1580190
                                                                                             |                                                                                                
tGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGttttt                                                                                            <  1:601278/100‑1 (MQ=255)
 gCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTa                                                                                           >  1:768363/1‑100 (MQ=255)
  cAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTaa                                                                                          >  1:54432/1‑100 (MQ=255)
         ccTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTg                                                                                   >  1:1179177/1‑100 (MQ=255)
          cTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGc                                                                                  <  1:967800/100‑1 (MQ=255)
          cTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGc                                                                                  >  1:1570119/1‑100 (MQ=255)
          cTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGc                                                                                  >  1:1570166/1‑100 (MQ=255)
            ggAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCAt                                                                                >  2:93337/1‑100 (MQ=255)
               aaaTACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATaaa                                                                             >  1:227450/1‑100 (MQ=255)
                 aTACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATaaaaa                                                                           <  1:110500/100‑1 (MQ=255)
                    cAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATgg                                                                        >  1:1009593/1‑100 (MQ=255)
                    cAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATgg                                                                        <  1:1355057/100‑1 (MQ=255)
                    cAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATg                                                                         >  1:608552/1‑99 (MQ=255)
                       cTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTg                                                                     <  2:437183/100‑1 (MQ=255)
                       cTGGTATGCCCATTCCACAAAAAACGGCTGACAAAGCTTTCAAACTGCCGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTg                                                                     <  1:619036/100‑1 (MQ=255)
                       cTGGTATGCCCATTCCACAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTg                                                                     <  2:1222494/100‑1 (MQ=255)
                       cTGGTATGCCCATTCACCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTg                                                                     <  2:1500827/100‑1 (MQ=255)
                         ggTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATGGCATAAAAATGGTTGAc                                                                   >  1:165910/1‑100 (MQ=255)
                           tATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACAt                                                                 <  2:492448/100‑1 (MQ=255)
                                  aTTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTg                                                          <  2:65477/100‑1 (MQ=255)
                                   ttCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGc                                                         <  1:1194141/100‑1 (MQ=255)
                                    tCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGcc                                                        <  1:1395220/100‑1 (MQ=255)
                                           aaaaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGcc                                                 <  1:1355130/100‑1 (MQ=255)
                                             aaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGt                                               <  2:1009587/100‑1 (MQ=255)
                                              aCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCAAAAAAATGGTTGACATAAACTTGCCATAAGCCGTg                                              <  2:195511/100‑1 (MQ=255)
                                                ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTggg                                            <  2:38885/100‑1 (MQ=255)
                                                ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTggg                                            <  2:1179172/100‑1 (MQ=255)
                                                ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTggg                                            <  1:93337/100‑1 (MQ=255)
                                                    gACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGc                                        <  1:1476299/100‑1 (MQ=255)
                                                      cAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGcac                                      <  1:874657/100‑1 (MQ=255)
                                                        aaGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAg                                    >  1:650609/1‑100 (MQ=255)
                                                          gCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGaa                                  <  2:286286/100‑1 (MQ=255)
                                                               cAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATc                             <  1:157054/100‑1 (MQ=255)
                                                                aaaCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCa                            >  1:858507/1‑100 (MQ=255)
                                                                aaaCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCa                            <  2:608550/100‑1 (MQ=255)
                                                                    tGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGa                        >  1:636250/1‑100 (MQ=255)
                                                                         gATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGAACAGAAAAATCAATGATACgg                   >  1:32058/1‑100 (MQ=255)
                                                                          aTGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAGAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGa                  <  2:768361/100‑1 (MQ=255)
                                                                           tGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGaa                 >  2:558231/1‑100 (MQ=255)
                                                                             aaaaTTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAtg               >  1:275733/1‑100 (MQ=255)
                                                                             aaaaTTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAtg               >  1:9747/1‑100 (MQ=255)
                                                                              aaaTTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAtgt              <  2:1082330/100‑1 (MQ=255)
                                                                                  ttCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtc          >  1:214420/1‑100 (MQ=255)
                                                                                    cACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtctc        >  2:1430837/1‑100 (MQ=255)
                                                                                     aCCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCa       <  1:56362/100‑1 (MQ=255)
                                                                                      cccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCAc      >  2:137968/1‑100 (MQ=255)
                                                                                      cccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCAc      >  1:842305/1‑100 (MQ=255)
                                                                                       ccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACt     <  2:435609/100‑1 (MQ=255)
                                                                                       ccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACt     <  2:1236396/100‑1 (MQ=255)
                                                                                        cacaTTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTa    >  2:735067/1‑100 (MQ=255)
                                                                                          caTTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATg  >  2:1293970/1‑100 (MQ=255)
                                                                                             |                                                                                                
TGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATG  >  NC_000913/1580001‑1580190

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