Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,474,138 (A)6→7 coding (15/594 nt) bdcR → transcriptional repressor for divergent bdcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,474,1321.A98.5% 214.2 / ‑3.9 67coding (9/594 nt)bdcRtranscriptional repressor for divergent bdcA
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (34/32);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

ACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGT  >  NC_000913/4474006‑4474264
                                                                                                                               |                                                                                                                                     
aCGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAGAAAACAAtc                                                                                                                           >  1:1108538/1‑139 (MQ=255)
  gACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAAtctc                                                                                                                         <  2:1435001/139‑1 (MQ=255)
  gACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAAtctc                                                                                                                         <  1:1170744/139‑1 (MQ=255)
       ccACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGtt                                                                                                                    >  2:1598957/1‑139 (MQ=255)
       ccACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGtt                                                                                                                    >  1:713387/1‑139 (MQ=255)
           cGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAg                                                                                                                <  1:1338579/139‑1 (MQ=255)
             aGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGgt                                                                                                              <  2:1108538/139‑1 (MQ=255)
                           cTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                                <  2:713387/139‑1 (MQ=255)
                              aCCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGc                                                                                             >  1:603050/1‑139 (MQ=255)
                                   tAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTg                                                                                        <  1:1598957/139‑1 (MQ=255)
                                     aaaGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      <  1:237725/139‑1 (MQ=255)
                                             cATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      >  2:974361/1‑131 (MQ=255)
                                             cATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      <  1:974361/131‑1 (MQ=255)
                                                      tCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATTGTCACTAAAAAAACAATCTcg                                                                                                                        <  2:1041125/88‑1 (MQ=255)
                                                      tCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTcg                                                                                                                        >  1:1041125/1‑88 (MQ=255)
                                                      tCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGc                                                                     >  2:1468367/1‑139 (MQ=255)
                                                      tCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGc                                                                     >  2:624663/1‑139 (MQ=255)
                                                        cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAgg                                                                                                               >  1:757313/1‑95 (MQ=255)
                                                        cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAgg                                                                                                               <  2:757313/95‑1 (MQ=255)
                                                        cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTg                                                                                        <  2:1668712/118‑1 (MQ=255)
                                                        cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTg                                                                                        >  1:1668712/1‑118 (MQ=255)
                                                              ttattaTGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCg                                                                       >  1:65006/1‑129 (MQ=255)
                                                              ttattaTGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCg                                                                       <  2:65006/129‑1 (MQ=255)
                                                                attatGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCtt                                                           >  1:1345456/1‑139 (MQ=255)
                                                                   atGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCttttt                                                        >  2:326055/1‑139 (MQ=255)
                                                                          aTGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaa                                                 <  2:171073/139‑1 (MQ=255)
                                                                           tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGgt                                                                                                              <  2:785940/77‑1 (MQ=255)
                                                                           tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGgt                                                                                                              >  1:785940/1‑77 (MQ=255)
                                                                           tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAGAGTGCTTTTTCACCaaaa                                                >  1:1167482/1‑139 (MQ=255)
                                                                           tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCttt                                                          >  2:1562865/1‑129 (MQ=255)
                                                                           tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCttt                                                          <  1:1562865/129‑1 (MQ=255)
                                                                           tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                                >  1:840910/1‑139 (MQ=255)
                                                                           tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                                >  1:1427860/1‑139 (MQ=255)
                                                                              aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTcg                                                                                                                        <  1:943048/64‑1 (MQ=255)
                                                                              aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTcg                                                                                                                        >  2:943048/1‑64 (MQ=255)
                                                                                tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGtt                                           <  2:1345456/139‑1 (MQ=255)
                                                                                  cAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGc                                                                        >  2:1712759/1‑108 (MQ=255)
                                                                                  cAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGc                                                                        <  1:1712759/108‑1 (MQ=255)
                                                                                     aaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      >  2:54704/1‑91 (MQ=255)
                                                                                     aaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      <  1:54704/91‑1 (MQ=255)
                                                                                       cTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTc                                                                                           <  1:1276957/84‑1 (MQ=255)
                                                                                       cTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAAGCCCGCGTTCCAGGTCGTCCCAGACGTTTCGCTc                                                                                           >  2:1276957/1‑84 (MQ=255)
                                                                                          ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGt                                 >  1:1636107/1‑139 (MQ=255)
                                                                                          ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGt                                 <  2:1427860/139‑1 (MQ=255)
                                                                                               aCTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGgt                                                                                                              <  1:173657/57‑1 (MQ=255)
                                                                                               aCTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGgt                                                                                                              >  2:173657/1‑57 (MQ=255)
                                                                                               aCTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAgtg                            >  2:309948/1‑139 (MQ=255)
                                                                                                    cTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaa                                                 <  2:147303/113‑1 (MQ=255)
                                                                                                    cTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaa                                                 >  1:147303/1‑113 (MQ=255)
                                                                                                    cTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCt                       >  1:129544/1‑139 (MQ=255)
                                                                                                         aaCAGGTTTTGTGATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGa                                                                                                     <  1:759154/56‑1 (MQ=255)
                                                                                                         aaCAGGTTTTGTGATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGa                                                                                                     >  2:759154/1‑56 (MQ=255)
                                                                                                         aaCAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCa                                                                                >  1:1039639/1‑77 (MQ=255)
                                                                                                         aaCAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCa                                                                                <  2:1039639/77‑1 (MQ=255)
                                                                                                           cAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTa                <  2:129544/139‑1 (MQ=255)
                                                                                                              gTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGc                        <  1:1020848/128‑1 (MQ=255)
                                                                                                              gTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTCTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGc                        >  2:1020848/1‑128 (MQ=255)
                                                                                                                    tAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGa            <  2:563824/134‑1 (MQ=255)
                                                                                                                    tAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGa            >  1:563824/1‑134 (MQ=255)
                                                                                                                     aaTGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCt      <  2:603050/139‑1 (MQ=255)
                                                                                                                        ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                                <  2:226609/46‑1 (MQ=255)
                                                                                                                        ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                                >  1:226609/1‑46 (MQ=255)
                                                                                                                        ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAg                              >  2:1690782/1‑112 (MQ=255)
                                                                                                                        ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAg                              <  1:1690782/112‑1 (MQ=255)
                                                                                                                        ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGa            <  2:170476/130‑1 (MQ=255)
                                                                                                                        ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGa            >  1:170476/1‑130 (MQ=255)
                                                                                                                         gTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGt  >  2:824595/1‑139 (MQ=255)
                                                                                                                               |                                                                                                                                     
ACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGT  >  NC_000913/4474006‑4474264

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