Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 584,262 (A)7→8 coding (582/750 nt) appY → global transcriptional activator; DLP12 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913584,2551.A89.2% 239.9 / 22.1 75coding (575/750 nt)appYglobal transcriptional activator; DLP12 prophage
Reads supporting (aligned to +/- strand):  ref base . (3/5);  new base A (41/25);  total (45/30)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTA  >  NC_000913/584156‑584347
                                                                                                      |                                                                                             
cGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGc                                                                                                >  2:846530/1‑100 (MQ=255)
cGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGc                                                                                                <  1:986738/100‑1 (MQ=255)
 ggAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCa                                                                                              <  2:806518/100‑1 (MQ=255)
     ttGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCaaaaa                                                                                          <  2:1632474/100‑1 (MQ=255)
      tGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCaaaaaa                                                                                        >  2:179403/1‑100 (MQ=255)
       gATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCaaaaaaa                                                                                       >  1:1486692/1‑100 (MQ=255)
        aTTTATACGAGTGAAAGTTTAATAAAAAAAAGCTTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCaaaaaaa                                                                                       >  1:228159/1‑99 (MQ=255)
        aTTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCCAAAAAAc                                                                                      <  2:1025106/100‑1 (MQ=255)
          ttATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCaaaaaaaact                                                                                    <  1:791840/100‑4 (MQ=255)
           tataCGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGAAACGTCATTTACTGAAATATTGAGAGATACTAGGATTAGGTATACaaaaaaaac                                                                                     >  1:189133/1‑96 (MQ=255)
            ataCGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCAt                                                                                  >  2:1764141/1‑100 (MQ=255)
              aCGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATa                                                                                 >  2:501408/1‑100 (MQ=255)
              aCGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATa                                                                                 <  2:1091852/100‑1 (MQ=255)
                gAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAAc                                                                               >  2:369536/1‑100 (MQ=255)
                 aGTGAAAGTTTAATAAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTAAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAAc                                                                               >  1:230420/1‑100 (MQ=255)
                 aGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACt                                                                              >  2:1592712/1‑100 (MQ=255)
                  gTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACtt                                                                             >  2:1257018/1‑100 (MQ=255)
                  gTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACtt                                                                             >  1:1749369/1‑100 (MQ=255)
                  gTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACtt                                                                             >  1:1390392/1‑100 (MQ=255)
                  gTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTc                                                                            <  2:800103/100‑1 (MQ=255)
                   tGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTc                                                                            >  1:793869/1‑100 (MQ=255)
                     aaaaTTTAATAAAAAAAAGATTAAGGGATGAAGGAAAGTCATTTACTGAAATATTGAGAGATACTAGGATGAAAAATAAAAAAAAAACTCATAACTTCaa                                                                          <  1:1592701/100‑1 (MQ=255)
                     aaaGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCa                                                                           >  2:1191614/1‑99 (MQ=255)
                       aGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAAc                                                                        <  2:179708/100‑1 (MQ=255)
                             ataaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTAt                                                                  >  1:1789198/1‑100 (MQ=255)
                             ataaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTAt                                                                  >  2:100306/1‑100 (MQ=255)
                             ataaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAACTCATAACTTCAAACTCTTaat                                                                 >  2:1724636/1‑97 (MQ=255)
                                aaaaaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTc                                                                >  1:1769070/1‑100 (MQ=255)
                                 aaaaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCt                                                               <  1:880753/100‑1 (MQ=255)
                                  aaaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTa                                                              <  2:1033969/100‑1 (MQ=255)
                                  aaaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAATAACTCATAGCTTCAAACTCTTATTc                                                                <  1:850676/98‑1 (MQ=255)
                                   aaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTAt                                                             >  1:1904689/1‑100 (MQ=255)
                                      aaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATAAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCaa                                                          <  2:1485630/100‑1 (MQ=255)
                                        gATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATg                                                        >  1:506822/1‑100 (MQ=255)
                                         aTTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTc                                                      >  2:1776144/1‑100 (MQ=255)
                                          ttAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTc                                                      >  2:1595581/1‑100 (MQ=255)
                                          ttAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTc                                                      >  1:1942630/1‑100 (MQ=255)
                                           tAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCg                                                     <  1:193264/100‑1 (MQ=255)
                                              gggATGAAGGAACGTCATTTACTGGAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAg                                                  >  2:986520/1‑100 (MQ=255)
                                                gATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAgac                                                >  1:1181635/1‑98 (MQ=255)
                                                gATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGcc                                                >  1:617360/1‑100 (MQ=255)
                                                gATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGcc                                                <  2:1670184/100‑1 (MQ=255)
                                                gATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGcc                                                <  2:1579722/100‑1 (MQ=255)
                                                 aTGAAGGAACATCATTTACTTAAATATTGAGATATTCTAGGATGAAAAAAAAAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGccc                                               <  1:407669/100‑1 (MQ=255)
                                                  tGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCa                                              <  2:1594810/100‑1 (MQ=255)
                                                   gAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAg                                             <  2:295198/100‑1 (MQ=255)
                                                    aaGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGa                                            <  1:519867/100‑1 (MQ=255)
                                                       gAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAt                                         >  2:63408/1‑100 (MQ=255)
                                                        aaCGTCATTTACTGAAATATTGAGAGATACTAGGATGAAGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtg                                        <  1:1791536/100‑1 (MQ=255)
                                                         acatCATTTACTGAAATATTGAGAGATACTAGGAAAAAAAAAAAAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGaaa                                          <  1:425863/94‑1 (MQ=255)
                                                          cGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtgtg                                      <  1:716675/100‑1 (MQ=255)
                                                          cGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtgtg                                      >  2:1824990/1‑100 (MQ=255)
                                                          cGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTgg                                     <  1:1055338/100‑1 (MQ=255)
                                                           gTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTgg                                     >  2:1238693/1‑100 (MQ=255)
                                                            tCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGc                                    <  2:1256/100‑1 (MQ=255)
                                                              aTTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCta                                  >  2:1659943/1‑100 (MQ=255)
                                                                 tACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATaa                               >  2:1821843/1‑100 (MQ=255)
                                                                 tACTGAAATATTGAGAGATACTAGGATGAAAAATAAAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCtata                                <  1:216027/99‑1 (MQ=255)
                                                                  aCTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAAc                              <  2:1749380/100‑1 (MQ=255)
                                                                   cTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACa                             >  2:277733/1‑100 (MQ=255)
                                                                   cTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACa                             <  2:91001/100‑1 (MQ=255)
                                                                   cTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCACAGAAATGTGGCTATAACa                             >  1:1058145/1‑100 (MQ=255)
                                                                    tGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAg                            <  2:591199/100‑1 (MQ=255)
                                                                     gAAATATTGAGAGATACTAGGATGAGGTATGCGAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGt                           >  1:564530/1‑100 (MQ=255)
                                                                     gAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGGAATGTGGCTATAACAGt                           >  2:1899066/1‑100 (MQ=255)
                                                                     gAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGt                           <  2:55084/100‑1 (MQ=255)
                                                                     gAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGt                           >  2:1682325/1‑100 (MQ=255)
                                                                     gAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGt                           >  2:178994/1‑100 (MQ=255)
                                                                     gAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCGTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGt                           >  1:1364871/1‑100 (MQ=255)
                                                                         tatTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACtt                       >  1:38072/1‑100 (MQ=255)
                                                                          atTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttc                      >  1:751394/1‑100 (MQ=255)
                                                                          atTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttc                      <  1:1606097/100‑1 (MQ=255)
                                                                           ttGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttca                     >  1:257402/1‑100 (MQ=255)
                                                                             gagagaTACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcata                   >  1:878619/1‑100 (MQ=255)
                                                                             gagagaTACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcata                   >  1:1541166/1‑100 (MQ=255)
                                                                                   tACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcat             >  1:713134/1‑100 (MQ=255)
                                                                                    aCTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcata            >  1:1467569/1‑100 (MQ=255)
                                                                                    aCTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcata            <  2:564766/100‑1 (MQ=255)
                                                                                    aCTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcata            <  2:1876485/100‑1 (MQ=255)
                                                                                     cTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcatat           >  2:1747015/1‑100 (MQ=255)
                                                                                        ggATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATAtgtg        <  1:1075271/100‑1 (MQ=255)
                                                                                            gAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCAtt    >  1:781674/1‑100 (MQ=255)
                                                                                              ggTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTa  >  1:1495276/1‑100 (MQ=255)
                                                                                              ggTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTa  >  2:958330/1‑100 (MQ=255)
                                                                                                      |                                                                                             
CGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTA  >  NC_000913/584156‑584347

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